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ABSTRACT 

Since the advent of SCUBA on the JCMT, a series of complementary surveys has re- 
solved the bulk of the far-infrared extragalactic background into discrete sources. This 
has revealed a population of heavily dust-obscured sources at high redshift (z > 1) 
undergoing an intense period of massive star- forming activity with inferred star forma- 
tion rates of several hundred to several thousand solar masses per year. Taken together, 
these existing surveys cover a total area of 460 sq. arcmin to a range of depths, but 
combining the results has hitherto been complicated by the fact that different survey 
groups have used different methods of data reduction and source extraction. In this pa- 
per we re-reduce and analyse all of the blank field surveys to date in an almost identical 
manner to that employed in the "SCUBA 8 mJy Survey" . Comparative source cata- 
logues are given which include a number of new significant source detections as well as 
failing to confirm some of those objects previously published. These new source cata- 
logues have been combined to produce the most accurate number counts to date from 
2 to 12.5 mJy. We find N(> 4mJy) = 620^19^ deg-^, N(> 6mJy) ^ 310ti|5g deg-^ 

N(> 8mJy) = 150l go deg"^ and N(> lOmJy) = 40l 20 deg"^ after correcting for 
the effects of incompleteness, flux-density boosting and contamination from spurious 
/ confused detections. Furthermore the cumulative number counts appear to steepen 
beyond Sg5o > SmJy, which could indicate an intrinsic turn-over in the underlying 
luminosity function placing an upper limit on the luminosity (and hence mass) of a 
high redshift galaxy. We have also investigated the clustering properties of the bright 
Ssso > 5mJy SCUBA population by means of 2-point angular correlation functions. 
We find a ~ 3.5cr excess of pairs within the first 100 arcsec over that expected from 
a Poisson distribution. Fits of a standard power-law of the form w(6') = K9~^ to the 
angular correlation functions for Ssso > 5 mJy are limited in accuracy by the small 
number of source detections but appear to be broadly consistent with that measured 
for EROs. Nearest-neighbour analyses further support strong clustering on arcmin 
scales, rejecting the null hypothesis that the distribution of the submm sources is ran- 
dom at the 95% confidence level for Sgso > 5 mJy, and at the 99% confidence level for 
Ssso > 7mJy. 

Key words: cosmology: observations - galaxies: evolution - galaxies: formation - 
galaxies: starburst - infrared: galaxies 
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1 INTRODUCTION 

Over the past seven years, a series of complementary deep 
850 /im surveys (eg. Smail et al. 1997, Hughes ct al. 1998, 
Blain et al. 1999, Barger et al. 1998, Barger et al. 1999, 
Bales et al. 2000, Scott et al. 2002, Cowie et al. 2002, Borys 
ct al. 2003, Webb ot al. 2003a) carried out using SCUBA 
(Holland et al. 1999) on the JCMT has resolved the bulk of 
the far-infrared (FIR) extragalactic background into discrete 
sources. These surveys vary in size and depth from ultra- 
deep surveys exploiting gravitational lensing from interven- 
ing clusters to study the very faintest submm sources (Smail 
et al. 1997, Cowie et al. 2002), small and deep blank field sur- 
veys such as the HDF (6 sq. arcmin to a uniform farms — 0.5 
mjy/beam; Hughes et al. 1998, Serjeant et al. 2003), through 
to moderate area and comparatively shallower blank field 
surveys such as the "SCUBA 8-mJy Survey" (a total of 250 
sq. arcmin to a uniform IcTrms — 2.5mJy/beam; Scott et al. 
2002). These surveys have revealed a population of heav- 
ily dust-enshrouded high-redshift {z > 1) morphologically- 
irregular galaxies. Although there is still considerable un- 
certainty in the fraction of SCUBA sources hosting either 
low luminosity or Compton-thick AGNs (Smail et al. 2002, 
Ivison et al. 2002) , deep X-ray observations with the Chan- 
dra and XMM-Newton telescopes have suggested that even 
when an AGN is present, it rarely dominates the far-infrared 
/ submillimctre emission from the galaxy (Prayer ct al. 1998, 
Alexander et al. 2003). Hence, if the thermal dust emission is 
dominated by reprocessed starlight, the inferred star forma- 
tion rates (SFRs) in the very brightest submillimctre galax- 
ies may be as high as several thousand solar masses per 
year. This is sufficient to form the most massive elliptical 
galaxies observed in the present-day Universe on timescales 
of ~ 1 Gyr, implying that the submillimctre sources may 
be the progenitors of todays massive ellipticals, and that 
the star-formation rate density in the early Universe was at 
least a factor of 2 larger than that derived from optical / 
UV observations (Stcidcl ot al. 1999). 

A key test of whether the bright submillimctre sources 
really are the progenitors of present-day massive spheroids 
is to measure their clustering properties. If they are indeed 
proto-massive ellipticals then they should be strongly clus- 
tered. This is an inevitable result of gravitational collapse 
from Gaussian initial density fluctuations; the rare high- 
mass peaks arc strongly biased with respect to the mass 
(eg. Benson et al. 2001). Each of the submillimctre survey 
consortia have performed their own reduction, source extrac- 
tion and simulations on their individual datasets, in order to 
study the nature of the SCUBA population in general. How- 
ever, each survey alone images only ~ 100 — 200 sq. arcmin 
of sky spread over several fields, resulting in discrepancies 
in the cumulative number counts at Ssso > 8 mjy by over a 
factor of 5 due to cosmic variance, and also potentially the 
effects of clustering. Furthermore, no individual survey has 
identified enough sources to make a significant measurement 
of the clustering properties, although tentative evidence was 
obtained for strong clustering on scales of 1-2 arcmin from 
the "SCUBA 8mJy Survey", the largest of these surveys 
undertaken to date (Scott et al. 2002). 

In this paper we combine the data from all of the blank 
field surveys completed up to the completion date of the 
"SCUBA 8 mJy Survey" . In order to do this we downloaded 



the raw data for the "Canada UK Deep Submillimctre Sur- 
vey (CUDSS)", the "Hawaii Flanking Fields Survey" and 
the "Hubble Deep Field (HDF) North Pencil Beam Survey" 
from the Canadian Astronomy Data Centre (CADC) archive 
and re-reduced it in an identical manner to that employed 
in the "SCUBA 8 mJy Survey" . This is decribed in Section 
2. The source extraction algorithm developed by Scott ot al. 
(2002) was used to identify significant submillimctre sources 
in each of the fields, and this method is described briefly in 
Section 3. In Section 4 we present and discuss the results 
of simulations carried out both in conjunction with the real 
data, and through the production of fully-simulated survey 
areas, for all of the survey fields. Comparative reduction 
methods and source lists down to a signal-to-noise ratio of 
3.00 are given in Section 5. These lists include some previ- 
ously unidentified significant sources but also CEist doubt on 
the reality of other published detections. Section 6 combines 
these new catalogues into a master source list, from which 
the most accurate cumulative number counts to date in the 
flux density range 2—12.5 mJy have been calculated. A num- 
ber of models are tested against these data points. In Section 
7, the new master catalogue has been used in an attempt to 
measure the clustering properties of the bright (> 5mJy) 
submillimctre sources, through 2-point angular correlation 
fuuctions and nearest-neighbour analyses. Section 8 outlines 
the resulting motivation and aims of the current blank field 
survey work being undertaken by the JCMT -|- SCUBA, and 
the conclusions from this work are presented in Section 9. 



2 DATA REDUCTION 

The raw data for the "Canada UK Deep Submillimctre Sur- 
vey (CUDSS)", the "Hawaii Flanking Fields Survey", and 
the "Hubble Deep Field Survey" (pencil beam only) were 
downloaded from the Canadian Astronomy Data Centre 
(CADC) archive. The data were then fully reduced using 
the same Interactive Data Language (IDL)-based reduction 
routines employed in the "SCUBA 8 mJy Survey" . This pro- 
cedure is fully described in Scott et al. (2002), but is sum- 
marised briefly below. 

Firstly, in order to take nodding into consideration, 
the off-position was subtracted from the on-position in the 
raw beam-switched data. The relative sensitivities of the 
bolometers, with respect to a reference bolometer (the cen- 
tral bolometers H7 and C14, for the 850 /im and 450 fim 
arrays respectively) were accounted for by multiplying by 
the standard flatfield values. The atmospheric opacity was 
measured wherever possible using a 6th order polynomial fit 
to the 225 GHz measurements from the Caltech Submillimc- 
tre Observatory (CSO), followed by a linear interpolation to 
850 and 450 /um using the relations given in Archibald et 
al. (2002a). On the occasions when the CSO opacity meter 
was out of service, or less than 7 reliable CSO observations 
were available within ±1 hour of each observation, a linear 
interpolation between successive skydip values was used to 
determine rssoMm and r450Mm. 

The next stage in the data reduction process was to 
remove spikes in the data resulting from cosmic-rays and bad 
bolometers, and to remove any residual sky emission. In this 
IDL reduction, deglitching and residual sky subtraction were 
undertaken by an iterative process, each iteration making 
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a temporal noise estimate and deglitching, followed by a 
spatial sky subtraction. There were no bright sources in any 
of submm survey fields that would have been significantly 
detected in any single jigglemap, let alone in sub-dividing 
the data-stream into shorter timescale chunks and so the 
procedure was as follows: 

• For each bolometer, noise estimates were made by fitting 
a Gaussian to the data-stream in chunks of 128 readout 
groups. 

• These time dependent noise estimates were then used to 
remove any spikes by performing a Scr clip on the data. 

• Using the fits to all of the bolometers in the array, a modal 
residual sky level was determined for each of these 128 read- 
out groups, and subtracted from these data. 

With each consecutive iteration, the deglitching process 
makes a harder cut. Noisy bolometers were assigned a low 
inverse variance weight in this way. In the case of the cali- 
bration data, however, the presence of a bright source would 
likely lead to over-enthusiastic clipping of the data, and so in 
this case a timeline without object signal was constructed. 
This was created by calculating the mean of the timestream 
data points recorded immediately before and after the read- 
out being considered, and subtracting this from the readout 
value. 

Flux conversion factors (FCFs) were deterniiucd by di- 
viding the flux density within the JCMT beam of a known 
calibration source by the measured peak voltage. Each of the 
individual jiggle-maps comprising the submm survey data 
were calibrated prior to producing the final coadded images 
using the gain value from whichever calibration source was 
taken closest in time. 

The final images were produced using an optimal noise- 
weighted drizzling algorithm (Fruchtor & Hook, 2002) with a 
pixel size of 1 sq. arcsoc. This is the same method as that em- 
ployed in the 'SCUBA 8-mJy Survey' (Scott ot al. 2002) and 
'Hubble Deep Field North' (Serjeant et al. 2003) data reduc- 
tions. Both output signal and noise maps were created, the 
signal in any one sq. arcscc pixel given by the noise- weighted 
average of the bolometer readouts at that position, and the 
corresponding noise value given by a noise-weighted average 
of the Gaussian fits to the readout histograms. Unlike a stan- 
dard shift-and-add technique which takes the flux density in 
each detector pixel and places it into the final map over an 
area equivalent to one detector pixel projected on the sky, 
drizzling takes the flux density and places it into a smaller 
area in the final map. Although this significantly reduces 
the signal-to-noise ratio in each pixel, this approach helps 
preserve information on small angular scales, provided that 
there are enough observations to fill in the resulting gaps. 
The area in the coadded map receiving the flux from one 
detector pixel is termed the 'footprint'. This method is an 
extreme example of drizzling in which the 'footprint' is se- 
lected to be as small as is practicable given the pointing 
errors invloved (termed the 'zero-footprint'), selected here 
to be one sq. arcsec. Unlike in the standard SURF reduc- 
tion, there is no intrinsic smoothing or interpolation between 
neighbouring pixels in this rebin procedure. Although there 
is some degree of correlation between pixels in the output 
zero- footprint signal maps in terms of the beam pattern, the 
corresponding pixel noise values represent individual mea- 
surements of the temporally varying sky noise averaged over 
the dataset integration time, at a specific point on the sky. 



and are hence statistically independent from their neigh- 
bours. In essence this method produces a very ovcrsamplod 
image with statistically independent pixels. A final 4CT-clip 
on the signal-to-noise was carried out to remove any remain- 
ing 'hot pixels'. A noise- weighted convolution with a beam- 
sized Gaussian point spread function (PSF) was used to pro- 
duce realistic smoothed maps of the survey areas whilst ac- 
counting for variable signal-to-noise between individual pix- 
els. 



3 SOURCE EXTRACTION 

The chopping-nodding mechanism of the telescope provides 
a valuable method of discriminating between real detections 
and spurious noise spikes in the data. With the exception of 
the Hubble Deep Field (HDF), each of the surveys used a 
single chop throw fixed in right ascension (RA), thus creat- 
ing negative sidelobes, half the depth of the peak flux den- 
sity, on either side of a real source. In the case of the HDF, 
this strategy was modified to use two chop throws fixed in 
RA, each chop throw used for approximately half of the to- 
tal integration time. This side-lobe signal can be recovered 
to boost the overall signal-to-noise ratio of a detection. 

For well-separated sources, convolving the images with 
the beam is formally the best method of source extraction 
(Eales et al. 1999, 2000, Serjeant et al. 2003). However, fol- 
lowing a careful examination of the reduced "8 mjy Survey" 
data (Scott et al. 2002), it became clear that some of the 
potential sources were partially confused. This was partic- 
ularly prominent in the map of ELAIS N2 where the nega- 
tive sidelobes of individual sources have overlapped and are 
therefore somewhat deepened relative to both source peaks. 
Consequently, in order to decouple any confused sidelobes a 
source-extraction algorithm was devised based on a simul- 
taneous maximum-likelihood fit to the flux densities of all 
potentially significant peaks in the maps. This is made feasi- 
ble by the independent data-points and errors yielded by the 
zero-footprint IDL-reduced maps. These peaks were identi- 
fied as any positive peak in the noise- weighted Gaussian con- 
volved signal maps. Using a peak-normalised beam-map as 
a source template (generated by binning together all of the 
observations of Uranus or CRL618 taken with the relevant 
chop throw), a basic model was constructed by centring a 
beam-map at the positions of every peak in the maps. The 
normalisation coefficients of each of the positioned beam- 
maps were then calculated simultaneously such that the fi- 
nal multi-source model provided the best description of the 
submm sky, as judged by a minimum fit- 

The fitting process is as follows. Suppose one considers 
a normalised beam-map B{x,y) as a source template and 
that at position (i, j) in the unconvolved image the signal 
is S(i,j) and the noise is N{i, 'j). If n peaks above a spec- 
ified flux threshold are located in the Gaussian-convolved 
image, one may construct a model to the unconvolved zero- 
footprint image such that beam-maps centred on each peak 
position are simultaneously scaled to give an overall best fit 
to the entire image. Using a minimum At as the maximum 
likelihood estimator then 

2 (x^S{i,j)-J2l^i"'kBk{x-i,y-j)\ 

^=2- mJr 
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Survey 


Total area 


Uniform noise 


1(7 rms noise 


Field 


/sq. arcmins 


area /sq. arcmins 


level mjy/beam 


ELAIb N2 


137 


113 


2.2 ± 0.7 


Lockman Hole wide area 


141* 


88 


2.7 ± 0.7 


Lockman Hole deep strip 




21 


1.8 ± 0.2 


CUDSS 03h wide area 


69* 


55 


1.8 ± 0.5 


CUDSS 03h deep area 




8 


1.1 ± 0.2 


CUDSS lOh 


10 


8 


1.3 ± 0.2 


CUDSS 14h 


61 


57 


1.5 it 0.3 


CUDbb 22n 


7 


5 


1.5 ± 0.3 


Hubble Deep Field 


10 


6 


0.6 ± 0.1 


SSA13 wide area 


72* 


45 


2.5 ±0.6 


SSA13 deep area 




8 


0.7 ±0.1 


SSA17 


24 


21 


1.6 ±0.5 


SSA22 


26 


21 


0.9 ±0.2 


Lockman Hole deep area 


11 


8 


0.8 ±0.1 



Table 1. Survey field areas and Icr rms noise levels in the regions of uniform noise, as given by the mean and standard deviation 
measured directly from the 14.5" FWHM Gaussian convolved noise maps. Total areas marked with * refer to fields composed of a small 
deep region within a wider shallower survey area and correspond to the full area of that entire field (ie. both shallow and deep). 



where Ofc is the best fit flux to the kth peak. 
Defining amk, an n x n matrix, as 



amk = 



Bk{x-i,y- j)Bm{x -i,y- j) 



and Prn, a vector of length n as 
o _sr S{i,j)Bm{x - i,y-j) 

then the best fit values of aj, are given by 

n 



(2) 



(3) 



(4) 



and the significances a{ak) of the peak detections are given 
by 



o'(afe) 



ak 



(5) 



Furthermore, this method can bo modified to deal with 
surveys which have used more than one chop throw or po- 
sition angle. The peaks are found in the same way as be- 
fore, by regridding all of the individual observations together 
(regardless of the chop throw or position angle used) and 
carrying out a noise-weighted smoothing with a beam-sized 
Gaussian. When conducting the fit, however, each partic- 
ular combination of chop throw and position angle is binned 
separately. If r different chop configurations have been used 
then the expression for becomes 




Sp{i,j) - J2k=i <^kBp,k{xp - i,yp - j) 
Npiijr 



(6) 



and 



amk=^iY^ 



=1 \ i,j 



Bp,k{xp - i,yp - j)Bm,k{xp - i,yp - j) 
Npiijy 



(7) 




Spii, j)Bp^rn(xp — i,yp — j) 

Npiijr 



(8) 



where the expressions for the best fit values of at and sig- 
nificances a(ak) are the same as before. 

For a fuU mathematical description of this source ex- 
traction algorithm, the reader is referred to Scott et al. 
(2002) and Mortier et al. (2005). 



and 



4 SIMULATIONS 

In order to assess the effects of confusion and noise on the 
reliability of the source-extraction algorithm, Monte Carlo 
simulations were carried out on all of the survey fields. The 
individual fields vary widely in size and depth from small, 
deep surveys covering a few sq. arcmin of sky down to the 
confusion limit (eg. the Hubble deep field), to wider, shal- 
lower surveys aimed at studying the most luminous sub- 
millimetre sources on scales of ~ 100 sq. arcmin (eg. the 
wide area Lockman Hole field from the "SCUBA 8-mJy Sur- 
vey"). The typical noise levels and areas of each of the fields 
are given in Table 1. The dependences of positional error, 
completeness and error in reclaimed flux density, on input 
source fiux density and noise in the maps, were determined 
by planting individual sources of known flux density into the 
real SCUBA maps. This has the advantage of testing the 
source-reclamation process against the real noise and con- 
fusion properties of the images, accounting for any cluster- 
ing in the faint background source population, for example. 
However, these simulations do not allow assessment of the 
level to which false or confused sources can contaminate an 
extracted source list. Therefore a number of fully-simulated 
images of the survey areas have also been created by assum- 
ing a reasonable 850 /itm source-counts model, derived from 
a best-fit power-law to the source counts given later in this 
paper. The results of analyzing these sets of simulations are 
discussed in the following two subsections. 
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Survey Field 


Noise Region 




a 


b 




ELAIb N2 


uniform 





17043 


2.3338 


0.37425 


Lockman Hole wide area 


uniform 





14800 


3.2690 


0.28300 


Lockman Hole deep strip 


uniform 





20896 


2.1064 


0.22527 


CUDS& 031i wide area 


uniform 





20081 


2.0013 


0.21037 


CUDSS 03h deep area 


uniform 





28937 


1.6500 


0.39179 


CUDSS lOh 


uniform 





22587 


1.3248 


0.42789 


CUDSS 14h 


uniform 





21016 


1.5467 


0.31834 


CUDSS 22h 


uniform 





26450 


1.9359 


0.12883 


Hubble Deep Field 


uniform 





39711 


0.6258 


0.15015 


SSA13 wide area 


uniform 





14627 


2.6614 


0.19068 


SSA13 deep area 


uniform 





22037 


0.8738 


0.13963 


SSA17 


uniform 





24248 


1.9077 


0.19140 




uniiorm 





21362 


U.4440 


O.o (69 / 


Lockman Hole deep area 


uniform 





24862 


0.8042 


0.12062 


ELAIS N2 

Lockman Hole wide area 


non-uni 
non-uni 


0.00588 
0.01388 


0.0030 
4.6093 


0.04114 
0.05206 



Table 2. Best fit values determined for a and b in equation 9, describing the percentage differential completeness against input source 
flux density for each of the survey fields and noise regions. 
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Rms noise mjy/beom 

Figure 1. Best-fit values for the parameter "a" as given in Table 
2, plotted against the Icr rms noise levels as determined from 
the uniform regions of the 14.5" Gaussian convolved noise images 
(values given in Table 1). The horizontal error bars show the 
standard deviation of the noise about the mean level. 
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12 3 

Rms noise mJy/beom 

Figure 2. Best-fit values for the parameter "b" as given in Table 
2, plotted against the Icr rms noise levels as determined from 
the uniform regions of the 14.5" Gaussian convolved noise images 
(values given in Table 1). The horizontal error bars show the 
standard deviation of the noise about the mean level. 



4.1 Simulations building on the real survey data 

A normalised beam-map, with the same chop throw and po- 
sition angle as that used in the real observations, was used 
as a source template. At flux density intervals of 0.5 mjy, 

spanning the entire range of flux densities for which real 
sources were recovered, fake sources were added into the 



unconvolved zero-footprint signal maps. This was done one 
fake source at a time, so as not to enhance significantly any 
existing real confusion noise within the image. The source- 
extraction algorithm was then re-run. This exercise was re- 
peated for 100 different randomly-selected positions on each 
image, at each flux density level, so that source reclamar 
tion could be monitored as a function of input flux density 
and positiou/noisc-lcvcl within the maps. The source recla- 
mation was deemed to have been successful if the source- 
extraction algorithm returned the fake source with signal- 
to-noise > 3.50 (a level selected as a compromise between 
recovering a reasonable number of sources and contamina- 
tion with spurious / confused sources - sec Subsection 4.2) 
within less than half a beam-width of the input position, 
but excluding from the analysis any fake sources which had 
fallen upon a position within half a beam- width of a brighter 
> S.OOcr peak already detected in the map. This is because 
the flux densities of the recovered sources within the real 
data span a broad range (~ 2 — 12 mJy), and it is not pos- 
sible to resolve two separate sources placed closer together 
than this - they would appear as one peak in the Gaussian- 
smoothed image. It is not realistic to consider, for example, 
a fake 2mJy source to have been successfully recovered if 
it lies almost on top of an 8mJy source already detected 
significantly in the map. Under this situation it is really the 
8mJy source already present in the image which is being 
recovered. Reversing the situation, however, the successful 
reclamation of a fake 8mJy source planted into the map in 
the near vicinity of an already significantly detected 2 mJy 
source (a possible scenario in the very deep images such as 
that of the HDF) would be included in the analysis because 
the fake source is making the dominant contribution to the 
combined fiux density. 

It is possible that future interferometers may resolve 
some of our point sources into multiple components; this 
is an inevitable caveat to any source count analysis. In the 
meantime however, the JCMT resolution provides an effec- 
tive working definition of point source for the current work. 

The Lockman Hole East field from the "SCUBA 8 mJy 
Survey" , the 03h field from the "Canada UK Deep Submil- 
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limetre Survey (CUDSS)", and the SSA13 field from the 
"Hawaii Submillimctrc Survey" , contain sections of map 
which are markedly deeper than the rest of the data. In the 
case of the Lockman Hole this was due to an early change in 
survey mapping strategy, and in the 03h and SSA13 fields 
this resulted from a deep pencil beam survey being incor- 
porated into the wider-area images. In each of these cases, 
separate sets of simulations were run on the deep and shal- 
lower sections of the fields. 

Additionally, regions of uniform and non-uniform noise 
were defined for each field (again treating the deep parts of 
the Lockman Hole, 03h and SSA13 fields as separate fields 
from the wider-area shallower part), using the "GAIA" tool 
to manually cut out a template of the uniform noise area 
using the Gaussian-smoothed noise maps. The deep pen- 
cil beam surveys, such as the HDF and CUDSS 22h field, 
are comprised of a stack of jiggle-map observations centred 
on one or two positions only, and so the non-uniform edge 
regions in these images are largely the result of undersam- 
pling from the bolometers on the outer ring of the array. The 
wider-area images, however, were built up from a series of 
jiggle-pointings, offset from each other by some fraction of 
an array width. Hence, the pointings forming the outer-most 
regions of the survey field lack the next consecutive set of 
integrations from what would have been the neighbouring 
pointing, resulting in a border of shallower (and hence nois- 
ier) observations. In Section 5, which discusses the various 
survey fields in detail, any sources recovered in these poorer 
noise regions have been marked with the term "edge" in the 
source list tables. 

4-1.1 Completeness 

The differential completeness is given by the percentage of 
sources recovered with signal-to-noise ratio > 3.50 at each 
input flux density level, and was found to be described well 
by the functional form 

differential completeness = 100(1 - e"''^''"''') (9) 

where x is the input flux density, and the values of a and 
b were determined by a minimised fit to the simulation 

results for each 850 nm survey field. The values of a and b 
determined from these fits arc given in Table 2. The plots 
for the percentage of sources recovered against input flux 
density for the individual survey fields may be found in Ap- 
pendix Al. 

The primary goal in allowing a fit of this nature was to 
obtain a best-fit description of the overall shape of the curve, 
rather than a detailed analysis of possible combinations of 
free-parameters 'a' and 'b' in space. However, even sim- 
ple plots of the best fit values of a and b against the la 
rms noise levels as measured from the beam-sized Gaussian 
convolved noise images (Figs. 1 and 2 respectively), show 
clear noise-dependent trends. The horizontal error bars re- 
flect the standard deviation of the noise values about the 
mean, in the uniform regions of the map. The best-fit val- 
ues of parameter 'a' show a general decrease with increasing 
rms noise levels, albeit with a fairly broad dispersion, par- 
ticularly between the deep pencil beam surveys such as the 
Hubble deep field and the SSA13 deep area field. This is 
likely a combination of being at the confusion limit (gener- 
ally high source density) and the variation in the number 



density of sources between these small area fields (cosmic 
variance and perhaps clustering effects also). The param- 
eter 'b' defines a lower flux density cut-off below which no 
sources are successfully recovered, and shows a much tighter 
correlation, increasing roughly linearly with the rms noise as 
b ~ 1.25 X noise. 

Unfortunately, the scatter in parameter 'a' with rms 
noise is too large to allow for a general differential complete- 
ness formula applicable to any survey field to be developed, 
based on these data. 

Figures Al, A2 and A3 show the completeness analysis 
for the uniform noise regions of the "SCUBA 8 m,Jy Survey" 
fields (ELAIS N2, Lockman Hole wide area and Lockman 
Hole deep strip). The error bars are given by the Poisson 
error on the rmmber of sources planted into the field in each 
noise region, and at each flux density. These fields have the 
largest shallow border regions of all the survey fields dis- 
cussed here, due to the survey strategy adopted to even out 
the noise (see Section 5.1). The corresponding completeness 
plots for the non-uniform regions of ELAIS N2 and the Lock- 
man Hole wide area are shown in Figs. A4 and A5 respec- 
tively. It is immediately obvious in comparing plots of uni- 
form and non-uniform noise that source recovery in the non- 
uniform edge regions is very much worse than in the fully- 
observed central areas, reaching at best ~ 10% at 15mJy 
as opposed to the ~ 90% in the uniform noise regions. The 
simulations carried out on the remaining smaller fields did 
not yield sufficiently good statistics in the non-uniform noise 
regions to allow any meaningful fit to be made, hence only 
plots for the uniform noise regions of the remaining fields 
have been presented (Figs A6 to A16). 

4- 1.2 Output versus Input Flux Density 

Using these simulations, it is also possible to determine the 

dependence of the mean output-to-input flux density ratio 
as a function of the input flux density, for those sources 
identified with signal-to-noise ratio > 3.50. This relation was 
found to be well described by the expression 

output flux density ^ ^^_dx _^ ^ 
input flux density 

where x is the input flux density, and the values of C, d and 
f were determined by a minimised At to the simulation 
results for each 850 /xm survey field and are given in Table 3. 
The plots of the ratio of output-to- input flux density against 
input flux density for the individual survey fields may be 
found in Appendix A2. 

The plots of mean output/input flux density ratio 
against input flux density are shown in Figs. A17 to A32. 
The error bars are the standard error on the mean. One of 
the first things to notice about the subsequent ratio plots, is 
that the effect of noise and confusion is to produce system- 
atic 'flux-boosting', the mean retrieved flux density always 
being greater than the input value. This effect is known as 
Eddington bias (Eddingtou 1913) and is apparent in any 
flux limited survey where a specific signal-to-noise thresh- 
old is employed. The presence of noise and confusion from 
the faint background population will vary the flux densi- 
ties with which a source of specified input flux density is 
retrieved. If, for example, one considers a very simple case 
of pure Gaussian noise on a fake source, the measured flux 
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Survey Field 


Noise Region 


C 


d 




f 




ELAIS N2 


uniform 


6.350 


0.4813 


1 


.0150 


1.11683 


Lockinaii Hole wide area 


uniform 


5.555 


0.3591 


1 


.0246 


1.08116 


Lockman Hole deep strip 


uniform 


6.543 


0.5313 


1 


.0154 


1.34420 


03h wide area 


uniform 


6.587 


0.5667 


1 


.0603 


1.05373 


03h deep area 


uniform 


13.575 


1.1842 


1 


.0537 


1.67181 


lOh 


uniform 


8.786 


0.9443 


1, 


.0439 


0.76535 


14h 


uniform 


8.086 


0.7701 


1, 


.0325 


1.71510 


22h 


uniform 


9.815 


0.8236 


1 


.1301 


1.51353 


Hubble Deep Field 


uniform 


12.106 


1.6358 


1 


.0513 


1.80567 


SSA13 wide area 


uniform 


4.606 


0.4312 


1, 


.0430 


1.52413 


SSA13 deep area 


uniform 


2.726 


0.8049 


1, 


.0516 


1.22079 


SSA17 


uniform 


6.769 


0.8624 


1, 


.0965 


1.60692 


SoAzz 


uniform 


8.262 


1.2241 


1, 


.0168 


1.23946 


Lockman Hole deep area 


uniform 


9.712 


1.5104 


1, 


.0422 


0.68671 


ELAIS N2 

Lockman Hole wide area 


non-uni 
non-uni 


3.686 
112.821 


0.2283 
0.6582 


1.1805 
1.2593 


1.84850 
3.39015 



Table 3. Best fit values determined for C, d and f in equation 10, describing the output to input flux density ratio against input flux 
density for each of the survey fields and noise regions. 



densities would be expected to have a symmetric distribu- 
tion about the actual source flux density, the exact charac- 
teristics of the distribution dependent on the level of noise 
applied. However, if a fixed signal-to-noise ratio is applied to 
the source extraction procedure, one will preferentially se- 
lect those sources which have been retrieved with a brighter 
flux density, as some of the fainter measured values will fail 
to make the signal-to-noiso cutoff'. Consequently, the mean 
retrieved flux density will always be larger than the input 
flux density. Applying the same noise characteristics to in- 
put sources of increasing brightness, the mean boosting ra^ 
tio is reduced, because only the larger negative fluctuations 
on the tail of the Gaussian noise distribution will allow the 
brighter sources to fall below the signal-to-noise threshold. 
For very bright sources the output-to- input flux density ratio 
approaches unity. Non-Gaussian noise and confusion will of 
course affect the distribution of the retrieved flux densities - 
in particular confusion of faint background sources may lead 
to a more asymmetric distribution, especially if the SCUBA 
population is found to strongly cluster. Simulations such as 
these, however, allow for an empirical numerical description 
on a field by field basis. 

Trends in the properties of parameters 'C, 'd' and 'f 
with rms noise are shown in Figs. 3, 4 and 5 respectively. 
Both parameters 'C and 'd' decrease with increasing rms 
noise. The decline is steep at low rms noise levels, but be- 
comes more shallow above farms — 1-5 — 2mJy. Parameter 
'd' shows a fairly tight correlation, however 'C shows too 
great a level of scatter to allow a general formula, based 
solely on rms noise, to be developed for the output-to-input 
flux density ratio. The parameter 'f represents the ratio of 
output-to-input flux density for very bright sources with a 
constant value of ~ 1 expected for all fields, regardless of 
noise level (the median value is in fact 1.04). 

It can also be readily seen from comparing Figs. A17 
and A18, with A20 and A21, that the level of flux-boosting 
is much greater in the non-uniform noise regions and with a 
much larger degree of scatter in the data points. For exam- 
ple, a source input to the ELAIS N2 or Lockman Hole fields 
("from the SCUBA 8-mJy Survey") would appear boosted 



on average by a factor of 1.2 — 1.3 if extracted from the uni- 
form noise regions. In the non-uniform regions though, the 

mean level of boosting is by a factor 2. Duo to the combi- 
nation of a poor level of retrieval and large flux boosting 
factors, any sources recovered in the non-uniform noise re- 
gions (marked as "edge" in Section 5) have been excluded 
from statistical analyses such as source counts, and cluster- 
ing measures etc. 

4-1.3 Positional Uncertainty 

The mean positional uncertainty in retrieving the fake 
sources was found to be well approximated by a linear de- 
pendence on the input flux density such that 

positional error = — gx + h (11) 

where x is the input flux density, the values of g and h for 
each 850 /^m survey field were determined by a minimised 
fit to the simulation results (given in Table 4), and the 
positional uncertainty is given in arcseconds. The plots of 
the mean positional uncertainty against input flux density 
for each of the individual fields are given in Appendix A3. 

Figs. 6 and 7 show the dependence of parameters 'g' 
and 'h' on rms noise. One might expect a general formula 
for positional error to depend on the ratio of input flux 
density to rms noise such that ggoneral oc Icr^s in this 
straight line desription. The data points are consistent with 
ggenerai X Ifrms ~ 0.35 arcscc, but the large scatter means 
that a simple straight line with negative gradient provides a 
similarly good description. The values of parameter 'h' have 
a median of ~ 4.65 between all the survey fields, and there 
is no obvious trend with rms noise. 

Figures A33 through to A48 show the mean positional 
error of the retrieved fake sources against input flux density. 
The error bars are the standard error on the mean. Again, 
one can see that the data from the "SCUBA 8 mjy Survey" 
fields (Figs. A33 to A37) show a greater scatter in the non- 
uniform noise regions than the central uniform noise parts, 
with a generally greater positional uncertainty in these bor- 
der areas 4 — 5 arcsec, as compared with 2 — 4 arcsec in the 
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Survey Field 


Noise Region 


g 


h 








arcsec/mjy 


arcsec 




ELAIb N2 


uniform 


0.17364 


4.6693 


1.45445 


Lockman Hole wide area 


uniform 


0.12235 


4.5791 


0.90109 


Lockman Hole deep strip 


uniform 


0.16592 


4.2357 


1.65490 


03h wide area 


uniform 


0.13772 


4.2672 


1.05487 


03h deep area 


uniform 


0.46026 


5.2313 


2.22813 


lOh 


uniform 


0.36036 


4.8814 


2.04834 


14h 


uniform 


0.20604 


4.2997 


0.52830 


22h 


uniform 


0.33656 


5.6634 


0.76516 


Hubble Deep Field 


uniform 


0.42609 


5.0274 


1.11186 


SSA13 wide area 


uniform 


0.11505 


4.2934 


0.76635 


SSA13 deep area 


uniform 


0.43070 


5.1689 


0.36683 


SSA17 


uniform 


0.37809 


5.1600 


0.65348 


SSA22 


uniform 


0.29406 


4.3191 


1.79829 


Lockman Hole deep area 


uniform 


0.33253 


4.2634 


0.54652 


ELAIS N2 


non-uni 


0.51966 


9.4885 


4.09028 


Lockman Hole wide area 


non-uni 


0.00000 


4.1870 


1.88165 



Table 4. Best fit values determined for g and h in equation 11, describing the mean positional error against input flux density for each 
of the survey fields and noise regions. 



uniform noise regions). The positional accuracy improves 

with higher flux density sources (and hence better signal-to- 
noise). Such estimates of positional error do not include any 
pointing errors arising whilst the data are being taken at the 
telescope. The pointing of the JCMT is known to be very 
accurate for such a large dish; typical pointing errors are less 
than 3 aresec, much less than the 850 /im beam size of 14.5 
arcsec FWHM. Several pointing problems have been discov- 
ered during the period in which these deep submillimetre 
surveys were undertaken. However, the fact that each map 
is built up from many shorter integration datasets limits 
the impact which these pointing inaccuracies can have on 
the final image. Overall, the uniform noise regions of all the 
survey areas suggest likely positional errors of 2 — 4 arcsec 
in source retrieval, arising from the effects of noise and con- 
fusion. 

Note that the positional uncertainty is known to be well- 
described by a Gaussian distribution in this data reduction 
methodology (e.g. Ivison et al. 2002, Serjeant et al. 2003, 
Coppin et al. 2006). 

4.2 Completely simulated maps 

In order to obtain constraints on the fraction of recovered 
850 ^m sources arising from confusion and noise, 100 simu- 
lated images of each of the survey fields were generated. In 
addition, these simulations were also used to estimate the 
integrated completeness and count correction factors. The 
assumed source counts were taken from the best fit of a 
simple power-law model (in the format first employed by 
Barger et al., 1999) to the differential counts given later in 
this paper, which were corrected for completeness and flux- 
boosting at the > 3.50(T level using the simulation results of 
Section 4.1. Specifically, the differential counts are given by: 

dN(S) _ No 
dS (a-|-S«) ' 

where No = 2.67 x 10*, a = 0.49 and a = 3.14, which 
predicts a total 850 /xm background of 3.8 x 10'* mjy deg~^, 



consistent with the value of 4.4 x 10*mJydeg~^ measured 

by Fixsen et al. (1998). A realistic model of the background 
counts was produced , by randomising tiie number of sources 
placed into each simulated field at 0.1 mjy intervals, from 
0.1 — 14.0 mJy, according to a normal distribution about the 
number expected. Each source was then allocated a random 
position and the whole image was convolved with the beam. 
The simulated field is initially created to be larger than the 
actual field, allowing for tiie negative sidelobes of sources 
centred off-field in the final image to appear in it. The clus- 
tering properties of the SCUBA population are, at present, 
not well characterised. Results presented in Section 7 suggest 
that at least the very brightest SCUBA sources (> 5mJy) 
are strongly clustered on arcminute scales, consistent with 
the idea that these objects are progenitors of present-day 
massive ellipticals, but there is insufficient blank field sur- 
vey data available to allow a realistic clustering component 
to be added into the selection of positions within the simu- 
lations. This means that these simulations which make the 
assumption of a random distribution (ie. no clustering) can 
only be used as a first approximation in determining the 
level of spurious / confused source contamination. 

Noise overlays were constructed by subtracting the full 
minimised fit model (i.e. the model representation of the 
full sky region comprised of the best-fit beam profiles to all 
of the peaks in the image) from the zero-footprint signal 
maps of the actual survey data. 

The subsequent simulation analyses were conducted at 
various signal-to-noise levels down to very low significance 
levels (> 1.50ct), even though the source catalogues pre- 
sented in Section 5 only reach > S.OOcr. This was to allow an 
assessment of whether sources recovered at > Scr by other 
submillimetre groups using different reduction and extrac- 
tion methods, but which were recovered at lower signal-to- 
noise in the analysis presented in Section 5, were likely to 
be real. 

At first glance one might think that this process would 

remove a significant amount of real Gaussian noise as well as 
faint sources from the signal image. However, one is not sim- 
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Rms noise mjy/beam 

Figure 3. Best-fit values for the parameter "C" as given in Table 
3, plotted against the Icr rms noise levels as determined from the 
uniform regions of the 14.5" Gaussian convolved noise images 
(values given in Table 1). The horizontal error bars show the 
standard deviation of the noise about the mean level. 



Figure 6. Best-fit values for the parameter "g" as given in Table 
4, plotted against the la rms noise levels as determined from 
the uniform regions of the 14.5" Gaussian convolved noise images 
(values given in Table 1). The horizontal error bars show the 
standard deviation of the noise about the mean level. 



Figure 4. Best-fit values for the parameter "d" as given in Table 
3, plotted against the la rms noise levels as determined from 
the uniform regions of the 14.5" Gaussian convolved noise images 
(values given in Table 1). The horizontal error bars show the 
standard deviation of the noise about the mean level. 



Figure 7. Best-fit values for the parameter "h" as given in Table 
4, plotted against the la rms noise levels as determined from 
the uniform regions of the 14.5" Gaussian convolved noise images 
(values given in Table 1). The horizontal error bars show the 
standard deviation of the noise about the mean level. 



Figure 5. Best-fit values for the parameter "f" as given in Table 
3, plotted against the la rms noise levels as determined from 
the uniform regions of the 14.5" Gaussian convolved noise images 
(values given in Table 1). The horizontal error bars show the 
standard deviation of the noise about the mean level. 



ply fitting a Gaussian to positive flux density peaks (wliicli 
would indeed remove a significant amount of the real noise) , 
but instead fitting tiie full beam profile which has the two 
negative sidelobes, each half the depth of the peak, to the 
signal image. Consequently, if a peak arose due to Gaus- 
sian noise rather than a real source, one would not find the 
accompanying sidelobes at the relevant position and of the 
right depth and so the best-fit normalisation of the beam 
profile at this position would be very close to zero, thus re- 
sulting in very little of the real noise being removed. 

The source-extraction algorithm was then re-run on 
these residual images to determine the number of sources 



which could be recovered from the noise overlays alone, at 
signal-to-noise thresholds of 1.50—4.00(7, spaced regularly at 
0.5(7 levels. Gaussian statistics predict that given the num- 
ber of beams in the 464 sq. arcmin of uniform noise, there 
are likely to be ~ 14 

noise peaks recovered at > S.OOcr and ~ 0.32 noise peaks 
recovered at > 4.00cr (in fact slightly less than this, as this 
calculation does not account for the recovery of the nega- 
tive sidelobes). The actual rmmbcrs recovered are 24 and 
1 at > 3.00a and > 4.00(7 respectively, comparing reason- 
ably well with the Gaussian estimates. Additional sources of 
high-significance (> 3.00(7) peaks in the residuals might be: 

1) A non-Gaussian component in the noise eg. microphonics. 

2) A poor fit of the model to the data in a small sub-region 
of the full dataset. 

3) Incomplete source removal, for example a faint source 
confused with a bright source such that only the brighter 
of the two sources could be identified by the presence of a 
peak. 

Future improvements to the source extraction algorithm 
will address points (2) and (3). Currently, however, a poor 
model-fit to sub-sections of the original map or the presence 
of any remaining real sources in the residual images will lead 
to an over-estimate of the level of spurious/confused source 
contamination, and so the results presented in subsequent 
tables may be considered an upper limit. 
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The final signal images were constructed by adding 
the unsmoothod noiso-ovcrlay to the simulated background 
counts, and trimming to the correct size and shape. The orig- 
inal zero-footprint noise maps were used as noise maps for 
the simulated images, and any "hot" pixels identified with 
signal-to-noise ratio > 4.00 were re-assigned large noise lev- 
els, as was done with the real data. The source extraction 
algorithm w£is then applied to each simulated image in an 
identical way to the actual survey maps. 

Simulated images of the 03h (CUDSS), SSA13 (Hawaii 
Survey) and Lockman Hole wide area (8mJy Survey) fields 
were created with the small deep regions combined into the 
wider area surveys, however, in the same way as the adding 
of one source into the real data and attempting to retrieve 
it (Section 4.1) the results for the deep and shallower areas 
were treated separately. The regions of uniform and non- 
uniform noise were also treated individually, as before. 

These simulations difiier slightly to those presented in 
Scott et al. (2002), and this is reflected in the results pre- 
sented for the "SCUBA 8mJy Survey" fields (ELAIS N2 
and the Lockman Hole East wide area field) . Combining the 
data from the various surveys improves the constraints on 
the 850 /xm source counts (as discussed in Section 6), and a 
steeper source counts model, fit to the combined counts, has 
been used to create these mock images. This leads to higher 
densities of fainter sources, and hence incresises the fraction 
of significant detections arising from confusion. The second 
difference is that a lower flux density cut-off of 14 m,]y has 
been applied, corresponding to the retrieved flux density of 
the brightest source detected in the uniform regions of any 
of the survey fields, as opposed to the 20mJy cut-off em- 
ployed in the earlier simulations. The presence of sources 
with artificially high input flux densities increases the frac- 
tion of objects retrieved with high signal-to-noise ratios, 
which in turn overestimates the integral completeness at 
a given signal-to-noise threshold, particularly in the non- 
uniform areas where the noise levels are higher. The third 
difference is the noise overlay added on to the background 
sources. In previous simulations this was created by rebin- 
ning the individual datasets with randomised bolometer as- 
trometry so as to smear out any sources present. This ap- 
proach was found to have problems in regions with several 
significant bright sources, which would become smeared to- 
gether on scrambling, creating a patch of excessive noise. For 
this reason, and in order to preserve the noise properties of 
the real data as far as possible, the residual signal maps 
were adopted as the overlaid noise. These residual maps are 
the difference between the pixel values of the actual un- 
convolved signal maps, and the best-fit model of the full sky 
region as constructed firom a series of idealised beam-profiles 
centred on every peak in the convolved signal image. Hence, 
the residual maps represent the excess noise levels superim- 
posed on top of the real data. The final difference is in the 
fiux densities of peaks in the convolved maps, considered as 
potential sources. In Scott et al. (2002), peaks identified at 
> 3 mjy were considered as possible sources and included in 
the source extraction matrix, whereas in these simulations 
all positive peaks were included in the maximum likelihood 
fit. 

The extracted sources were each identified with the 

brightest input source, located within 8 arcsec of the re- 
trieved peak position. Regions of uniform or non-uniform 



noise were assigned according to whether the input position 
lay within the uniform noise cutouts. This raises the pos- 
sibility of a source located very close to the uniform/non- 
uniform boundary being input and assigned one noise area, 
but extracted a few arcseconds away under a different noise 
classification. In these circumstances, both locations were 
allocated the input position classification so as not to un- 
derestimate the completeness. 

The tabulated results also reflect a broad range of flux 
density thresholds, some of which arc not of particular in- 
terest to every field, for example one would not expect to 
recover a 2 mJy source in a field where the rms noise lev- 
els are > 2 mJy. These values were included to allow trends 
with flux density in the various quantified properties to be 
identified. The tabulated results of these simulations may be 
found in Appendix B. 

The results of the integral completeness analyses are 
given in Tables Bl to B16, for the uniform noise regions 
of each individual field, and the non-uniform noise regions 
of the "8mJy Survey" fields. Flux density thresholds of 
2—10 mJy at 2 mJy intervals, and signal-to-noise thresh- 
olds of 1.50 — 4.00(7 at 0.50cr intervals were considered. The 
quoted errors are the la Poisson error on the number of 
sources input to the fields. Comparing the values given in Ta- 
bles Bl and B3 with B2 and B4, there is a marked contrast in 
the fraction of sources successfully retrieved in the uniform 
and non-uniform noise regions. For example, at Ssso > 8 mJy 
and a significance of 3.50(t the uniform noise regions of the 
Lockman Hole East and ELAIS N2 are 70 — 75% complete, 
whereas the non-uniform noise regions are only 10 — 15% 
complete. As the fiux density threshold drops to reach the 
faint limit at which significant (> S.OOcr) sources can still be 
detected in the uniform noise regions of the respective sur- 
veys (corresponding to an integral completeness of 65—75%), 
there is a drop of 5—10% in the fraction of sources recovered 
for every increase of 0.50a in the signal-to-noise threshold. 
The estimated completeness percentages in the small deep 
surveys (such as the lOh and SSA13 deep fields) should be 
considered as lower limits due to their small area in relation 
to the beam size. The undersampling in the jiggle pattern af- 
fects the data taken up to a beam-width into the field. These 
elevated noise levels in turn can mask the identification of a 
potential source located up to a further half a beam-width 
into the map, despite that region being fully sampled and 
hence uniform noise. A greater proportion of the field area 
in a small SCUBA map is affected by this problem than in a 
wider-area field and may artificially increase the number of 
bright sources which fail to be recovered. The error bars on 
the integral completeness measurements at bright flux den- 
sities (8 or lOmJy) are also generally large (±20%) in the 
smaller fields due to a low source density and hence fewer 
bright sources being entered into the pencil-beam maps. 

Tables B17 - B32 give the percentage count correction 
factors, above a specific flux density level, and at a given 
signal-to-noise ratio threshold, for the uniform noise regions 
of each individual field, and the non-uniform noise regions of 
the "8 mJy Survey" fields. Count correction factors at 2 mJy 
intervals from Ssso > 2mJy to Ssso > lOmJy, and signifi- 
cances of better than 1.50— 4.00it at 0.50cr intervals were con- 
sidered. The first point of note is that the percentage count 
corrections in the non-uniform noise regions of the "SCUBA 
8mJy Survey" fields (Tables B18 and B20) are in most cases 
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greater than 100% for the 3.50 and 4.00 significance levels at 
all flux densities, reflecting the poor levels of completeness in 
these regions. The converse, however, is true in the regions 
of uniform noise, implying that the effects of flux-boosting 
(discussed in Subsection 4.1.2) have a stronger effect on the 
source counts than incompleteness. In the Lockman Hole 
East and ELAIS N2 flclds, the correction factors applied to 
the bright counts (8 or lOmJy) is 30 - 40% at S/N > 3.00, 
50 - 60% at S/N > 3.50, and 70 - 90% at S/N > 4.00, the 
trend in signal-to-noise ratio indicative of an increase in the 
contamination of spurious/confused sources in the raw cat- 
alogues as the significance threshold is lowered. The count 
corrections become less severe with decreasing noise levels. 
The intermediate depth and sized areas such as the 03h wide 
area field ("CUDSS"), Lockman Hole deep strip ("8mJy 
Survey"), and SSA17 and SSA22 ("Hawaii Survey") require 
a 70 — 80% correction to the raw counts at S/N > 3.50, and 
80 - 90% at S/N > 4.00. The deep single-pointing SCUBA 
maps such as the lOh and 22h fields ("CUDSS") and the 
Lockman Hole deep field and SSA13 deep field ( "Hawaii Sur- 
vey" ) require no count-correction in the 8 — 10 nijy range, 
and a 70 — 90% correction at Sgso > 4 mjy to Sgso > 6 mjy, 
for applied signal-to-noise thresholds of 3.00 or higher. At 2 
or 3 mJy, the very faintest source flux density levels accessi- 
ble in a blank fleld survey due to the confusion limit being 
reached, the count correction factor again exceeds 100%. 
This is because the extraction of sources becomes less com- 
plete as confusion worsens, and begins to offset and even 
exceed the effect of flux-boosting which dominated at the 
brighter end of the counts. These simulations show that the 
raw source counts will be an overestimate of the true source 
counts at the faint end, and this will be readdressed in Sec- 
tion 6 where the 850 fjxn source counts are considered in 
greater detail. 

The final property of the survey data investigated by 
these simulations is the relationship of the output-to-input 
fiux densities of the sources, with signal-to-noise ratio. The 
results for the uniform noise regions of each individual field, 
and the non-uniform noise regions of the "8 mJy Survey" 
fields are given in Tables B33 to B48, for significances in the 
range > 1.50cr to > 4.00a at 0.50a intervals. The first quan- 
tity to be considered was the fraction of "sources" recov- 
ered above a specific signal-to-noise threshold which could 
be attributed to noise only, by running the source extrac- 
tion algorithm purely on the residual signal maps with no 
background counts added in (Column 2 in the tables). Each 
source recovered from the simulated images was then classi- 
fied according to the relation between the output and iden- 
tified input fiux density. The classes were: 

1) Fainter. The retrieved fiux density was fainter than the 
input source with which it had been identified (Sin > Sout). 

2) Within error bars. The input fiux density lay within the 
farms error bars of the retrieved value (Sout — errout < Sin < 
Sout + errout). 

3) Boosted. The input fiux density was less than the lower 
error boundary on the output value, but was still within a 
factor of 2 of the measured flux density (Sout/2 < Sin < 
Sout — errout). 

4) Spurious / confused. The fltted flux density to the peak 
could not be identified with a source in the input cata- 
logue, located within 8 arcsec and a factor 2 in brightness 



(Sin < Sout/2). 

The percentage of sources classified as (1), (2), (3) and (4) 
are given in Columns 3, 4, 5 and 6 respectively. As a point 
of note, the peaks identified in the residual signal image do 
not just contribute to the confused / spurious fraction, but 
may affect any of the classifications of source to some ex- 
tent. In all cases, the fraction of sources which are recovered 
at a fainter fiux density than they were input is < 10%. 
In the uniform noise regions, ~ 65 — 70% of the sources 
recovered with S/N > 4.00 may be identified as boosted 
or within the error bars. In the wider-area shallower sur- 
veys, this number drops to 55 — 65% for S/N > 3.50 and 
~ 40-55% for S/N > 3.00. The decline with signal-to-noise 
ratio is less marked in the fields with lower rms noise levels, 
and in fact remains approximately constant at the 65 — 70% 
level in the deep single-pointing SCUBA fields. The num- 
ber of sources categorised as boosted is approximately the 
same as the number for which the identified input object 
fell within the extracted farms error bars. The confused / 
spurious fraction of sources in the non-uniform regions of 
the "8mJy Survey" fields are markedly higher than their 
uniform region counterparts, even cutting at high signal-to- 
noise levels. The Lockman Hole East is the worst of the two, 
with ~ 90% of the recovered "sources" unidentified with an 
input source at least 1/2 as bright, even at > 4.00(7. The 
ELAIS N2 field is not quite as severe, but still ~ 40% of the 
recovered > 4.00a peaks fall into the spurious / confused 
category, rising to ~ 90% at > 3.00a. This affect arises due 
to redundancy i.e. the number of times a region of sky has 
been observed. The central regions of the '8 mJy Survey' 
fields have a higher redundancy and hence according to the 
"Central Limit Theorem" the noise in those areas which 
have received the full integration time will be more Gaus- 
sian. The noise in the undersampled perimeter regions of 
the wide area maps which have not received the full integra- 
tion time, however, will be less Gaussian and lead to a lower 
source detection reliability for a given signal-to-noise thresh- 
old in these areas. This casts severe doubt on the reality of 
any > 3.00a objects identified in the high noise regions near 
the edge of these maps. 

The simulations imply that up to 20 — 30% of the 

> 4.00a peaks may the result of confusion or noise, increas- 
ing to 30-40% at > 3.50a. and - 30-60% at > 3.00a. The 
fields with the highest noise levels (Lockman Hole wide area 
field, ELAIS N2, and the SSA13 wide area field) show the 
highest levels of contamination and the deepest surveys the 
least at both S/N > 3.00 and S/N > 3.50, however, there is 
no obvious trend of spurious / confused fraction with noise 
at > 4.00a. We reiterate that without particularly tight con- 
straints on the number density of the faint SCUBA popular- 
tion, and without knowing the clustering properties, these 
quantities should be considered a rough guideline only. 

The opposing effects of increasing the completeness of 
a catalogue by dropping to lower signal-to-noise thresholds, 
while at the same time also introducing a larger fraction of 
spurious / confused sources suggests that setting a cut-off of 

> 3.50a is a good compromise for selecting SCUBA sources 
to follow-up. Unless otherwise stated, subsequent analyses 
in Sections 6 and 7 are based on the S/N > 3.50 lists. 

Simulations of a similar nature to these have been car- 
ried out by Eales et al. (2000), who used their raw 14h field 
counts to produce 5 simulations of this field. They reported 
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an integral completeness of ~ 90% at the Ssso > 3mJy 
based on a 3a catalogue, which is rather higher than that 
implied from this analysis (integral completeness ~ 60% at 
Ssso > 4mJy based on objects retrieved at > S.OOcr). This 
discrepancy may be explained as a combination of 2 effects. 
The first is in the source counts model used. Eales et al. 
(2000) used the raw 850 /irri counts from the 14h field to 
create the simulated images. They also, however, reported 
the flux-boosting effect when comparing their output and 
input catalogues, quoting a median factor of 1.44, albeit 
with a large scatter about this value. This means that the 
input source counts will on average have overestimated the 
real input counts by a factor of 1.44, and this may in turn 
have affected the completeness estimate. The second, and 
probably dominant point, is that the two source extraction 
mechanisms differ, leading to different source lists and dif- 
ferent significances for those sources common to both. Ap- 
poximately 50% of the sources reported as > 3a by Eales 
et al., were recovered at < S.OOu using the simultaneous 
maximum-likelihood fit described in Section 3. If one in- 
stead compares with the completeness values at > 2.00(7 
or > 2.50(7 from these simulations, the completenss level is 
closer to 70%, which is slightly more inkeeping with that of 
Eales et al. (2000). 

Hughes & Gaztanaga (2000) have also performed more 
sophisticated Monte Carlo simulations of SCUBA surveys, 
as part of a broader investigation into the effects of confusion 
and noise on submillimetre surveys in general, for both exist- 
ing and planned instruments and telescope facilities. They 
incorporated a clustering component into the simulations 
by employing an N-body simulation with 200^ particles in 
a 600h~'^ Mpc box, produced with the same matter power 
spectrum as that measured for the APM Survey galaxies 
(Gaztanaga & Baugh, 1998, and references therein), which 
was then replicated to cover the total extent of the sur- 
vey. The local IRAS 60 /im luminosity function (Saunders 
et al. 1990) was interpolated to longer wavelengths, and a 
model of pure luminosity evolution of the form (1 -I- z)^ for 
< z < 2.2, (1 + 2.2)'' for 2.2 < z < 6 and an exponen- 
tial cutoff for z > 6 was assumed to account for the star 
formation histories. An Arp 220 SED was assumed through- 
out. The simulations of Hughes & Gaztanaga (2000) find the 
same effects of confusion and noise on the expected number 
counts as those simulations presented here. Their mock im- 
ages of comparable size and depth to the "SCUBA 8 mjy 
Survey" fields (Lockman Hole East wide area and ELAIS 
N2) imply that a count-correction factor of ~ 50% for a 
S/N > 3 catalogue at Ssso > 8m,Iy is required, which is 
in keeping with the quantities presented in Tables B17 and 
B19, again implying a large quantity of either boosted or 
spurious / confused sources. This affects the brighter flux- 
limited surveys more, due to the steep decline in galaxy 
counts and therefore increased contribution of the measured 
counts due to noise. 

5 COMPARISON WITH EXISTING SOURCE 
LISTS 

5.1 The SCUBA 8 mJy Survey 

The "SCUBA 8mJy Survey" (Lutz et al. 2001, Scott et al. 
2002, Fox et al. 2002, Ivison et al. 2002, Almaini et al. 2003) 




5 0-5 
ARC MINUTES 

CENTRE: R.A. 16 36 48.85 DEC +41 01 48.5 J2000 



Figure 8. The 850 iim image of the ELAIS N2 field, smoothed 
with a beam-size Gaussian (14.5 aresee FWHM). Tlie numbered 
circles highlight those sources found at a significance of > 3.00. 
The labelling corresponds to the numbers in Table 5. 




5 0-5 
ARC MINUTES 

CENTRE: R.A. 10 52 8.82 DEC +57 21 33.8 J2000 



Figure 9. The 850 /jja image of the Lockman Hole East field, 
smoothed with a beam-size Gaussian (14.5 arcsec FWHM). The 
numbered circles highlight those sources found at a significance 
of > 3.00. The labelling corresponds to the numbers in Table 6. 
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RA 


DEC 


Ssso 


S/N 


Noise 


Previous 


Prev. 


Sep. 




(J2000) 


(J2000) 


/mjy 




Region 


Reference 


S/N 


/ arcsec 


01 


16:37:04.29 


+41:05:30.9 


11.1 ± 1.7 


8.54 


central 


S02 (N2.01) 


8.59 


0.9 


02 


"1/^ r'o /^o 

16:36:58.62 


+41:05:24.9 


11.0 ± 1.9 


6.92 


central 


C^riri /ATri rkri\ 

b02 (N2.02) 


6.27 


1.3 


03 


16:36:58.18 


+41:04:39.9 


9.5 ± 1.8 


6.02 


central 


S02 (N2.03) 


5.86 


2.1 


04 


16:36:50.04 


+40:57:33.0 


8.3 + 1.8 


5.20 


central 


S02 (N2.04) 


5.18 


0.5 


05 


16:36:39.36 


+40:56:38.9 


9.1 + 2.4 


4.14 


central 


S02 (N2.07) 


4.07 


1.0 


06 


16:37:04.25 


+40:55:44.9 


9.2 + 2.4 


4.12 


central 


S02 (N2.06) 


4.13 


0.6 


07 


16:36:35.66 


+40:55:56.9 


8.4 + 2.2 


4.10 


central 


S02 (N2.05) 


4.16 


1.3 


08 


16:37:02.50 


+41:01:22.9 


6.1 + 1.7 


4.00 


central 


S02 (N2.12) 


3.65 


0.1 


09 


16:37:07.97 


+40:59:30.9 


6.0 ± 1.6 


3.94 


central 








10 


16:36:51.37 


+41:05:06.0 


6.5 ± 1.8 


3.90 


central 


S02 (N2.17) 


3.50 


0.3 


11 


16:36:22.41 


+40:57:04.8 


9.2 ± 2.6 


3.84 


edge 


S02 (N2.09) 


3.76 


0.2 


12 


16:37:07.46 


+41:02:36.9 


6.2 ± 1.7 


3.84 


central 


S02 {N2.30) 


3.13 


0.5 


13 


16:36:49.34 


+41:04:17.0 


7.7 ± 2.2 


3.73 


central 


S02 (N2.20) 


3.48 


0.7 


14 


16:36:58.78 


+40:57:32.9 


5.0 ± 1.4 


3.71 


central 


S02 (N2.08) 


3.82 


0.2 


15 


16:36:44.48 


+40:58:38.0 


7.3 ±2.1 


3.66 


central 


S02 (N2.11) 


3.67 


0.2 


16 


16:36:48.81 


+40:55:54.0 


5.5 ± 1.6 


3.65 


central 


b02 (JN2.10j 


3.69 


0.1 


17 


16:36:31.25 


+40:55:46.9 


6.4 + 1.9 


3.64 


central 


S02 (N2.13) 


3.56 


0.6 


18 


16:37:04.27 


+41:01:06.9 


6.4 + 1.9 


3.62 


central 








19 


16:36:19.68 


+40:56:22.7 


11.2 + 3.3 


3.55 


edge 


S02 (N2.14) 


3.55 


0.4 


20 


16:37:10.10 


+41:00:16.8 


5.1 ± 1.5 


3.54 


central 


S02 (N2.15) 


3.52 


1.1 


21 


16:36:59.41 


+40:59:57.9 


8.1 ± 2.5 


3.50 


central 








22 


16:37:19.47 


+41:01:37.7 


12.4 ± 3.8 


3.46 


edge 


S02 (N2.23) 


3.45 


0.5 


23 


16:36:27.90 


+40:54:03.9 


13.2 ± 4.1 


3.42 


edge 


S02 (N2.22) 


3.46 


0.1 


24 


16:36:48.27 


+41:03:52.0 


7.0 ± 2.2 


3.25 


central 


S02 (N2.27) 


3.25 


0.3 


25 


16:36:34.50 


1 Ar\ f^'T r\ 

+40:57:23.9 


6.4 ± 2.0 


3.30 


central 








26 


16:36:57.11 


+40:59:36.0 


7.6 ± 2.4 


3.28 


central 








27 


16:37:12.23 


+41:02:57.8 


5.1 ± 1.6 


3.27 


central 








28 


16:37:10.54 


+41:00:48.8 


4.6 ± 1.5 


3.26 


central 


S02 (N2.21) 


3.47 


0.5 


29 


16:36:26.89 


+41:02:22.8 


10.9 ± 3.5 


3.24 


edge 


b02 (JN2.24) 


3.43 


0.2 


30 


16:36:33.96 


+41:01:36.9 


11.8 ± 3.8 


3.22 


edge 








31 


16:36:39.79 


+41:00:33.9 


5.9 ± 2.0 


3.14 


central 


b02 (JN2.32) 


3.07 


0.2 


32 


16:36:36.51 


+41:05:17.9 


11.8 ± 4.0 


3.13 


edge 








33 


16:36:24.07 


+40:59:34.8 


10.4 ± 3.5 


3.12 


central 


OrkO /TVTO 

b02 (N2.29) 


3.14 


0.4 


34 


16:36:53.66 


+41:00:49.0 


7.2 ± 2.4 


3.11 


central 








35 


16:36:44.83 


+40:56:51.0 


5.7 ± 1.9 


3.10 


central 


S02 (N2.35) 


3.02 


0.3 


36 


16:36:45.44 


+41:04:52.0 


5.6 ± 1.9 


3.08 


central 








37 


16:36:35.90 


+41:01:37.9 


9.3 ± 3.2 


3.06 


central 


S02 (N2.19) 


3.49 


0.1 


38 


16:36:32.02 


+41:00:04.9 


9.2 ± 3.1 


3.06 


central 


S02 (N2.26) 


3.26 


0.2 


39 


16:36:53.14 


+41:03:46.0 


5.7 ± 2.0 


3.05 


central 








40 


16:37:08.62 


+41:04:48.9 


4.3 ± 1.5 


3.05 


central 








41 


16:.36:25.57 


+41:00:35.8 


7.8 ± 2.7 


3.04 


central 








42 


16:36:52.87 


+41:02:52.0 


3.9 ± 1.3 


3.03 


central 


S02 (N2.36) 


3.00 


0.3 


43 


16:36:27.00 


+40:58:14.8 


6.8 ± 2.4 


3.02 


central 


S02 (N2.34) 


3.05 


0.2 


44 


16:37:01.81 


+41:06:22.9 


6.1 ± 2.1 


3.01 


edge 








45 


16:37:09.13 


+41:01:59.9 


4.3 + 1.5 


3.01 


central 










lo:jo:5U.40 


+40:58:54- 


4-6 ± 1 .6 


2.90 


central 


b02 (iN2.33j 


3.06 


0.2 




16:36:28.21 


+41:01:41.9 


6.8 + 2.5 


2.88 


central 


S02 (N2.31) 


3.07 


1.5 




16:36:47.21 


+41:08:48.0 


4.0 ±1.9 


2.14 


central 


S02 (N2.28) 


3.24 


0.1 




16:36:18.34 


+40:59:11.7 


10.5 ± 5.2 


2.07 


edge 


S02 (N2.25) 


3.37 


0.5 




16:36:51.99 


+41:05:54.0 


6.8 ±3.6 


1.93 


edge 


S02 (N2.16) 


3.51 


4.0 




16:36:11.36 


+40:59:25.6 


13.3 ± 11.8 


1.13 


edge 


S02 (N2.18) 


3.49 


0.6 



Table 5. 850 /im source list for the ELAIS N2 field of the "SCUBA 8mjy Survey". Sources are marked in Fig 8. 



is divided between two fields; the Lockman Hole East (cen- 
tred at RA 10:52:08.82, DEC +57:21:33.8) and ELAIS N2 
(centred at RA 16:36:48.85, DEC +41:01:48.5). These fields 
both lie in regions of low galactic cirrus emission (Schlegel, 
Finkbeiner & Davis 1998), and have a vast quantity of multi- 
wavelength data available for follow-up studies. They were 
also selected to coincide with deep Infrared Space Observa- 



tory (ISO) surveys at 6.7, 15, 90 and 175 /im (Lockman Hole 
East - Elbaz et al. 1999, Kawara et al. 1998; ELAIS N2 - 
Oliver et al. 2000, Serjeant et al. 2000, Efstathiou et al. 2000, 
Rowan- Robinson et al. 2004). This section describes the final 
results of the "SCUBA 8 mJy Survey" , currently the largest 
of the deep 850 /im blank field surveys to be completed. The 
Lockman Hole East map covers ~ 140 sq. arcmin of sky of 
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which ~ 120 sq. arcmin of the map have an rms noise level 
of <785() — 2.7 mjy/bcam and the remaining 20 sq. arcmin 
form a deeper strip in the centre of the image to an rms 
noise level of csso — 1-8 mjy/beam. The ELAIS N2 field is 
also ~ 140 sq. arcmin in size and reaches an rms depth of 
IT850 — 2.2 mJy/beam. 

More up-to-date 850 /im source lists for the "SCUBA 
8mJy Survey" than those given in Scott et al. (2002) are 
presented here. There are two reasons for the differences be- 
tween this list and that given in Scott et al (2002). Firstly, 
~ 20 sq. arcmin of additional data were taken after the pub- 
lication of the original paper, most of it in the ELAIS N2 
field, and this has now been added into the survey maps. Sec- 
ondly, rather than doing a fit to only those peaks identified 
at > 3mJy in the smoothed images, the simultaneous max- 
imum likelihood source extraction algorithm was applied to 
every peak in the maps. 

Column 1 gives the source number in order of de- 
creasing signal-to-noise ratio as derived from the new IDL- 
reduction and simultaneous maximum-likelihood source ex- 
traction method, and corresponds to the labelling of the 
circled sources in Figs. 8 and 9. Columns 2 and 3 give the 
right ascension and declination of the source in J2000 coordi- 
nates. Column 4 gives the simultaneously fitted 850 yitm flux 
densities of the sources. The error includes a 10% calibration 
error combined in quadrature. Column 5 gives the measured 
signal-to-noise ratio of the source from the simultaneously- 
fitted model. Column 6 defines the noise region in which 
the source was found; 'deep' corresponds to the deep linear 
strip at the centre of the Lockman Hole East field, 'central' 
corresponds to the parts of the map which have seen the 
full integration time (outside of the deep area), and 'edge' 
corresponds to the rather noisier regions near the perime- 
ter which have not seen the full integration time. Column 

7 gives any previous reference to the 850 /im source. Refer- 
ence S02 is an abbreviation for Scott et al. (2002). Column 

8 gives the previously recorded signal-to-noise ratio where 
applicable, and Column 9 gives the distance between the 
listed and previously referenced positions. The table listings 
given in italics correspond to previously referenced sources 
with S/N > 3.00, which did not meet this criterion in this 
analysis. 

Table 5 shows that the top 7 sources originally detected 
at > 4.00(7 in the ELAIS N2 field (Scott et al. 2002) remain 
robustly detected at this signal to noise threshold. One can 
also see that the top 15 sources presented in Scott et al. 
(2002) are confirmed with S/N > 3.50. Similarly, Table 6 
shows that the top 10 sources presented in Scott et al. for 
the Lockman Hole East field are confirmed with S/N > 4.00 
and the original top 15 sources have all been detected at 
> 3.50cT. 

5.2 The Canada-UK Deep Submillimetre Survey 

The "Canada UK Deep Submilhmetre Survey (CUDSS)" 
(Eales et al. 1999, Lilly et al. 1999, Gear et al. 2000, 
Eales et al. 2000, Webb et al. 2003a), covers a total of 
~ 130 sq. arcmin over 4 regions of sky, selected to coin- 
cide with areas observed in the "Canada-France Redshift 
Survey (CFRS)" (Lilly et al. 1995). The 03-Hour field is 
composed of a deep pencil beam area (~ 8 sq. arcmin in 
size, (T850 = l.lmjy/beam; Eales et al. 1999, Lilly et al. 
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Figure 10. The 850^4111 image of the CUDSS 03-Hour field, 
smootlicd with a beam-size Gaussian (14.5 arcsec FWHM). Tlie 
numbered circles higliligtit those sources found at a significance 
of > 3.00. The labelling corresponds to the numbers in Table 7. 
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Figure 11. The 850 /xm image of the CUDSS 10-Hour field, 
smoothed with a beam-size Gaussian (14.5 arcsec FWHM). The 
numbered circles highlight those sources found at a significance 
of > 3.00. The labelling corresponds to the numbers in Table 8. 
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RA 


DEC 


Ssso 


S/N 


Noise 


Previous 


Prev. 


Sep. 




(J2000) 


(J2000) 






Region 


Reference 


S/N 


/ arcsec 


01 


10:52:01.33 


+57:24:43.3 


9.6+ 1.6 


7.68 


deep 


S02 


(LH.Ol) 


8.10 


0.6 


02 


10:52:38.21 


+57:24:35.1 


11.0 + 2.3 


5.33 


central 


S02 


(LH.02) 


5.22 


0.9 


03 


10:51:58.39 


+57:18:00.3 


7.9 + 1.7 


5.31 


deep 


S02 


(LH.03) 


5.06 


1.0 


04 


10:52:04.05 


+57:25:29.3 


7.8 + 1.8 


4.88 


deep 


S02 


(LH.04) 


5.03 


1.4 


05 


10:52:30.39 


+57:22:13.2 


10.8 ± 2.6 


4.58 


central 


S02 


(LH.06) 


4.50 


2.1 


06 


10:52:22.71 


+57:19:32.3 


14.0 ± 3.4 


4.52 


central 


S02 


(LH.09) 


4.20 


0.3 


07 


lU:oi:oi.o4 


+O7:2o:o5.2 


8.0 + 2.0 


4.40 


deep 


S02 


(LH.07) 


4.50 


n A 
U.4 


08 


10:51:59.60 


+57:24:21.3 


4.8 ± 1.2 


4.26 


deep 


S02 


(LH.08) 


4.38 


3.3 


09 


10:51:59.26 


+57:17:18.3 


8.0 ± 2.1 


4.12 


deep 


S02 


(LH.05) 


4.57 


0.4 


10 


10:51:42.39 


+57:24:45.1 


12.0 ± 3.2 


4.04 


central 


S02 


(LH.IO) 


4.18 


0.1 


11 


10:51:53.82 


+57:18:47.3 


6.7+ 1.8 


4.01 


deep 


S02 


(LH.27) 


3.38 


0.3 


12 


10:52:16.78 


+57:19:23.3 


10.1 ± 2.8 


3.80 


central 


S02 


(LH.17) 


3.55 


0.3 


13 


10:51:33.57 


+57:26:41.0 


9.5 + 2.7 


3.73 


central 


S02 


(LH.13) 


3.69 


0.2 


14 


10:52:07.77 


+57:19:07.3 


5.6+ 1.6 


3.71 


deep 


S02 


(LH.12) 


4.01 


0.6 


15 


lU:51:oU.4o 


+0 r :20:o ( .9 


12.0 + 3.5 


3.70 


edge 


S02 


(LH.ll) 


4.01 


1.1 


16 


10:52:36.37 


+57:25:15.1 


5.7+ 1.7 


3.67 


central 










17 


10:52:04.30 


+57:27:01.3 


9.2 ± 2.7 


3.64 


central 


S02 


(LH.14) 


3.61 


2.3 


18 


10:52:05.67 


+57:20:53.3 


4.6 ± 1.4 


3.52 


deep 


S02 


(LH.22) 


3.61 


0.4 


19 


10:52:24.58 


+57:21:19.3 


11.7 + 3.5 


3.50 


central 


S02 


(LH.15) 


3.60 


0.3 


20 


10:52:27.18 


+57:22:21.2 


12.1 + 3.7 


3.40 


central 










21 


10:52:37.67 


+57:20:30.1 


9.8 + 3.1 


3.38 


edge 


S02 


(LH.20) 


3.51 


0.3 


22 


10:51:46.97 


+57:24:51.2 


7.1 + 2.2 


3.38 


central 


S02 


(LH.23) 


3.48 


0.3 


23 


10:52:28.07 


+57:25:09.2 


5.8 + 1.8 


3.37 


central 


S02 


(LH.16) 


3.56 


10.4 


24 


10:52:03.94 


+57:20:07.3 


4.0 + 1.3 


3.36 


deep 


S02 


(LH.31) 


3.24 


0.4 


25 


10:51:48.36 


+57:21:48.2 


12.1 + 3.8 


3.34 


central 










26 


10:52:09.87 


+57:20:40.3 


9.0 + 2.8 


3.33 


central 


S02 


(LH.34) 


3.16 


0.4 


27 


10:52:34.57 


+57:20:02.1 


10.0 + 3.2 


3.29 


central 


S02 


(LH.28) 


3.31 


0.3 


28 


10:51:42.89 


+57:24:12.1 


11.0 + 3.6 


3.26 


central 


S02 


(LH.24) 


3.47 


0.1 


29 


10:52:01.71 


+57:19:16.3 


4.0 + 1.3 


3.23 


deep 


S02 


(LH.21) 


3.50 


0.3 


30 


10:51:59.99 


+57:20:39.3 


4.2 + 1.4 


3.21 


deep 


S02 


(LH.32) 


3.22 


0.3 


31 


10:52:27.28 


+57:19:06.2 


7.7 + 2.5 


3.20 


central 


S02 


(LH.26) 


3.31 


0.3 


32 


10:52:34.51 


+57:25:34.1 


4.1 ± 1.4 


3.18 


central 










33 


10:52:10.86 


+57:24:13.3 


8.4 + 2.8 


3.17 


central 










34 


10:51:33.81 


+57:19:29.0 


7.8 + 2.6 


3.12 


central 


S02 


(LH.33) 


3.20 


0.1 


35 


10:51:55.77 


+57:23:12.3 


4.0 + 1.3 


3.12 


deep 




I^Ijxi. lo ) 


3.55 


0.6 


ou 


^c\■^^ ■9'^ 9Q 

lU.u-L.zo.zy 


_i_c:7.9n.on o 
T^O ( . ZU.OO.O 


1 q Q -I- 4 7 


o.uy 


edge 










37 


10:52:29.57 


+57:26:20.2 


4.7+1.6 


3.09 


central 






3.54 


1.6 


38 


10:52:41.14 


+57:21:47.1 


14.3 + 4.9 


3.06 


edge 










39 


10:52:03.45 


+57:16:54.3 


7.9 + 2.7 


3.06 


central 


S02 


(LH.36) 


3.00 


0.5 


40 


10:51:52.56 


+57:22:29.2 


4.4 + 1.5 


3.06 


deep 












10:52:42.21 


+57:18:28.0 


10.0 ± 3.5 


2.99 


edge 


S02 


{LH.30) 


3.25 


0.1 




10:52:16.43 


+57:25:04.3 


5.0 ±1.8 


2.95 


central 


S02 


(LH.29) 


3.30 


2.7 




10:52:36.03 


+57:18:20.1 


10.5 ± 3.8 


2.91 


edge 


S02 


(LH.25) 


3.46 


0.3 




10:51:57.61 


+57:26:03.3 


6.5 ± 2.3 


2.90 


central 


S02 


(LH.35) 


3.02 


0.3 



Table 6. 850 /im source list for the Lockman Hole East field of the 

1999), cmboddod in a widor-aroa, shallower map covering 
an additional 55 sq. arcmin with a typical rms noise level of 
(T850 = 1-8 mjy/beam (Webb et al. 2003a). The 14-Hour field 
(Bales et al. 2000) is similar in size, the uniform noise re- 
gion covering approximately 57 sq. arcmin, but to a slightly 
deeper uniform noise level of aggo = 1.5mJy/beam. The 
10-Hour and 22-Hour fields (Bales et al. 1999, Lilly et al. 
1999) are small in area, each having a uniform noise region 
of ~ 7 sq. arcmin, with rms noise levels of ctsbo = 1.3 and 
1.5mJy/beam respectively. 

The data were originally reduced independently by the 
Cardiff and Toronto groups, using the standard SURF proce- 
dures. They made additional attempts to improve the qual- 



"SCUBA 8 mjy Survey" . Sources are marked in Fig 9. 

ity of the final map, firstly by allowing the residual sky 
removal to be a linear function of position (i.e. a planar 
fit was applied rather than a D.C. oflfset), and secondly by 
examining the Fourier-transform of each bolometer's mea- 
sured signal to search for non-white noise profiles (although 
Bales et al. (2000) reported this produced negligible im- 
provements to the final regridded images). The chop throw 
in all cases was fixed at 30 arcsec in an east-west direction as 
observed on the sky. Source extraction was carried out by 
convolving a normalised template of the full beam-profile, 
constructed from the many observations of Uranus taken 
thoughout the lifetime of this survey, with the raw survey 
maps. The "CUDSS" team used a method of noise modelling 



© 0000 RAS, MNRAS 000, 000-000 



16 S.E. Scott et al. 





RA 


DEC 


Ssso 


S/N 


Noise 


Previous 


Prev. 


Sep. 




(J2000) 


(J2000) 


/mjy 




Region 


Reference 


S/N 


/ arcsec 


01 


03:02 


43.84 


+00:09:52.6 


7.0 + 1.4 


5.90 


central 


W03 (03h.l9) 


3.2 


1.8 


02 


03:02 


36.04 


+00:08:16.6 


3.4 + 0.8 


4.89 


deep 


W03 (03h.06) 


5.4 


1.3 


03 


03:02 


42.84 


+00:07:57.6 


3.6 + 0.9 


4.61 


deep 


W03 (03h.O2) 


6.7 


3.9 


04 


03:02 


40.77 


+00:09:20.6 


6.9 + 1.7 


4.51 


central 








05 


03:02 


31.17 


+00:08:18.6 


4.8 + 1.2 


4.31 


central 


W03 (03h.03) 


6.4 


5.1 


06 


03:02 


56.57 


+00:08:08.6 


7.0 + 1.8 


4.21 


central 


W03 (03h.24) 


3.0 


3.5 


07 


03:02 


47.31 


+00:09:21.6 


/I T _1_ 1 Q 

4.7 ± 1.0 


4.18 


central 








08 


03:02 


53.51 


+00:07:52.6 


9.0 ± 2.3 


4.16 


central 








09 


03:02 


44.51 


+00:06:53.6 


5.6 ± 1.5 


4.14 


central 


W03 (03h.04) 


6.2 


2.2 


10 


03:02 


44.31 


+00:08:16.6 


4.0 ± 1.0 


4.14 


central 


W03 (03h.05) 


5.8 


5.3 


11 


03:02 


44.64 


+00:06:37.6 


6.0 ± 1.6 


4.03 


central 


W03 (03h.01) 


7.4 


3.4 


12 


03:02:29.04 


+00:09:03.6 


5.3 + 1.5 


3.87 


central 








13 


03:02:41.51 


+00:10:46. 


9.8 i 2.8 


3.81 


central 








14 


03:02:35.64 


+00:06:09.6 


7.9 ± 2.2 


3.79 


central 


W03 (03h.07) 


5.3 


2.2 


15 


03:02:27.84 


+00:06:45.6 


6.0 ± 1.8 


3.63 


central 


W03 (03h.l5)* 


3.5 


7.8 
















W03 (03h.27)* 


3.0 


13.5 


16 


03:02:53.11 


+00:09:41.6 


4.8 ± 1.4 


3.50 


central 


W03 (03h.20) 


3.2 


3.3 


17 


03:02 


25.31 


+00:10:17.6 


7.2 + 2.2 


3.49 


central 








18 


03:02 


40.97 


+00:06:44.6 


3.0 + 0.9 


3.48 


deep 








19 


03:02 


53.51 


+00:06:23.6 


7.3 + 2.2 


3.47 


central 


W03 (03h.23) 


3.0 


10.9 


20 


03:02 


25.97 


+00:09:06.6 


6.3 + 1.9 


3.47 


central 


W03 (03h.l4) 


3.5 


3.0 


21 


03:02 


58.17 


+00:06:07.6 


13.9 + 4.3 


3.44 


edge 








22 


03:02 


52.44 


+00:08:58.6 


4.3 + 1.4 


3.32 


central 


W03 (03h.l0) 


4.5 


1.4 


23 


03:02 


43.51 


+00:10:51.6 


9.0 + 2.9 


3.23 


central 








24 


03:02 


48.24 


+00:08:03.6 


3.9 + 1.3 


3.15 


central 








25 


03:02 


45.11 


+00:09:53.6 


4.0 + 1.3 


3.14 


central 










03:02 


52.97 


_i_nn-i 1 -91 


6 2 + 21 


o.uo 


central 


vvuo ^uon.j.±^ 




]^ I 


27 


03:02 


30.51 


+00:08:50.6 


4.8 + 1.6 


3.07 


central 








28 


03:02 


40.31 


+00:11:42.6 


6.5 + 2.2 


3.04 


central 








29 


03:02 


55.37 


+00:09:49.6 


5.3 + 1.8 


3.03 


central 








30 


03:02 


35.71 


+00:12:05.6 


7.2 + 2.5 


3.00 


central 










03:02 


35.91 


+00:09:57.6 


3.7 ± 1.3 


2.97 


central 


W03 (03h.l3) 


3.8 


4.4 




03:02 


38.77 


+00:10:27.6 


5.2 ± 1.9 


2.80 


central 


W03 (03h.l2) 


4.0 


1.9 




03:02 


26.24 


+00:06:18.6 


4.8 ± 1.9 


2.67 


central 


W03 (03h.08) 


5.0 


2.3 




03:02 


26.11 


+00:08:17.6 


3.8 ± 1.5 


2.66 


central 


W03 (03h.21) 


3.1 


3.4 




03:02 


32.77 


+00:10:20.6 


4.0 ±1.7 


2.51 


central 


W03 (03h.l8) 


3.3 


5.8 




03:02 


31.37 


+00:10:33.6 


4.5 + 2.0 


2.33 


central 


W03 (0311.17) 


3.4 


5.2 




03:02 


34.97 


-f-UU. Utf. 10. 


J .4 ±1.0 


Z.oo 


central 


VVUcS yjoii.zo) 


o.U 


A ^ 
4.0 




03:02 


28.31 


+00:10:19.6 


3.5 ±1.7 


2.09 


central 


W03 (03h.09) 


4.6 


8.9 




03:02 


35.44 


+00:08:51.6 


2.3 ±1.1 


2.08 


central 


W03 (03h.l6) 


3.4 


9.5 




03:02 


39.17 


+00:06:14.6 


2.2 ±1.3 


1.77 


central 


W03 {03h.22) 


3.1 


12.5 




03:02 


38.11 


+00:11:10.6 


2.9 + 2.4 


1.22 


central 


W03 (03h.25) 


3.0 


8.2 



Table 7. 850 /im source list for the 03h field of the "Canada-UK Deep Submillimetre Survey" . Sources are marked in Fig 10. 



RA 


DEC 


S850 


S/N 


Noise 


Previous 


Prev. 


Sep. 


(J2000) 


(J2000) 


/mJy 




Region 


Reference 


S/N 


/ arcsec 


01 10:00:36.86 


+25:14:56.9 


3.3 + 1.1 


3.12 


central 


E99 (lOh.B)* 


5.0 


3.4 












E99 (lOh.C)* 


3.8 


15.9 












E99 (lOh.D)* 


3.6 


15.0 


10:00:38.12 


+25:14:51.9 


2.8 ±1.0 


2.79 


central 


E99 (lOh.A) 


6.1 


2.1 



Table 8. 850 /im source list for the lOh field of the "Canada-UK Deep Submillimetre Survey". Sources are marked in Fig 11. 
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RA 


DEC 


Ssso 


S/N 


Noise 


Previous 


Prev. 


Sep. 




(J2000) 


(J2000) 


/mjy 




Region 


Reference 


S/N 


/ aJTcsec 


01 


14:17:40.03 


+52:29:07.0 


8.5 ± 1.4 


7.62 


central 


EOO (14h.01) 


10.1 


2.1 


02 


14:17:51.86 


+52:30:32.0 


6.0 ± 1.1 


6.01 


central 


EOO (14h.02) 


6.3 


2.1 


03 


14:18:00.61 


+52:28:20.0 


7.2 ± 1.5 


5.45 


central 


EOO (14h.03) 


5.4 


3.6 


04 


14:17:43.21 


+52:28:16.0 


5.7+1.4 


4.57 


central 


EOO (14h.04) 


5.3 


2.0 


05 


14:18:07.51 


+52:28:22.9 


5.8 + 1.4 


4.42 


central 


EOO (14h.05) 


4.5 


2.3 


06 


14:17:38.05 


+52:32:50.0 


4.9 + 1.4 


3.90 


central 








07 


14:18:09.39 


+52:32:02.9 


5.4 + 1.5 


3.76 


central 








08 


14:18:12.91 


+52:33:21.9 


9.3 + 2.7 


3.71 


edge 








09 


14:17:56.03 


+52:32:59.0 


3.7+1.1 


3.59 


central 








10 


14 


17:36.07 


+52:33:15.0 


4.1 ± 1.2 


3.47 


central 








11 


14 


17:42.22 


+52:30:31.0 


3.7 ± 1.1 


3.46 


central 


EOO (14h.l8) 


3.0 


4.5 


12 


14 


17:45.61 


+52:33:23.0 


3.9 ± 1.2 


3.37 


central 








13 


14 


17:46.93 


+52:29:20.0 


4.2 ± 1.3 


3.31 


central 








14 


14 


17:25.02 


+52:30:41.9 


9.0 ± 2.9 


3.25 


edge 


EOO (14h.l7) 


3.3 


4.5 


15 


14 


17:35.11 


+52:28:53.0 


4.5 ± 1.4 


3.25 


central 








16 


14 


18:03.26 


+52:32:29.0 


3.7 ± 1.2 


3.20 


central 








17 


14 


17:47.25 


+52:32:36.0 


3.3 ± 1.1 


3.19 


central 


EOO (14h.ll) 


3.5 


2.4 


18 


14 


17:42.66 


+52:30:04.0 


4.5 ± 1.5 


3.14 


central 








19 


14 


17:35.53 


+52:32:11.0 


3.2 + 1.1 


3.14 


central 








20 


14 


18:08.71 


+52:28:00.9 


4.4 + 1.5 


3.13 


central 


EOO (14h.09) 


4.1 


4.1 


21 


14 


18:03.68 


+52:29:33.9 


3.1 + 1.0 


3.05 


central 


EOO (14h.l0) 


3.5 


5.0 


22 


14 


17:48.13 


+52:32:51.0 


3.4 + 1.2 


3.04 


central 








23 


11 


17:13.12 


+52:32:10.0 


3.9 ± 1.3 


3.01 


(.■(~nl ral 










u 


17:41.57 


+52:28:27.0 


3.7 tt 1.3 


2.96 


central 


EOO (14n.l3) 


3.4 


3.9 




u 


18:11.68 


+52:30:05.9 


5.7 ± 2.0 


2.90 


central 


EOO (14n.l9) 


3.0 


2.5 




H 


18:12.44 


+52:29:13.9 


6.3 ± 2.4 


2.77 


central 


EOO (14h.l6) 


3.7 


6.3 




U 


18:09.28 


+52:31:01.9 


4.1 ±1.6 


2.61 


central 


EOO (14h.l4) 


3.3 


6.0 




U 


18:01.61 


+52:30:28.0 


3.3 ± 1.3 


2.56 


central 


EOO (14h.08) 


4.0 


13.8 




u 


17:56.45 


+52:29:12.0 


3.0 ±1.2 


2.50 


central 


EOO (14h.06) 


4.2 


5.2 




u 


18:05.21 


+52:28:55.9 


3.3 ± 1.4 


2.41 


central 


EOO (14h.l2) 


3.4 


0.9 




u 


17:29.53 


+52:28:17.9 


4.2 ±1.9 


2.26 


central 


EOO (14h.l5) 


3.1 


2.4 




u 


18:01.60 


+52:29:44.0 


1.8 ±1.0 


1.76 


central 


EOO (14h.07) 


3.2 


6.8 



Table 9. 850 /im source list for the 14h field of the "Canada-UK Deep Submillimetre Survey". Sources are marked in Fig 12. 





RA 


DEC 


Ssso 


S/N 


Noise 


Previous 


Prev. 


Sep. 




(J2000) 


(J2000) 


/mJy 




Region 


Reference 


S/N 


/ arcsec 


01 


22:17:59.18 


+00:17:36.9 


5.6 + 1.3 


4.62 


central 








02 


22:17:59.18 


+00:18:22.9 


4.8 + 1.3 


4.14 


central 








03 


22:17:55.58 


+00:17:36.9 


3.6 + 1.1 


3.36 


central 









Table 10. 850 /im source list for the 22h field of the "Canada-UK Deep Submillimetre Survey". Sources are marked in Fig 13. 



which is in effect quite similar to the way in which the noise 
maps wore created in the IDL reduction of the "SCUBA 
8 mJy Survey" . They began with the basic assumption that 
the noise on any bolometer was independent of the noise 
on every other bolometer, and then measured the standard 
deviation of the intensities for each bolometer in units of 
one hour (the length of each CUDSS pointing). Artificial 
data were then created by replacing the real data with the 
output of a Gaussian random-number generator with the 
same standard deviation as the real data, re-running the 
sky subtraction and clipping routines from the SURF pack- 
age to account for any non-Gaussian nature arising from 
these processes, and finally rescaling the mock data so that 
it had the same standard deviation as the real data. In total. 



1000 simulated maps were generated, each of which was con- 
volved with the beam template as were the real data. The 
final noise maps were produced by measuring the standard 
deviation of these convolved maps, pixel by pixel. 

As previously stated in Section 3, the method of con- 
volving the raw images with the normalised point spread 
function (PSF) is formally the best method of source- 
extraction, provided that the sources are all well separated 
from one another. It does, however, run into difficulties when 
dealing with partially confused sources, a problem which is 
likely to be fairly common, as inferred from the ELAIS N2 
image (Section 5.1) and the clustering analyses of Section 7. 
In the "SCUBA 8 mJy Survey" the problem of confusion was 
tackled by means of a ma^ximum-likelihood fit of the beam 
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6 4 2 -2 -4 -6 

ARC MINUTES 

CENTRE: R.A. 14 17 49.06 DEC +52 30 17.5 J2000 



Figure 12. The 850^01 image of the CUDSS 14-Hour field, 
smoothed with a beam-size Gaussian (14.5 arcsec FWHM). The 
numbered circles highlight those sources found at a significance 
of > 3.00. The labelling corresponds to the numbers in Table 9. 



(/I 

V- 




2 1 0-1-2 

ARC MINUTES 

CENTRE: R.A. 22 17 57.14 DEC +00 17 45.4 J2000 



Figure 13. The 850 //m image of the CUDSS 22-Hour field, 
smoothed with a beam-size Gaussian (14.5 arcsec FWHM). The 
numbered circles highlight those sources found at a significance 
of > 3.00. The labelling corresponds to the numbers in Table 10. 



template to all potential sources simultaneously using the 
raw data. Ealcs ot al. (2000) instead addressed the prob- 
lem of confusion by attempting a deconvolution with the 
CLEAN algorithm (Hogbom 1974). They created an initial 
list of possible sources based on the beam-convolved signal 
image divided by the Gaussian generated noise image, then 
iteratively CLEANed the raw data in boxes centred on the 
positions of the potential sources. For each source the in- 
formation from CLEAN was then used to remove all other 
possible sources from the raw image, before again carrying 
out a convolution with the beam template on this new map, 
and dividing by the noise to measure the signal-to-noise rar 
tio. 

Comparative source lists for detections with S/N > 3.00 
from the two different reductions and source extraction pro- 
cedures arc given in Tables 7-10 (see also Figs 10 - 13). 
Columns 1-5, 8 and 9 are analogous to those of Table 5. 
Column 6 defines the noise region in which the source was 
found; 'deep' corresponds to the deep pencil-beam surveys 
which constitute part of a wider-area and somewhat shal- 
lower image, 'central' corresponds to the parts of the map 
which have seen the full integration time (outside of the deep 
area), and 'edge' corresponds to the rather noisier regions 
near the perimeter which have not seen the full integration 
time. Column 7 gives any previous reference to the 850 
source. Reference E99 is an abbreviation for Eales et al. 
(1999), EOO is an abbreviation for Eales et al. (2000), and 
W03 is an abbreviation for Webb et al. (2003a). The pres- 
ence of a * indicates that a previous reduction found more 
than one source here, whereas in this reduction only one was 
found. 

As can be seen from comparing the two different re- 
duction and source extraction procedures, a majority of the 
most highly significant objects (> 4.00cr) identified in the 
initial CUDSS analyses are also recovered here, but as one 
considers detections at decreasing signal-to-noise ratios, the 
two catalogues increasingly diverge. The comparative results 
for the 03-Hour field presented in Table 7 show that the top 
6 sources detected by Webb et al. (2003a) are robustly con- 
firmed at the > 4.00(7 level while the top 7 are recovered 
above a signal-to-noise threshold of 3.50. We find all but 3 
of their 12 > 4.00cr sources at better than > 3.00cr. In ad- 
dition wc have detected a further 3 previously unpublished 
sources at better than > 4.00c7. Dropping down the list to 
lower significances, however, the resulting source catalogues 
of significant (> 3.00cr) detections are actually markedly dif- 
ferent, in particular only a handful of objects arc common 
to both lists for 3.00 < S/N < 4.00. That is not to say that 
all detections under a signal-to-noise ratio of 3.00 are spuri- 
ous, simply that there is a rapidly increasing probability of 
contamination from fake sources on decreasing the signal-to- 
noise threshold. A few of those objects identified at > 3.00cr 
by the original analysis of Webb et al. (2003a) fall only just 
short of this criterion in this new analysis adding some cred- 
ibility to their reality. In the 14-Hour field, the top 5 sources 
identified by Ealcs et al. (2000) are also securely recovered as 
the top 5 detections in this analysis with S/N > 4.00. Again, 
on dropping to lower signal-to-noise ratios the two cata- 
logues diverge. Only a further 5 of the remaining 14 objects 
detected by Eales et al. (2000) at > 3.00(t are recovered by 
this criterion in this maximum-likelihood analysis, although 
a few of these objects fall only just below this threshold. For 
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RA 


DEC 


SsBO 


S/N 


Noise 


Previous 


Prev. 


Sep. 


(J2000) 


(J2000) 


/mjy 




Region 


Reference 


S/N 


/ circsec 


01 10:34:02.05 


+57:46:27.1 


4.9 ±0.9 


6.45 


central 


B99 (LH.l) 


5.1 


2.1 


02 10:33:55.80 


-57:15:10.1 


2,0 = 0.7 


3.70 


(■('111 ral 








10:33:5542 


+57:47:38.1 


2.0 ± 0.8 


2.59 


central 


B99 (LH.2) 


3.3 


10.6 



Table 11. 850 fim source list for the Loclsman Hole Field of the "Hawaii Flanking Fields Survey". Sources are marked in Fig 14. 





RA 


DEC 


Ssso 


S/N 


Noise 


Previous 


Prev. 


Sep. 




(J2000) 


(J2000) 


/mJy 




Region 


Reference 


S/N 


/ arcsec 


01 


13:12:31.82 


+42:44:28.6 


3.6 + 0.7 


6.52 


deep 


B99 (SSA13.1) 


4.7 


3.4 


02 


13:12:13.94 


+42:37:00.7 


10.2 + 2.2 


5.19 


central 








03 


13:12:27.56 


+42:45:01.5 


2.8 + 0.6 


5.08 


deep 


B99 (SSA13.2) 


3.8 


6.0 


04 


13:12:19.93 


+42:39:30.7 


6.1 + 1.6 


4.22 


central 








05 


13:12:08.51 


+42:38:19.7 


7.0 + 1.9 


4.03 


central 








06 


13:12:25.00 


+42:39:56.7 


7.8 + 2.1 


4.01 


central 


B99 (SSA13.6) 


3.4 


1.3 



07 13:12:25.82 +42:39:38.7 

08 13:12:17.66 +42:42:51.7 

09 13:12:22.29 +42:45:00.7 



8.7 + 2.5 3.75 central 
7.4 + 2.1 3.66 central 
2.7 + 0.8 3.65 deep 



10 


13: 


12:13.94 


+42:39:49.7 


5.9 + 1.8 


3.43 


central 


11 


13: 


12:31.26 


+42:40:22.7 


8.1 + 2.5 


3.41 


central 


12 


13: 


:12:05.79 


+42:38:52.7 


10.2 + 3.2 


3.36 


central 


13 


13; 


12:14.40 


+42:43:33.7 


10.1+3.2 


3.34 


central 


14 


13; 


12:27.28 


+42:41:54.7 


6.6 + 2.1 


3.33 


central B99 (SSA13.7) 


15 


13; 


12:04.88 


+42:37:51.7 


10.2 + 3.4 


3.18 


central 


16 


13: 


12:33.71 


+42:40:22.6 


7.9 + 2.6 


3.18 


central 


17 


13: 


12:28.99 


+42:40:14.7 


7.1 + 2.4 


3.14 


central 


18 


13; 


12:11.22 


+42:42:20.7 


5.9 + 2.0 


3.14 


central 


19 


13; 


:12:27.27 


+42:38:59.7 


9.1 + 3.1 


3.03 


central 



13:12:18.66 


+42:38:25.7 


4.0 ± 1.4 


2.99 


central 


B99 (SSA13.5) 


3.3 


12.8 


13:12:25.65 


+42:43:48.5 


1.6 ± 0.6 


2.69 


deep 


B99 (SSA13.3) 


3.2 


1.6 


13:12:05.95 


+42:44:37.7 


1.5 ±3.2 


0.49 


edge 


B99 (SSA13.9) 


3.4 


10.1 


n/a 


n/a 


n/a 


n/a 


bad bol. 


B99 (SSA13.8) 


3.5 


n/a 


n/a 


n/a 


n/a 


n/a 


edge 


B99 (SSA13.4) 


3.3 


n/a 



Table 12. 850 /im source list for the SSA13 Field of the "Hawaii Flanking Fields Survey". Sources are marked in Fig 15. 





RA 


DEC 


Ssso 


S/N 


Noise 


Previous 


Prev. 


Sep. 




(J2000) 


(J2000) 


/mJy 




Region 


Reference 


S/N 


/arcsec 


01 


17:06:37.03 


+43:55:31.8 


3.2 + 1.0 


3.51 


central 


B99 (SSA17.3) 


3.7 


2.0 


02 


17:06:29.53 


+43:55:08.8 


4.0+1.2 


3.36 


central 








03 


17:06:25.08 


+43:57:40.8 


5.6 + 1.8 


3.33 


central 


B99 {SSA17.1) 


4.2 


2.0 


04 


17:06:32.86 


+43:54:05.8 


5.1 + 1.7 


3.16 


central 


B99 (SSA17.4) 


3.6 


3.4 




17:06:25.55 


+43:54:39.8 


3.1 ± 1.4 


2.35 


central 


B99 (SSA17.2) 


3.9 


0.6 




17:06:20.37 


+43:54:09.8 


2.6 ± 2.6 


1.02 


edge 


B99 (SSA17.5) 


3.1 


5.8 



Table 13. 850 ^lm source list for the SSA17 Field of the "Hawaii Flanking Fields Survey". Sources are marked in Fig 16. 



those potential sources common to both catalogues, the com- 
bination of the IDL-rcduction and simultaneous maximum- 
likelihood source extraction algorithm provide much more 
conservative values of the signal-to-noise ratio, in most cases 
by 1 — 2(7. At this stage it is not possible to say which of 
the two reduction algorithms is the more accurate. The deep 
radio imaging of the "SCUBA 8 mJy Survey" fields (Ivison 
et al. 2002) has provided quite a stringent test of the IDL- 
reduction and the simultaneous maximum-likelihood source 



extraction algorithm technique, robustly detecting 50% of 

the bright SCUBA sources uncovered in the "8mJy Sur- 
vey", and less significantly detecting a further 20% of the 
objects, suggesting a 30% upper limit on the contamina- 
tion from spurious / confused sources. A similar analysis of 
the CUDSS fields, however, is unlikely to be very informa- 
tive regarding a comparison of the two independent source 
lists, since the CUDSS fields are smaller and deeper than the 
"8 mJy Survey" fields, designed with the aim of studying less 
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RA 


DEC 


Ssso 


S/N 


Noise 


Previous 


Prev. 


Sep. 




(J 2000) 


(J 2000 J 


/mjy 




Region 


Reierence 


S/N / 


arcsec 




01 22:17:33.96 


+00:13:53.4 


4.7 + 0.8 


6.97 


central 


B99 (SSA22.1) 


6.9 


2.2 




02 22:17:35.03 


+00:15:36.4 


2.7 + 0.7 


4.29 


central 


B99 (SSA22.2) 


5.3 


5.1 




03 22:17:31.23 


+00:16:07.4 


3.1 + 0.8 


4.04 


central 










04 22:17:41.56 


+00:16:04.4 


3.7+1.0 


3.96 


central 


B99 (SSA22.5) 


3.1 


7.0 




05 22:17:40.90 


+00:14:56.4 


3.0 + 0.8 


3.92 


central 










06 22:17:35.96 


+00:15:56.4 


2.5 ± 0.7 


3.77 


central 


B99 (SSA22.3) 


4.0 


13.1 




07 22:17:37.36 


+00:16:21.4 


3.6+1.0 


3.75 


central 










08 22:17:29.30 


+00:13:57.4 


2.6 + 0.7 


3.66 


central 










09 22:17:19.96 


+00:15:25.4 


8.1 + 2.7 


3.22 


edge 










10 22:17:33.43 


+00:16:13.4 


2.4 + 0.8 


3.11 


central 










11 22:17:33.76 


+00:15:42.4 


1.8 + 0.6 


3.10 


central 


B99 (SSA22.4) 


o.o 


O.O 




12 22:17:41.03 


+00:13:32.4 


3.2 + 1.1 


3.08 


central 








S50 fim source list for the SSA22 Field of the "Hawaii Flanking Fields Survey" . Sources 


are marked 


in Fig 17. 




RA 


DEC 


Ssso 


S/N 


Noise 


Previous 


Prev. 


Sep. 




(J2000) 


(J2000) 


/mJy 




Region 


Reference 


S/N 


/ arcsec 


01 


12:36:52.01 


+62:12:27.0 


5.5 + 0.7 


12.13 


central 


H98 (HDF.l) 


14.0 


2.2 














Serj03 (HDF.l) 


15.3 


1.6 


02 


12:36:56.59 


+62:12:06.0 


3.6 + 0.6 


7.06 


central 


H98 (HDF.2) 


5.4 


2.3 














Serj03 (HDF.2) 


7.6 


2.6 


03 


12:36:53.16 


+62:13:55.0 


2.7 + 0.6 


4.79 


central 


Serj03 (HDF.8) 


3.5 


0.8 


04 


12:36:50.58 


+62:13:18.0 


2.1 + 0.5 


4.35 


central 


H98 (HDF.4)* 


4.6 


2.6 














H98 (HDF.5)* 


4.2 


o o 
y.y 














Serj03 (HDF.4&5) 5.1 


n/a 


05 


12:36:44.87 


+62:11:40.0 


2.9 + 0.8 


3.87 


central 








06 


12:37:03.17 


+62:13:04.0 


2.8 + 0.8 


3.66 


central 










12:36: U- 00 


+62:13:09.0 


1.2± 0.6 


2.13 


central 


H98 (HDF.3) 


5.0 


7.5 














Serj03 (HDF.3) 


2.1 


2.9 




n/a 


n/a 


n/ a 


n/a 


edge 


Serj03 (HDF.6) 


3.8 


n/a 




n/a 


n/a 


n/ a 


n/a 


edge 


Serj03 (HDF.7) 


3.7 


n/a 



Table 15. 850 /um source list for the "Hubble Deep Field Survey". Sources are marked in Fig 18. 



bright SCUBA sources in the range 3mJy < Ssso < 6mJy, 
and hence only the very brightest (and most significant of 
the CUDSS sources) arc likely to be detected in a radio im- 
age even at a depth of ai.4GHz = 5 /xJy/beam. 



5.3 The Hawaii Submillimetre Survey 

The "Hawaii Flanking Fields Survey (HFFS)" (Barger et 
al. 1998, Barger, Cowie & Sanders 1999, Barger, Cowie & 
Richards 2000), covers a total of ~ 110 sq. arcmin over 4 
regions of sky. The Lockman Hole deep field is a small pen- 
cil beam map (Barger et al. 1998), covering ~ 8 sq. arcmin 
of sky to an rms noise level of asso = 0.8mJy/beam. The 
SSA13 field is composed of a deep pencil beam area (~ 8 sq. 
arcmin in size, crsso = 0.7mJy/beam; Barger et al. 1998), 
embedded in a wider-area, shallower map covering an addi- 
tional 45 sq. arcmin with a typical rms noise level of ugso = 
2.5mJy/bcam (Barger, Cowie & Sanders 1999). The SSA17 
and SSA22 fields (Barger, Cowie & Sanders 1999) both have 
regions of uniform noise covering approximately 20 sq. ar- 
cmin of sky, the SSA17 field to usso = 1.6mJy/beam, and 
the SSA22 field to ctsso = 0.9 mjy/beam. 

The data were originally reduced using the standard 



SURF pipeline, and source extraction was carried out by 
convolving the signal maps with the beam. In order to de- 
termine the absolute noise levels, Barger, Cowie & Sanders 
(1999) first eliminated any significant sources (estimated 
to be at the > 2.8ct level) by subtracting appropriately 
normalised versions of the beam profile. They then placed 
beam-sized apertures at random positions on the residual 
signal map, using the standard deviation between the en- 
closed pixels to estimate the noise. These values were then 
used to iteratively adjust the normalization of the variance 
array values, until the dispersion of the signal-to-noise ratio 
was approximately one. 

Comparative source lists for detections with S/N > 3.00 

from the two different reductions and source extraction pro- 
cedures are given in Tables 11-14 (see also Figs 14 - 17). 
Columns 1-6, 8 and 9 are analogous to those of Tables 7 to 
10. 

In one case, a source previously identified by Barger, 
Cowie & Sanders (1999) has been marked with 'bad boL' 
as this region appears to have been observed with a bad 
bolometer, hence the "source" is most likely to be an arte- 
fact of surf's attempt to interpolate between neighbour- 
ing areas of good quality data. Column 7 gives any previous 
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Figure 14. The 850 image of the Loekman Hole field from 
the Hawaii Flanking Fields Survey, smoothed with a beam-size 
Gaussian (14.5 arcsec FWHM). The numbered circles highlight 
those sources found at a significance of > 3.00. The labelling 
corresponds to the numbers in Table 11. 



Figure 16. The 850 /im image of the SSA17 field from the 
Hawaii Flanking Fields Survey, smoothed with a beam-size Gaus- 
sian (14.5 arcsec FWHM). The numbered circles highlight those 
sources found at a significance of > 3.00. The labelling corre- 
sponds to the numbers in Table 13. 




Figure 15. The 850 image of the SSA13 field from the 
Hawaii Flanking Fields Survey, smoothed with a beam-size Gaus- 
sian (14.5 arcsec FWHM). The numbered circles highlight those 
sources found at a significance of > 3.00. The labelling corre- 
sponds to the numbers in Table 12. 
© 0000 RAS, MNRAS 000, 000-000 



Figure 17. The 850 jum image of the SSA22 field from the 
Hawaii Flanking Fields Survey, smoothed with a beam-size Gaus- 
sian (14.5 arcsec FWHM). The numbered circles highlight those 
sources found at a significance of > 3.00. The labelling corre- 
sponds to the numbers in Table 14. 
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reference to the 850 jim source. Reference B99 is an abbre- 
viation for Bargcr, Cowio & Sanders (1999). 

The two reduction and source-extraction algorithms 
again produce quite different results. There is no clear trend 
in the variations between the independent measurements 
of signal-to-noise for those sources detected significantly in 
both analyses; in some cases the SURF reduction yields a 
higher signal-to-noise estimate whereEis for other sources the 
situation is reversed. The IDL-based reduction and max- 
imum likelihood algorithm has identified nearly twice as 
many peaks at the > S.OOu level as the original catalogue 
of Barger, Cowie & Sanders (1999), including all of their 
original > 4.00cr objects. The comparative results for the 
Lockman Hole deep field given in Table 11 show that of the 
two sources published in Bargcr, Cowic & Sanders (1999), 
their most significant source is confirmed at > 6(7. The most 
discrepant of the fields is the SSA13 field (Table 12), which 
has particulaly uneven noise when compared to all of the 
other survey fields and it is most likely because of this that 
only 4/9 of the original detections could be identified at 
S/N > 3.00. In the deep part of the SSA13 field, the top two 
sources are confirmed at better than 5.00cr in this analysis. 
However, in the wider shallower area only one further source 
listed in Barger, Cowie & Sanders (1999) could be confirmed 
at S/N > 4.00. Our re-analysis of the SSA17 field (Table 
13) found only one of the original sources with S/N > 3.50 
and with only 3/5 of the objects listed in Barger, Cowie & 
Sanders at > S.OOct. In the SSA22 field, all 5 of the > 3.00a 
detected by Barger, Cowie & Sanders (1999) are recovered, 
with firm agreement on the top 2 sources at better than 
4.00ct in both catalogues. 




2 1 0-1-2 

ARC MINUTES 

CENTRE: R.A. 12 36 51.66 DEC +62 12 56.5 J2000 



Figure 18. The 850 fim image of the Hubble Deep Field North, 
smoothed with a beam-size Gaussian (14.5 arcsec FWHM). The 
numbered circles highlight those sources found at a significance 
of > 3.00. The labelling corresponds to the numbers in Table 15. 



5.4 The Hubble Deep Field Submillimetre Survey 

The Hubble Deep field is the deepest of the submillimetre 
surveys, covering approximately 6 sq. arcmin of sky down 
to the confusion level of usso = 0.5 mjy/beam. It differs 
slightly from the other SCUBA surveys in that two different 
chop throws of 30" and 45", both fixed in celestial coordi- 
nates, were applied to the observations, each for approxi- 
mately half of the integration time. Both SURF and IDL 
reductions have previously been carried out on this field 
(Hughes ct al. 1998 and Serjeant et al. 2003 respectively), 
but the maximum-likelihood simultaneous-fitting algorithm 
has not been used previously. 

Comparative source lists for detections with S/N > 3.00 
from the three different reductions and source extraction 
procedures are given in Table 15 (see also Fig 18). Column 
1 gives the source number corresponding to the labelling on 
Fig. 18, in order of decreasing signal-to-noise. Columns 2 
and 3 give the right ascension and declination of the source 
in J2000 coordinates. Column 4 gives the simultaneously fit- 
ted 850 /iin flux densities of the sources. The error includes 
a 10% calibration error combined in quadrature. Column 5 
gives the measured signal-to-noise ratio of the source from 
the simultaneously fitted model. Column 6 defines the noise 
region in which the source was found; 'central' corresponds 
to the parts of the map which have seen the full integration 
time (outside of the deep area), and 'edge' corresponds to 
the rather noisier regions near the perimeter which have not 
seen the full integration time. Column 7 gives any previous 
reference to the 850 jj,m source. Reference H98 is an abbrevi- 



ation for Hughes et al. (1998), and Serj03 is an abbreviation 
for Serjeant et al. (2003). The presence of a * indicates that 
Hughes ct al. (1998) deconvolved two sources here, whereas 
Serjeant et al. (2003) and my own reduction extracted only 
one. Column 8 gives the previously recorded signal-to-noise 
ratio where applicable, and column 9 gives the distance be- 
tween the listed and previously referenced positions. The 
table listings given in italics correspond to previously ref- 
erenced sources with S/N > 3.00, which did not meet this 
criteria in this analysis. 

The analysis of Serjeant et al. (2003) combined the 
deep jigglemap imaging of the Hubble Deep Field North 
with additional photometry which both increased the map 
size slightly and improved the sampling at some positions 
around the edges. Since photometry measurements require 
pre-selcction of a known object's position this would bias the 
number of sources recovered in this region of sky and so this 
data has not been included in the analysis presented here. 
Consequently the sources labelled as HDF.6 and HDF.7 pre- 
sented in Serjeant et al. (2003) which do not peak within the 
main jigglemap imaging have not been re-identified in this 
analysis although they have been recovered in the "HDF-N 
Supermap" of Borys et al. (2003). However, two new sources 
which have eluded previous analyses, have been recovered at 
the S/N > 3.50 level and all of the original detections, with 
the exception of HDF.3 in the Hughes et al. (1998) analysis, 
have been confirmed with S/N > 4.00. 
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Flux density Raw 850 fj,m source counts Corrected 850 /im source counts 
/mjy N{> 5)deg-2 N{> 5)deg-2 

550 
1120 
320 
620 
200 
370 
120 
200 
110 
190 
90 
150 
70 
120 
70 
110 
60 
100 
50 
80 
50 
60 
50 
60 
40 
60 
30 
40 
30 
30 
30 
30 
20 
20 
10 
10 
10 
10 
10 
10 
10 
10 
10 
10 



2.0 


2880 ± 310 


39201 


2.5 


1680 ± 180 


2140t 


3.0 


1080 ± 120 


1250l 


3.5 


830 ± 100 


680t 


4.0 


700 ± 


90 


620l 


4.5 


700 ± 


90 


490l 


5.0 


540 ± 


70 


380t 


5.5 


500 ± 


70 


330l 


6.0 


420 ± 


60 


310l 


6.5 


340 ± 


60 


230l 


7.0 


300 ± 


50 


180t 


7.5 


260 ± 


50 


180l 


8.0 


210 ± 


40 


150l 


8.5 


160 ± 


40 


lOOl 


9.0 


130 ± 


30 


70t 


9.5 


90 ± 


30 


60t 


10.0 


70 ± 


20 


40t 


10.5 


60 ± 


20 


20+ 


11.0 


30 ± 


20 


10+ 


11.5 


20 ± 


10 


10+ 


12.0 


20 ± 


10 


10+ 


12.5 


10 ± 


10 


10+ 



Table 16. The 850 //m source counts per square degree based on sources with S/N > 3.50 in both survey maps, and excluding those 
detected in the non-uniform noise regions. Column 1 gives the flux density and column 2 the cumulative raw counts per square degree 
with the Poisson error. Column 3 gives the cumulative corrected counts per square degree, the upper error corresponding to the Poisson 
error, and the lower error accounting for both the Poisson error and the presence of spurious sources based on the simulation data. 



6 SOURCE COUNTS 

Previous measurements of the cumulative number counts 
have differed quite markedly between the various surveys 
(eg. Blain at al. 1999, Barger, Cowie & Sanders 1999, Eales 
et al. 2000, Scott et al. 2002), particularly at bright flux 
densities (> 5mJy). This is not surprising given the small 
area of sky observed by each individual survey and the low 
number density of bright sources. In addition to cosmic vari- 
ance, if the SCUBA sources have a tendency to cluster on 
scales of a few arcminutes (for which evidence is given in 
the next section) the problem is further compounded since 
it is more likely that a mapping a blank field of only several 
tens of square arcminutes in size will fall in an area of lower 
source density. Here, the sources detected with a signal-to- 
noise ratio > 3.50 from our rcanalysis of the CUDSS, Hawaii 
and HDF surveys have been combined with those identified 
in the 8mJy Survey, to produce the most accurate num- 
ber counts to date, from 2 — 12.5 mJy at 0.5 mJy intervals. 
The regions of non-uniform noise towards the edge of the 
maps and any sources they contained were excluded. The 
simulations described in Section 4.1 were used to correct 
for the effects of flux-density boosting and incompleteness 
on a field by field basis. Estimating the level of contami- 
nation from spurious / confused sources, however, requires 
the generation of fully simulated images. The accuracy of 
such images is hampered by the lack of knowledge regarding 
the clustering properties of the SCUBA population down to 
the faintest flux density levels. The simulations described in 
Section 4.2 imply that the fraction of spurious / confused 
sources (defined as having no input source at least half as 



bright as the output flux density) at a significance of > 3.50a" 
is ~ 30%. This number is in line with the upper limit placed 
on the confused / spurious fraction from deep radio follow- 
up of the "8mJy Survey" fields (Ivison et al. 2002). The raw 
and corrected cumulative number counts are given in Table 
16 along with la Poisson error bars (calculated from the 
square root of the number of detections on which the source 
count was based). The estimated 30% fraction of spurious / 
confused sources has been combined in quadrature with the 
lower Poisson error bar in the corrected number counts. 

Figure 19 shows the corrected versus the raw number 
counts at each fiux density, together with a dotted line mark- 
ing the locus of where the raw and corrected counts take the 
same value. An increase in the cumulative number counts 
along either axis correponds to a decrease in the flux den- 
sity threshold. One can see immediately that at brighter flux 
densities, the effect of boosting is stronger than incomplete- 
ness and hence the real number density of sources above 
a specified flux density threshold is lower than the directly 
measured value. This may also be implied from an exam- 
ination of the raw data alone. The area of sky mapped in 
the "SCUBA 8 mJy Survey" is almost identical to the area 
of sky mapped by the other deeper surveys combined. At 
bright fiux densities the deeper surveys are essentially com- 
plete and so one can compare the number of detections above 
a specified signal-to-noise ratio and flux density threshold 
between the shallower and deeper maps. For a retrieved 
SssoMm > lOmJy and S/N > 3.50, the "SCUBA 8mJy Sur- 
vey" identified 8 sources, whereas in all the other surveys 
only 1 source was found to satisfy this criteria. Similarly, for 
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Figure 19. Cumulative number counts corrected for the eflfects of boosting and incompleteness versus the raw number counts. The 
dotted line marks the case where the raw and corrected cumulative numbers counts are the same. The error bars are as given in Table 
16. 



a retrieved Sgso^m > 8mJy and S/N > 3.50, the "SCUBA 
8 mjy Survey" identified 20 objects, whereas only a total of 
5 were found in the other surveys. Given the small area of 
sky considered, part of the discrepancy between the number 
of bright sources recovered in the "SCUBA 8 mJy Survey" 
and the other deeper blank field surveys can be explained 
by small number statistics but it seems highly unlikely that 
this is the sole cause. Moving towards fainter flux density 
thresholds, the increased source density makes incomplete- 
ness more of a problem and around the confusion limit of 
~ 3 m,Jy (where the density of sources is ~ 1000 per square 
degree) the raw and corrected cumulative number counts are 
approximately the same. 

One of the simplest ways of describing the number 
counts is by carrying out a simple parametric fit of a power- 
law model to the differential number counts as a function 
of flux density. The differential number counts at a specific 
flux density were determined by the difference in cumula- 
tive counts between the two values on either side of the data 
point, divided by the change in flux density, which approx- 
imates to a measure of the gradient of a tangent to the cu- 
mulative number counts curve at that point. The complete- 
ness and boosting corrected number counts were used for 



this procedure. Following Barger, Cowie & Sanders (1999), 
a model of the form 



n(S): 



No 



dN(S) 

dS ~(a-hS«) 



(13) 



was fitted to the differential data points by means of a min- 
imised method, where the values of No, a and a were al- 
lowed to vary freely. The best fit values for these parameters 
were No = 2.67 x 10**, a = 0.49 and a = 3.14, predicting a to- 
tal 850 nm background of 3.8 x 10'' mJy deg~^, mid-way be- 
tween the 850 /um extragalactic values of 3.1 x 10* mJy deg~^ 
and 4.4 x lO" mJy deg-^ as measured from COBE-FIRAS 
by Puget et al. (1996) and Fixsen et al. (1998) for the lower 
and upper values respectively. The values of the fitted pa- 
rameters are in fact very close to the values originally de- 
termined by Barger, Cowie & Sanders (1999) who found 
No = 3.0 X 10", a = 0.5 and a = 3.2. Figure 20 shows a plot 
of the differential number counts against fiux density. The 
solid diamonds represent data from the combined blank field 
survey re-analysis only, with la Poisson error bars on the 
y-axis and a flux density range of 1 mJy on the x-axis, corre- 
sponding to the change in flux density between the two data 
points on either side of the point at which the differential 
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Figure 20. A plot of differential number counts vs. flux density. The solid diamonds represent data from the combined blank field survey 
re-analysis only, with la Poisson error bars on the y-axis, and the flux density range on which the differential count is based marked 
as an error bar on the x-axis. The open triangles represent data points from the lensing surveys of Blain et al. (1999) and Cowie et al. 
(2002). The solid curve is a best fit parametric model of the form n(s) = ^^^g„ ^ . 



counts has been determined. The open triangles represent 
data points from the lensing surveys of Blain et al. 1999 
and Cowie et al. 2002, and the solid curve is the best fit 
parametric model. 

Fig. 21 shows the corrected cumulative number counts, 
as well as a scries of models predicting the 850 fiin source 
counts for an adopted = 0.3, = 0.7 cosmology, in line 



with the most recent cosmological parameter results from 
the supernova cosmology project (Perlmutter et al. 1999), 
X-ray clusters (Allen, Schmidt & Fabian 2002a, 2002b), and 
WMAP observations of the cosmic microwave background 
(Spcrgcl et al. 2003). A value of Ho = 67kms"^Mpc"^ 
was assumed throughout. The black lines are parameterised 
models, the solid lines following a simple power-law descrip- 
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Flux/mJy 

Figure 21. A plot of cumulative number counts vs. flux density, along with a series of models assuming an Qm = 0.3, Qa = 0.7 
cosmology. The solid diamonds represent the cumulative source counts from the completeness and boosting corrected counts derived 
from the combined blank field survey re-analysis only, with error bars as given in Table 16. The crosses and round-edged squares represent 
data points from the lensing surveys of Blain et al. (1999) and Cowie et al. (2002) respectively, while the stars indicate the results of and 
Borys et al. (2003). The various models are described in the main text. 



tion of the data points and are the same for both figures. The 
dashed hnes represent work by Rowan-Robinson (2000,2001) 
who generated best-fit models of the measured far-infrared 
through to submillimetre number counts, as constrained by 
the extragalactic background, for a number of assumed cos- 
mologies. In the remaining models (red, blue and green 
curves), values of dust temperature Td and emissivity index 



13 based on optically thin greybody emission, were allowed 
to vary where necessary so that as good a fit as possible 
could be obtained for the = 0.3, Q,a = 0.7 cosmology. 

The solid black line is the integrated dN/dS best-fit 
parametric model given in equation 13, providing a good 
description of the cumulative number counts for Ssso < 
lOmJy. Rowan-Robinson (2001) has also taken a param- 
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eterised approach to modelling the infrared through to sub- 
millimctrc counts and backgrounds, in this case by us- 
ing multiwavelength observational data in those wavelength 
regimes to place constraints on models made up of four spec- 
tral components: infrared cirrus, an M82-like starburst, an 
Arp220-like starburst and an AGN torus. The model as- 
sumes that the evolution of the star formation rate manifests 
itself as pure luminosity evolution. The models, consistent 
with infrared and submillimetre counts and backgrounds, 
showed a flat star formation rate from « = 1 — 3, in agree- 
ment with other studies of the star formation history such as 
the HDF (Hughes et al. 1998). The most striking difference 
between other modelling work and Rowan-Robinson (2001) 
is the dominant role of the cirrus component at submillime- 
tre wavelengths. The black dashed line in Fig. 21 shows the 
model of Rowan-Robinson (2001) which appears to be the 
best-fit of those presented here to the cumulative 850 /im 
source count data points over the whole flux density range. 

The red lines on Fig. 21 represent a pure luminosity 
evolution of the form (1 + z)"^ out to a threshold rcdshift, 
and constant thereafter. The solid line began with the 60 /im 
luminosity function of Saunders et al. (1990), interpolated 
to 850 /um assuming an optically thin greybody with a single 
dust temperature and emissivity index, to describe the far- 
infrared / submm dust SED. The same values for Td and ,5 
were assumed in calculating the K-correction for consistency. 
The dashed line began with the 850 fim luminosity function 
measured by Dunne et al. (2000). A dust temperature of 
40 K and emissivity index of /3 = 1.3 appeared to best satisfy 
the data points, using threshold redshifts of z = 2.0 and 
z = 1.5 for the Saunders et al. (1990) and Dunne et al. 
(2000) based luminosity functions respectively. This simple 
description of luminosity evolution appeared to work well 
for the interpolated Saunders et al. (1990) 60 pm luminosity 
function for Sgso < 8 mjy, but predicted too many sources 
brighter than this. Conversely, using the Dunne et al. (2000) 
850 /xm luminosity function, a good fit to the data points 
was found for Ssso > 2 mJy but overpredicted the number 
of sources fainter than this by a factor of 2-3. 

The solid and dashed blue lines on the cumulative num- 
ber density plot also assume pure luminosity evolution of 
the interpolated 60 /xm Saunders et al. (1990) and directly 
measured 850 /im Dunne et al. (2000) luminosity functions 
respectively. In this case, however, the luminosity evolution 
g{z) takes a more realistic form, which is fully compatible 
with models of cosmic chemical evolution and naturally in- 
cludes a peak in the evolution function (Jameson et al. 1999, 
Small et al. 2002): 



g{z) = {1 + «)^''^sech^[bln(l + z) — c]cosh^c 



(14) 



where best fits values for the parameters b and c based on 
multiwavelength far-infrared to submillimetre counts and 
constraints on the probable redshift distribution of the 
850 /um population (Small et al. 2002) axe b = 2.2 ± 0.1 and 
c = 1.84 ± 0.1. The interpolated 60 nm luminosity function 
combined with this description of pure luminosity evolution 
fits the data well for Ssso < 8 mJy assuming an optically thin 
greybody description of the thermal dust emission with pa- 
rameters Td = 39K and l3 = 1.3. The data points across the 
whole flux density range are fltted well using the Dunne et 
al. (2000) 850 fiin luminosity function assuming dust emis- 
sion parameters of Td = 30 K and /3 = 1.0, although it 



should be noted that this combination of dust temperature 
and emissivity index does not describe the spectral energy 
distribution well for any known local galaxies based on the 
"SCUBA Local Universe Galaxy Survey (SLUGS)" (Dunne 
et al. 2000). 

The solid green line on Fig. 21 uses the 60 /xm luminosity 
function of Saunders et al. (1990) combined with luminosity 
evolution taking the form (1 -|- up to a threshold redshift 
^thresh, and (1 -|- z)~^ thereafter (Chapman et al. 2002a). 
They find a best fit of athresh = 2.6 for the transitional red- 
shift, based on current constraints on the redshift distribu- 
tion as well as the source counts. This particular 850 /im 
counts model has a rather shallower gradient than the oth- 
ers and could only be made to match number counts of this 
combined re-analysis over the fiux density range ^ 3 — 7mJy 
for Td = 36 K and a very extreme /3 = 2.0. Possibly a 
higher value of the dust emissivity index could produce a 
steeper counts model, however, such high values are not ob- 
served in in any known class of objects. No satisfactory fit 
could be obtained beginning with the Dunne et al. (2000) 
850 /im luminosity function. The proposed Chapman et al. 
(2002a) luminosity evolution does not seem consistent with 
the blank field survey source counts, although the number 
counts derived from the lensing survey of Blain et al. 1999 
are much shallower, predicting rather higher number den- 
sities of bright sources, and could plausibly be fit by this 
scenario for less extreme values of beta. 

One particularly interesting feature of this source 
counts re-analysis is the apparent steepening of the cumula- 
tive number counts beyond Sgso > 8 mJy. This could at least 
in part be due to small number statistics of the brightest 
sources, however, if real has some important implications. 
Firstly, it could indicate an intrinsic turn-over in the un- 
derlying luminosity function. In turn, this would suggest a 
very interesting upper limit on the luminosity of a high red- 
shift galaxy, perhaps refiecting an upper limit on the over- 
all mass of the system. This could place useful constraints 
on galaxy formation theories. Secondly, a steepening of the 
bright source counts could make the SCUBA population 
much more prone to the effects of gravitational lensing. For 
weaJi lensing scenarios, the ratio of the number of observed 
sources brighter than a flux density threshold to the true 
number of sources brighter than that flux density threshold 
is given by 



Nobs(> S) 

Ntrue(> S) 



(15) 



where /u is the magnification amplitude and 7 is the slope 
of the cumulative source counts. For 1 < Ssso < 8mJy the 
counts slope 7 ~ 2.5, whereas for Ssso > 8mJy this in- 
creases to 7 ~ 5.5. No attempt has been made to correct the 
850 /xm source counts for the effects of lensing, but there is 
some evidence to suggest that some blank field sources have 
been gravitationally lensed. Almaini et al. (2003) found a 
strong cross-correlation in the ELAIS N2 field between the 
distribution of SCUBA and Chandra sources even though 
the coincidence of detections in the submillimetre and X-ray 
wavelengths was only ~ 5%. One proposed explanation for 
this effect is that the SCUBA and Chandra sources trace the 
same large scale structure at high redshift {z > 1), however, 
there also appears to be a similarly strong cross- correlation 
between the low redshift 7— band sources and the higher red- 
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shift SCUBA and Chandra sources in this field (Almaini et 
al. 2005) suggesting an alternative explanation may be valid: 
the SCUBA and Chandra detections may have been mag- 
nified in certain regions of the field by the presence of high 
mass density structure, as traced by the /—band imaging, 
at 2 ~ 0.5. Chapman et al. (2002b) have also pointed out 
that some submillimetre sources have apparent counterparts 
which are optically bright galaxies, I < 21.5, lying at modest 
redshifts, z < 1. This could of course be explained by these 
counterparts being correct, reflecting a population of galax- 
ies which are detected as very cold, luminous submillimetre 
sources, however, a second explanation for such systems is 
that the optically bright galaxy is a foreground object acting 
as a gravitational lens, amplifying the more distant SCUBA 
galaxy. The detection of luminous molecular CO emission 
at the rcdsliift of the optically bright galaxy would provide 
a powerful test to distinguish between these two scenarios. 
If the latter explanation of gravitational lensing is correct. 
Chapman et al. (2002b) estimate that up to 3 — 5% of the 
> 10 mjy submillimetre sources detected in blank field sur- 
veys could be gravitationally amplified by foreground galax- 
ies. 

Gravitational lensing by clusters of galaxies can of 
course be used to study the fainter submillimetre sources, 
and this technique has been successfully applied by Blain et 
al. (1999) and Cowie et al. (2002) (the data points marked 
by the crosses and curved-edged squares on Fig. 21 respec- 
tively). The agreement between the source counts derived 
from the two cluster-lensing surveys is extremely good at 
fiux densities below Sgso ^ 2 mJy, and there is a smooth 
transition at this point between the faint number counts 
from the cluster-lensing surveys and the brighter source 
counts from our combined reanalysis of the blank field sur- 
veys. However, the number counts derived by Blain et al. 
(1999) at 4 and 8mJy are significantly higher than the com- 
bined blank field survey counts presented in Table 16. The 
most likely explanation for this discrepancy is that the num- 
bers quoted by Blain et al. (1999) suffer badly from small 
number statistics at bright flux densities, their lensing sur- 
vey being composed of a number of small fields covering a 
total area of sky of only a few tens of square arcminutes. 
In comparison, the total area of sky observed by the "8 mJy 
Survey", "CUDSS", "HDF" and "Hawau Survey" is 460 sq. 
arcmin, an order of magnitude larger. If the bright SCUBA 
population does strongly cluster on arcminute scales, as im- 
plied by the evidence presented in Section 7, this could also 
affect the bright end of the number counts in small area 
surveys such as Blain et al. (1999). 

The data points marked on Fig. 21 by open stars are 
determined from the 165 sq. arcmin "Hubble Deep Field 
North SCUBA Super-map" presented by Borys et al. (2003). 
The rms noise level of of this composite map ranges from 
cTrms — 0.5 — 7mJy / beam (Borys et al. 2003). In general 
the counts of Borys et al. (2003) are slightly higher than 
those determined from this re-analysis, although they agree 
to well within the error bars for Ssso < lOmJy. At brighter 
fiux densities, the two sets of cumulative number counts be- 
gin to diverge, those of Borys et al. (2003) continuing to fol- 
low an approximately power-law decline with increasing flux 
density whereas those presented in this paper steepen quite 
markedly. There are several possible reasons for this discrep- 
ancy. Firstly both of these analyses suffer from small number 



statistics at the brightest flux densities (Ssso > 10 mJy) and 
so cosmic variance may be an issue particularly if these ob- 
jects tend to cluster on arcminute scales. Secondly Borys et 
al. had problems calibrating their wide area scan-map data 
which represents the shallowest regions of the super-map. 
If they have over-estimated their flux conversion factors by 
~ 20% the two sets of number count measurements agree 
well. Finally 4/5 of the sources recovered by Borys et al. 
at S/N > 3.5 and Ssso > 16mJy have rms noise levels of 
Crms > 5mJy/beam suggesting that these may have been 
recovered in non-uniform edge regions for which the simula- 
tions presented in section 4.2 have shown to have very poor 
source reliability. This would obviously increase the cumula- 
tive number count measurements at all flux densities fainter 
than this. 

Overall, the increased accuracy of this new source 
counts analysis will allow the evolutionary nature, as well 
as the dust and star-forming properties of the submillime- 
tre population to be studied in much more detail than has 
been possible before. This will be particularly useful when 
combined with knowledge of the redshift distribution of sub- 
millimetre sources, which is slowly being built up at the 
present. 



7 CLUSTERING 

If the bright 850 fim sources are indeed the progenitors 
of massive elliptical galaxies then they should be strongly 
clustered, an inevitable result of gravitational collapse from 
Gaussian initial density fluctuations since the rare high-mass 
peaks are strongly biased with respect to the mass. There is 
a great deal of evidence to support the presence of this bias 
at high redshift. The correlations of Lyman-break galaxies 
at 2 ~ 3 (Steidel et al. 1999) are almost identical to those 
of present-day field galaxies, even though the mass must 
have been much more uniform at early times. Furthermore, 
the correlations increase with UV luminosity (Giavalisco & 
Dickinson 2001) reaching scale lengths of ro ^ 7.5h~^ Mpc - 
approximately 1.5 times the present-day value. In the case of 
luminous proto-ellipticals an even stronger bias is expected 
since one is selecting not just massive galaxies but those that 
have collapsed particularly early in order to generate the old- 
est stellar populations. This is suggested by studies of the 
local Universe which have shown that early-type galaxies are 
much more clustered than late- type galaxies (eg. Guzzo et al. 
1997, Willmer et al. 1998), and more recently by the findings 
of Daddi et al. (2000) who have investigated the clustering 
properties of extremely red objects (EROs). They detect a 
strong clustering signal of the EROs which is about an or- 
der of magnitude larger than the clustering of if— selected 
field galaxies, and also report a smooth trend of increasing 
clustering amplitude with increasing R — K colour, reaching 
ro ^ llft~^ Mpc for i? — if > 5. These results axe probar- 
bly the strongest evidence to date that the largest fraction 
of EROs is composed of ellipticals at 2 > 1. There are al- 
ready some hints of strong clustering in the bright sub-mm 
population from the discovery of a strong excess of bright 
SCUBA sources axound high-redshift AGN (Ivison et al. 
2000, Stevens et al. 2003). Blain et al. (2004) have also com- 
bined a sample of 73 submillimetre galaxies over seven fields 
for which they had spectroscopic redshifts, finding tentative 
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evidence for strong clustering of submillimetre galaxies with 
a correlation length ro ~ (6.9 ± 2.1)h~^ Mpc. 



7.1 Angular correlation functions 

The angular correlation function w(^) is the projection of 
the spatial function on the sky and is defined in terms of 
the joint probability 5P of finding two galaxies separated by 
an angular distance 6 with respect to that expected for a 
random distribution 



6P = N^[i + w{e)]snidn2 



(16) 



where 5Qi and 5ii2 are elements of solid angle, and N is the 
mean surface density of objects. If w(6) = the distribution 
is homgeneous. A positive w(0), therefore, corresponds to an 
over-density of sources separated by distance 9. 

There are a variety of possible estimators for w{9) as 
a function of pair-count ratios. Following Landy & Szalay 
(1993), we have adopted the estimator: 



w{e) 



(DD - 2DR + RR) 
RR 



(17) 



in which the variance is minimised to almost Poisson level. 
DD is the number of distinct data pairs in the real im- 
age within a bin covering a specified range of 0, DR is the 
number of cross-pairs between the real and mock catalogues 
within the same range of 0, and RR is the number of random- 
random pairs. DR and RR are normalised with respect to 
the total number of data-data pairs from the real image. 

Scott et al. (2002) made the first attempt at measuring 
an angular correlation function for the bright SCUBA pop- 
ulation. Their results provided tentative evidence for strong 
clustering of the submillimetre sources on scales of 1-2 ar- 
cmin but proved inconclusive due to the small number of 
sources recovered in the "8 mjy Survey" alone. More recent 
attempts to measure the clustering strength have relied on 
combining catalogues of submillimetre sources with Lyman- 
break galaxies (Webb et al. 2003b) or Chandra sources (Al- 
maini et al. 2003) and assuming that even though the ob- 
jects were rarely coincident, the two populations identified 
at difi^erent wavelengths were tracing the same large scale 
structure. These results further supported the view that the 
bright SCUBA galaxies cluster strongly on arcminute scales, 
however, in both cases the combined correlation functions 
were dominated by the larger number of Lyman-break galax- 
ies or Chandra sources. 

Further SCUBA galaxy clustering measurements, some 
tentatively detecting clustering, have been made (Blain et 
al. 2004, Borys et al. 2004, Webb et al. 2003a). 

Here, the 2-point angular correlation function analysis 
combines the sources identified in the "8 mJy Survey" fields 
with those objects detected in this re-reduction of the other 
blank field surveys, effectively creating a master catalogue 
with sources from double the area of the "8 mJy Survey" 
alone. 

A catalogue of randomly placed fake sources was cre- 
ated for each of the survey fields, with the number of fake 
objects contained in each field's random catalogue chosen 
to be directly proportional to the area of the image (rmm- 
ber in mock catalogue = 0.01 x Arca(/arcsec^)), so as to 
refiect a uniform number density across the sky. The wider 



area fields which dominate this analysis are largely signal-to- 
noise limited even at bright flux densities, and hence there is 
an increasing probability of finding significant detections in 
deeper regions of the maps. Although the positions were al- 
located randomly, the Gaussian convolved noise maps were 
used to weight the number density of sources across the 
image, since a larger density of sources above a specified 
signal-to-noise threshold would be expected in regions of 
lower noise. In practice, this meant dividing the full image 
into a series of sub-images, 20 arcsec by 20 arcsec in size, 
and calculating the mean noise level in each of these grid 
sections. For sources brighter than ~ 5 mJy, the rmmber of 
sources expected above a constant signal-to-noise threshold 
increases approximately as N(> Sssojim) oc Sggg^^^, and con- 
sequently the relative number of sources in each sub-image 
scales as (noise)"^'^. The fake source positions were then 
allocated according to a Poisson distribution, masking any 
positions which were covered by a negative sidelobe accom- 
panying a significant source in the real survey data. If the 
number of fake sources included in the catalogue is very 
much larger than the number recovered from the real data 
(as is the case here), this is in essence equivalent to choosing 
positions by combining the results from source extraction on 
a series of fully simulated survey fields. Detections lying in 
non-uniform regions of the maps in both the real and fake 
catalogues were omitted from the correlation function calcu- 
lations since many of the objects recovered in the real data 
in these high noise regions were likely to be spurious (see 
Section 4.2). 

Although combining all existing blank field survey data 
together does increase the number of sources on which a cor- 
relation function can be based, the numbers are still fairly 
small. Above a signal-to-noise ratio of 3.50, there are a to- 
tal of 53 sources detected brighter than 5mJy, of which 51 
may be used in a correlation function analysis (the other 2 
are lone sources in two of the deep single SCUBA pointing 
fields). This number increases to 104 sources brighter than 
5mJy at > 3.00a, of which 100 may be used in a correlation 
function analysis. However, this will also increase the con- 
tribution from spurious sources in the catalogue and may 
therefore dilute the measured signal. Figs. 22 and 23 show 
the directly measured angular correlation data points for 
detections brighter than 5 mJy, with la Poisson error bars, 
and for signal-to-noise thresholds of > 3.50 and > 3.00 re- 
spectively. The bin-size used in both plots is 29 arcsec (twice 
the beam-size). Rather surprisingly, there is little change in 
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Figure 22. 2-point angular correlation function for sources 
brighter than 5mjy, detected at a significance of > 3.50a, over 
all of the survey fields. The error bars are 1<t Poisson errors. The 
solid power-law line indicates the correlation function found by 
Daddi et al. (2000) for EROs with R - K > 5 and K < 18.5, 
and the dashed power-law line indicates the correlation function 
found by Giavalisco et al. (1998) for Lyman break galaxies at 
z ~ 3. The dot-dash line shows the best-fit power-law to the data 
points. The vertical dotted line indicates the size of the JCMT 
beam at 850 /im. 



Figure 24. 2-point angular correlation function for sources 
brighter than 6mjy, detected at a significance of > S.OOir, over 
all of the survey fields. The error bars are la Poisson errors. The 
solid power-law line indicates the correlation function found by 
Daddi et al. (2000) for EROs with R - K > 5 and K < 18.5, 
and the dashed power-law line indicates the correlation function 
found by Giavalisco et al. (1998) for Lyman break galaxies at 
2 ~ 3. The dot-dash line shows the best-fit power-law to the data 
points. The vertical dotted line indicates the size of the JCMT 
beam at 850 /im. 
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Figure 23. 2-point angular correlation function for sources 
brighter than 5mjy, detected at a significance of > S.OOo", over 
all of the survey fields. The error bars are Ict Poisson errors. The 
solid power-law line indicates the correlation function found by 
Daddi et al. (2000) for EROs with R - K > 5 and K < 18.5, 
and the dashed power-law line indicates the correlation function 
found by Giavalisco et al. (1998) for Lyman breaJj galaxies at 
z ~ 3. The dot-dash line shows the best-fit power-law to the data 
points. The vertical dotted line indicates the size of the JCMT 
beam at 850 /im. 



Figure 25. 2-point angular correlation function for sources 
brighter than 7mjy, detected at a significance of > 3.00a-, over 
all of the survey fields. The error bars are la Poisson errors. The 
solid power-law line indicates the correlation function found by 
Daddi et al. (2000) for EROs with R - K > 5 and K < 18.5, 
and the dashed power-law line indicates the correlation function 
found by Giavalisco et al. (1998) for Lyman break galaxies at 
z ~ 3. The dot-dash line shows the best-fit power-law to the data 
points. The vertical dotted line indicates the size of the JCMT 
beam at 850 /tm. 
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the proportion of excess pairs between corresponding data 
points in the first 100 arcscc of the two plots. The Gaussian- 
convolved noise images did not imply that there was any- 
thing unusual about the survey data in regions of higher 
than average 3.00cr source density, indicating that elevated 
noise where there is a concentration of significant detections 
is unlikely to be the cause. There are two likely explanations 
as to why the signal from the > 3.50cr catalogue is not di- 
luted in the > S.OOcr plot. Firstly, the simulation results in 
Section 4.2, which made estimates of the fraction of spuri- 
ous / confused sources at various signal-to-noise thresholds, 
have been pessimistic. It is diflicult to generate accurate 
realisations of the SCUBA maps on which to make esti- 
mates of contamination in a source catalogue due to con- 
fusion and noise without knowing the clustering properties 
of the SCUBA population in the first place, and the fact 
that there is very little difi'erence in the proportion of excess 
pairs on scales of ~ 1 arcmin between the two signal-to-noise 
datasets may be suggesting that most of the 3.00 — 3.49(7 
sources are real. Perhaps a more likely scenario, however, is 
that the clustering properties of the SCUBA population are 
strong over a significant range of fiux densities, and thus de- 
tections arising from the confusion of faint objects are more 
likely to be found on the same scale and in the vicinity of the 
real bright SCUBA sources. The major difference between 
the > 3.00 and > 3.50a plots is the size of the Poisson error 
bars. As one might expect, the error bars are a factor of 
~ \/2 smaller in Fig. 23, reflecting the increase by a factor 
of 2 of the number of sources considered, increasing the sig- 
nificance of the number of excess pairs contained in the first 
3 bins (14.5 - 101.5 arcsec) from 2.5 to 3.7a. 

Cutting at higher flux densities decreases the number 
of sources available for producing a 2-point angular corre- 
lation function even further, thus using the smaller "safer" 
catalogue of sources above a sigual-to-noise threshold of 3.50 
leads to tentative but non-significant measures of the clus- 
tering strength above the noise level. Instead, the analysis 
has been repeated using the full > 3.00a source lists for flux 
density cut-offs of > 6 and > 7mJy shown in Figs. 24 and 
25 respectively, pushing the available data to its limits. It 
should be noted that these flux density thresholds are the 
raw values as measured directly from the maps and should 
therefore not be taken as "absolute" because of the boost- 
ing effects described in Section 4.1. However, it is fair to 
say that cutting at Ssso > 7mJy defines a set of objects 
with generally higher star formation rates than cutting at 
Sgso > 5 mjy and so Figs. 23 to 25 can still be used to look 
for any trends in clustering strength with increasing flux 
density. For comparison, the number of excess pairs within 
~ 100 arcsec of each other in the real data when compared 
to a random distribution is significant at the 3.3a level for 
Ssso/jm > 6mJy, and at the 3.5a level for Sgso/jm > 7mJy. 

For all three fiux density limits, the measured corre- 
lation functions indicate a clustering strength much larger 
than that measured by Giavalisco et al. (1998) for Lyman- 
break galaxies at 2: ~ 3 (the dashed lines in the figures), 
even given the rather large error bars on the SCUBA data 
points. Could the apparent difference in clustering strength 
simply be due to projection effects over redshift space? 
The Lyman-break technique uses colour selection to identify 
high-redshift galaxies through multi-band imaging across 
the 1216A fine and the 912A Lyman break. At 2; > 2.5 the 



Lyman limit is redshifted far enough into the optical win- 
dow to be observable in broad-band ground-based photome- 
try. By placing filters on either side of the redshifted Lyman 
limit one can find high-redshift objects by their strong spec- 
tral breaks. Giavalisco et al. (1998) used a custom photo- 
metric system, U„GR (Steidel & Hamilton 1993) optimized 
for selecting Lyman-break galaxies with 2; ~ 3. By the na- 
ture of this method, 90% of the galaxies they used were 
confined to the redshift range 2.6 < z < 3.4, with none at 
z < 2.2. The redshift range of the bright SCUBA popula- 
tion used in these calculations of angular correlation func- 
tions is much more uncertain, but Ivisou et al. (2002) have 
suggested a median redshift of « = 2.4 based on the radio-to- 
submillimetre spectral indices (Carilli & Yun, 1999 & 2000), 
with inferred redshifts spanning the range of 2 ~ 1 — 4. 
Therefore, unless the bright SCUBA population occupies a 
very much narrower redshift band than implied by the radio- 
to-submillimetre spectral indices the dilution of the angular 
clustering signal by projection over redshift cannot be the 
reason for the large difference in clustering strength. 

Instead, this contrast in clustering properties implies 
that the Lyman-break galaxies and bright SCUBA sources 
are sampling two different stages or mass domains in galaxy 
formation (see also Barger, Cowie & Richards 2000, Webb et 
al. 2003b). The stronger clustering exhibited by the bright 
850 yum sources suggests that these objects are tracing the 
rarest high- mass peaks of the Gaussian initial density fluctu- 
ations and are the progenitors of the most massive ellipticals, 
whereas the weaker clustering of the Lyman-break galaxies 
indicates a weaker bias with respect to mass, detecting the 
formation of smaller disk or bulge systems. The strength 
of clustering, however, is consistent with that measured by 
Daddi et al. (2000) for extremely red objects (EROs) with 
K < 18.5 and R — K > 5 (the solid lines in the flgures), 
perhaps suggesting an evolutionary sequence from SCUBA 
source to ERG. The numerical coincidence between the co- 
moving number densities of EROs and that estimated for 
bright SCUBA sources is also in line with this idea (Scott 
et al. 2002). 

Attempts have been made to fit a standard power- 
law describing the angular correlation function to the data 
points shown in Figs. 22 to 25. The fit is slightly complicated 

by the fact that the global number density is unknown and 
must be estimated from the sample to hand. Hence, 

(l-hWt.ue(6')) = (l+WobB(e)) (^^)' (18) 

where w(&)truo is the true correlation function, w(6')obs is 
the observed correlation function, ii is the mean rmmber 
density of sources and < n > is the global number density of 
sources. One can rewrite the observed mean number density 
n in terms of a pertubation 6n on the global number density 
< n > (ie. ii n > ±<5n), such that in averaging over a 
set of sky areas the size of the sample 



(l + wt,uc(e)) = (l+Wobs(e))(l + a2) 



(19) 



where a^ is the rms surface density variation. If the real 
correlation functions take the power- law form Wtrue(^) = 



A^-^ = (S)-"'^then 



WobB(6') 



{Ae- 



1 + 



(20) 
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Flux density 
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A(0/deg) 
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<t2 
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1 circsec 


> 5 mjy 


> 3.50 


5.00 > 


< 10^2 


0.70 


0.386 


49.9 


> 5 mJy 


> 3.00 


1.82 > 


< 10^2 


0.89 


0.230 


39.9 


> 6 mJy 


> 3.00 


2.49 > 


< 10^3 


1.36 


0.280 


43.8 


> 7mJy 


> 3.00 


1.31 X 10-3 


1.56 


0.265 


46.7 



Table 17. Best-fit parameters to the angular correlation functions shown in Figs. 22 to 25. Column 1 gives the flux density cutoff as 
measured from the raw images and Column 2 gives the signal-to-noise cutoff for source detection. Colum 3 gives the amplitude A, Column 
4 gives the slope of the power-law 5, Column 5 gives the integral constraint <t^ and Column 6 gives the characteristic angular separation 
So. 



All SCUBA Surveys, S>5mJy, S/N>3.0 



50 100 150 200 

B /arcseconds 



250 



300 



30 



All SCUBA Surveys, S>6mJy, S/N>3.0 



50 100 150 200 

B /arcseconds 



250 300 



Figure 26. Nearest-neighbour analysis for sources brighter than 
5 mJy, detected at a significance of > S.OOcr over all of the sur- 
vey fields. The vertical dotted line indicates the size of the JCMT 
beam at 850 /Ltm. The solid line shows the distribution of nearest- 
neighbour pairs for the actual dataset, whereas the dashed line 
shows the expected nearest-neighbour histogram for the same sur- 
face density of sources when distributed randomly. The probabil- 
ity that the two distributions are the same is < 0.05. 



Figure 27. Nearest-neighbour analysis for sources brighter than 
6mjy, detected at a significance of > 3.00(t over all of the sur- 
vey fields. The vertical dotted line indicates the size of the JCMT 
beam at 850 /um. The solid line shows the distribution of nearest- 
neighbour pairs for the actual dataset, whereas the dashed line 
shows the expected nearest-neighbour histogram for the same sur- 
face density of sources when distributed randomly. The probabil- 
ity that the two distributions are the same is < 0.05. 



One can then either determine the values of A, 5 and by 
allowing them all to be free parameters in a minimised 

fit, or alternatively one can rewrite the integral constraint 
in terms of the amplitude multiplied by a constant (ie. a = 
A X C) and estimate C by doubly integrating an assumed 
true w(^) over the field area O, 



AC^ 



^41 



w(6»)daidn2 



(21) 



This can be done numerically using the random-random cor- 
relation, such that 



C = 



ENrr(6')6i- 



(22) 



ENrr(e) 

Both of these methods have been attempted for all four 

sets of data points, the free parameter fit proving to be the 
more successful of the two. The results are given in Table 
17. 

The large error bars and scatter in the data points be- 
yond the first 3 fairly robust data point measurements (out 

to ^ 100 arcsec) means that the best-fit parameters for a 
power-law 2-point angular correlation function are not well 



constrained. However, it is interesting to note that in all 
4 cases the characteristic angular scale length 0o indicates 
strong clustering of the bright SCUBA population on scales 
of ~ 40 — 50 arcsec. 

The indication is that a much larger SCUBA survey, 
approaching 0.5 square degrees (4 times the area covered 
here) is required to obtain a meaningful measure of the cor- 
relation power-law slope and scale length. This is discussed 
briefly in Section 8 with more detailed accounts given in van 
Kampen et al. (2005) and Mortier et al. (2006). 



7.2 Nearest neighbour statistics 

An alternative method for measuring the strength of cluster- 
ing is a nearest-neighbour analysis. This procedure measures 

the distribution of the separations between each source and 
its closest neighbour as compared to what one would expect 
from a random distribution, and is sometimes more informa- 
tive in deciding whether sources are clustered when dealing 
with small datasets like these. 

The nearest-neighbour distributions for the bright 
SCUBA sources were measured from the S/N > 3.00 
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Figure 28. Nearest-neighbour analysis for sources brigiiter than 
7mjy, detected at a significance of > S.OOcr over all of the sur- 
vey fields. The vertical dotted line indicates the size of the JCMT 
beam at 850 //m. The solid line shows the distribution of nearest- 
neighbour pairs for the actual dataset, whereas the dashed line 
shows the expected nearest-neighbour histogram for the same sur- 
face density of sources when distributed randomly. The probabil- 
ity that the two distributions are the same is < 0.01. 



datascts, for flux density thresholds of Sgso > 5, > 6, and 
> 7mJy as measured directly from the raw maps (i.e. no 
corrections for flux-boosting effects have been applied). One 
can readily weight any uncharacteristically deep areas of an 
image for calculating a 2-point angular correlation function 
by creating a large mock catalogue of randomly positioned 
sources, equivalent to a combined catalogue from conducting 
many simulations, but which has the local number density 
of detections in any given region weighted according to the 
number counts and the noise as described in the previous 
Section. A nearest neighbour distribution, however, requires 
mock images with the same source number density as the 
actual survey fields. In this analysis 100 mock catalogues 
for each individual field were generated, the positions al- 
located randomly according to a Poisson distribution and 
constrained to lie within the boundaries of the original sur- 
vey data. This was done separately for the three flux density 
thresholds, each fake source list containing the same num- 
ber of > 3.00(7 sources as the real image. Histograms of 
the number of nearest neighbours against nearest-neighbour 
separation are shown in Figs 26 - 28 for flux density cut- 
offs of 5, 6 and 7mJy respectively, using bins the size of 
the SCUBA beam (14.5 arcsec). The histograms of the real 
nearest-neighbour distributions are shown by solid lines, and 
the distributions from combining the mock source catalogues 
arc illustrated by the dashed lines, normalised according to 
the total number of real source pairings. 

In all three plots there is a clear excess of bright 850 /xm 
sources separated from their nearest neighbour at ~ 50 arc- 
sec. This is most noticeable in the Ssso > 7mJy histogram 
(Fig. 28) where there are 51 nearest neighbour pairs com- 
pared to 36 expected within the first ~ 100 arcsec (a sta- 
tistical excess at the 2.5a level), and 34 nearest-neighbour 
pairs compared to 16 oxpoctod within the first ~ 50 arcsec 
(a statistical excess at the 4.5c7 level). 



Flux density 


D 




Number 


P(D) 


threshold 




1 arcsec 


sources 




> 5 mjy 


0.1587 


82.6 


100 


< 0.05 


> 6mjy 


0.1659 


65.5 


78 


< 0.05 


> 7mJy 


0.2898 


65.5 


61 


< 0.01 



Table 18. Results of applying a Kolmogorov-Smirnov test to the 
cumulative fraction of nearest-neighbour pairs. Column 1 gives 
the flux density cutoff as measured from the raw images. Column 
2 gives the maximum absolute difference between the observed 
and expected normalised cunmlativc distributions, D, and the 
nearest-neighbour separation at which this occurs is given in Col- 
umn 3. Column 4 gives the probability that the two distributions 
are the same. 



A Kolmogorov-Smirnov test is a simple method for test- 
ing the probability that two distributions are identical. The 
test statistic, D, is defined as the maximum absolute differ- 
ence between an observed (So) and an expected (Se) nor- 
malised cumulative distribution, in this case applied to the 
cumulative fraction of nearest-neighbour pairs with angular 
separation: 



D : 



|So(e)-Se(e)| 



(23) 



If the measured D value exceeds a critical value when com- 
pared to the known sampling distribution for D appropriate 
to the number of data points, then the two distributions may 
be rejected as being the same at that level of significance. 

As shown in Table 18, the distributions of > 3.00cr 
sources brighter than > 5 and > 6 mJy were rejected as be- 
ing consistent with a random distribution at the 95% con- 
fidence level, and the sources brighter than > 7mJy were 
inconsistent with a random distribution at better than the 
99% confidence level. Overall this again suggests that the 
bright SCUBA population is strongly clustered on scales of 
~ 1 arcmin. 



8 CONCLUSIONS AND FUTURE WORK 

Deep blank-field surveys conducted with SCUBA on the 
JCMT have successfully resolved ~ 30 — 50% of the far- 
infrared oxtragalactic background into discrete sources down 
to the confusion limit of Ssso ^ 2 — 3mJy, with deeper 
surveys making use of gravitational lensing from interven- 
ing massive clusters probing the very faintest submillime- 
tre sources. This has revealed a population of heavily dust- 
enshrouded galaxies at high redshift (z > 1) undergoing a 
burst of massive stax-forming activity. The work presented 
here is an investigation into the nature of the most luminous 
850 /im sources (Sgso > 5mJy), with particular considera- 
tion of their link with the formation and evolution of the 
most massive elliptical galaxies visible in the present-day 
Universe. 

Combining the data from the "SCUBA 8 mJy Survey" 

with other existing blank field surveys, re-reduced and anal- 
ysed in an identical manner, approximately doubles the 
area of sky observed with SCUBA analysed in a homoge- 
neous manner. The combined datasets allow the determina- 
tion of the most accurate source counts to date for Ssso > 
2mJy. The differential source counts for Sgso < 8mJy are 
well satisfied by a power-law of the form dN(> S)/dS = 
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No/(a + S") where the best-fit values of the free param- 
eters wore No = 2.67 x lO", a = 0.49 and a = 3.14. If 
one assumes this is true over the full range of flux densities 
from the very faint to the very bright, this predicts a to- 
tal 850 /um background of 3.8 x 10* mjy deg~^, mid-way be- 
tween the 850 nm extragalactic values of 3.1 x 10* mJy deg~^ 
and 4.4 x 10* mJy deg^^. However, there are some indi- 
cations that the source counts steepen for Ssso > 8mJy. 
This may simply be a result of small number statistics, but 
is possibly indicative of a high-mass cutoff of the SCUBA 
population. A steepening of the source counts would also 
make the bright SCUBA population more prone to the ef- 
fects of gravitational lensing. More physical models of the 
source counts, beginning with a local 850 /itm luminosity 
function, are able to reasonably successfully model the num- 
ber counts using only strong pure luminosity evolution with 
redshift, and assuming a single temperature optically thin 
greybody to represent the thermal spectral energy distri- 
bution of the dust with characteristic dust temperatures in 
the range 30 < Tj < 50 K, and emissivity index in the range 
1.0 < /3 < 2.0. 

Ttentative indications that the bright SCUBA popu- 
lation clusters on scales of 1 — 2 arcmin were first iden- 
tified in the "SCUBA 8mJy Survey" (Scott et al. 2002), 

but these results proved inconclusive due to the small num- 
ber of significant detections. However, the full source cata- 
logue derived from combining all of the survey fields con- 
tains a sufficient number of significant detections to place 
more meaningful constraints on the clustering properties of 
the bright submillimetre population. Measurements of angu- 
lar correlation functions and nearest neighbour statistics for 
Ssso > 5, 6, and 7 mJy, imply strong clustering on scales of 
■~ 40 — 50 arcsec at a significance level of 3.5 — 4(7. Attempts 
have been made to fit the angular correlation function data 
points with a traditional power law of the form w(^) = A0~^ , 
although the small number of sources detected even by com- 
bining all of the existing blank field surveys together mean 
that the values of A and S arc not well constrained. How- 
ever, the data points are in fairly good agreement with the 
clustering strength measured for extremely red objects by 
Daddi et al. (2000), suggesting a strong link between the 
SCUBA and ERO populations. 

The results from this work have provided a major part 
of the motivation for a new wider-field extragalactic submil- 
limetre survey, which has now been underway for three years 
at the JCMT. The aim of this survey is to map 0.5 square de- 
grees of sky, ^ 7 times the total area of the "SCUBA 8mJy 
Survey" to a comparable depth of (Jrms — 2.5 mJy /beam 
at 850 /im. The "SCUBA Half Degree Extragalactic Sur- 
vey (SHADES)" is split over two regions of sky; the Lock- 
man Hole East (a continuation of the 8mJy Survey), and 
the Subaru-XMM Deep Field. Based on the 850 /im source 
counts, a survey of this nature is expected to yield 200 — 300 
sources detected at the > 3.5ct level. This survey also has 
associated deep radio, mid-infrared, near-infrared, optical 
and X-ray imaging. Early example results from SHADES 
are discussed by Dunlop (2005), while full description of 
the SHADES science goals, design and data analysis tech- 
niques is given in Mortier et al. (2005). In van Kampen et al. 
(2005) a variety of theoretical models describing galaxy for- 
mation and evolution are presented with particular emphasis 
on their predictions regarding the redshift distribution and 



clustering properties of the SCUBA population. The first 
results on the sub-millimetre rmmber counts from SHADES 
will be presented by Coppin et al. (2006), while the results 
of a search for radio and mid-infrared identifications in deep 
VLA and Spitzer data will be presented by Ivison et al. (in 
preparation). 
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APPENDIX A: SIMULATIONS BUILDING ON 
THE REAL SURVEY DATA 

Al Differential completeness plots 

The following plots show the percentage of sources recovered 
with signal-to-noise ratio > 3.50 against input flux density 
level for each of the survey fields. Those fields composed 
of a deep pencil beam survey within a wider-area shallower 
survey have had these two components treated separately. 
The plotted curves are the best-fit solutions to the empirical 
functional form 



differential completeness = 100(1 — e 



a(x— b) 



(Al) 



where x is the input flux density, and the values of a and 
b were determined by a minimised fit to the simulation 
results for each 850 yum survey field. The values of a and b 
determined from these fits are given in Table 2. 




Input flux density /mJy 

Figure Al. Percentage of sources recovered against input flux 
density, for the uniform noise regions of the ELAIS N2 field from 
the "SCUBA 8mjy Survey". 




nput flux density /mJy 

Figure A2. Percentage of sources recovered against input flux 

density, for the uniform noise regions of the Lockman Hole East 
wide area field from the "SCUBA 8 mJy Survey" . 
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Figure A3. Percentage of sources recovered against input flux 
density, for the uniform noise regions of the Lockman Hole East 
deep strip from the "SCUBA 8mjy Survey". 
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Figure A4. Percentage of sources recovered against input flux 
density, for the non-uniform noise regions of the ELAIS N2 fleld 
from the "SCUBA 8mjy Survey". 
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Figure A5. Percentage of sources recovered against input flux 
density, for the non-uniform noise regions of the Lockman Hole 
East wide area fleld from the "SCUBA 8 mJy Survey" . 




nput flux density /mJy 

Figure A6. Percentage of sources recovered against input flux 
density, for the uniform noise regions of the 03 hour wide area 
field from the "CUDSS". 




nput flux density /mJy 

Figure A7. Percentage of sources recovered against input flux 
density, for the uniform noise regions of the 03 hour deep area 
from the "CUDSS". 
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Figure A8. Percentage of sources recovered against input flux 
density, for the uniform noise regions of the 10 hour field from 
the "CUDSS". 
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Figure A9. Percentage of sources recovered against input flux 
density, for the uniform noise regions of the 14 hour field from 
the "CUDSS". 
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Figure AlO. Percentage of sources recovered against input flux 
density, for the uniform noise regions of the 22 hour field from 
the 'CUDSS". 
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Figure All. Percentage of sources recovered against input flux 

density, for the uniform noise regions of the SSA13 wide area field 
from the "Hawaii Submm Survey". 
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Figure A12. Percentage of sources recovered against input fiux 
density, for the uniform noise regions of the SSA13 hour deep area 
from the "Hawaii Submm Survey" . 
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Figure A13. Percentage of sources recovered against input flux 
density, for the uniform noise regions of the SSA17 field from the 
"Hawaii Submm Survey". 
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Figure A14. Percentage of sources recovered against input flux 
density, for the uniform noise regions of the SSA22 field from the 
"Hawaii Submm Survey". 
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Figure A15. Percentage of sources recovered against input fiux 
density, for the uniform noise regions of the Lockman Hole deep 
area from the "Hawaii Submm Survey". 
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Figure A16. Percentage of sources recovered against input flux 
density, for the uniform noise regions of the Hubble Deep Field. 

A2 Output versus Input Flux Density 

The following plots show the ratio of output to input flux 
density of sources recovered with signal-to-noise ratio > 3.50 
against input flux density level for each of the survey fields. 
Those fields composed of a deep pencil beam survey within 
a wider-area shallower survey have had these two compo- 
nents treated separately. The plotted curves are the best-fit 
solutions to the empirical functional form 

output flux density ^ Qg-dx _|_ ^ /^2) 
input flux density 

where x is the input flux density, and the values of C, d and 

f were determined by a minimised fit to the simulation 
results for each 850 /xm survey field and are given in Table 
3. 
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Figure A19. The ratio of output to input flux density against 
input flux density, for the uniform noise regions of the Lockman 
Hole East deep strip from the "SCUBA 8 mJy Survey" . 
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Figure A20. The ratio of output to input flux density against 
input flux density, for the non-uni noise regions of the ELAIS N2 
fleld from the "SCUBA 8mjy Survey". 
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Figure A17. The ratio of output to input flux density against 
input flux density, for tlic uniform noise regions of the ELAIS N2 
field from the "SCUBA 8mJy Survey". 
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Figure A21. The ratio of output to input flux density against 
input flux density, for the non-uni noise regions of the Lockman 
Hole East wide area field from the "SCUBA 8 mJy Survey" . 
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Figure A18. The ratio of output to input flux density against 
input flux density, for the uniform noise regions of the Lockman 
Hole East wide area field from the "SCUBA 8 mJy Survey" . 
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Figure A22. The ratio of output to input flux density against 
input flux density, for the uniform noise regions of the 03 hour 
wide area fleld from the "CUDSS". 
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Figure A23. The ratio of output to input flux density against 
input flux density, for the uniform noise regions of the 03 hour 
deep area from the "CUDSS". 
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Figure A24. Percentage of sources recovered against input flux 
density, for tire uniform noise regions of tlie 10 iiour field from 
the "CUDSS". 
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Figure A27. The ratio of output to input flux density against 
input flux density, for the uniform noise regions of the SSA13 
wide area field from the "Hawaii Submm Survey" . 
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Figure A28. The ratio of output to input flux density against 
input flux density, for the uniform noise regions of the SSA13 
hour deep area from the "Hawaii Submm Survey" . 
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Figure A25. The ratio of output to input flux density against 
input flux density, for the uniform noise regions of the 14 hour 
field from the "CUDSS". 
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Figure A26. The ratio of output to input flux density against 
input flux density, for the uniform noise regions of the 22 hour 
field from the 'CUDSS". 
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Figure A29. The ratio of output to input flux density against 
input flux density, for the uniform noise regions of the SSA17 field 
from the "Hawaii Submm Survey" . 
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Figure A30. The ratio of output to input flux density against 
input flux density, for the uniform noise regions of the SSA22 field 
from the "Hawaii Submm Survey" . 
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Figure A31. The ratio of output to input flux density against 
input flux density, for the uniform noise regions of the Lockman 
Hole deep area from the "Hawaii Submm Survey" . 
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Figure A32. The ratio of output to input flux density against 
input flux density, for the uniform noise regions of the Hubble 
Deep Field. 



A3 Positional Uncertainty 

The following plots show the mean positional uncertainty of 
sources recovered with signal-to-noise ratio > 3.50 against 
input flux density level for each of the survey fields. Those 
fields composed of a deep pencil beam survey within a 
wider-area shallower survey have had these two components 
treated separately. The mean positional uncertainty in re- 
trieving the fake sources was found to be well approximated 
by a linear dependence on the input flux density such that 



positional error = — gx -I- h 



(A3) 



where x is the input flux density, the values of g and h for 
each 850 nm survey field were determined by a minimised 
fit to the simulation results (given in Table 4), and the 
positional uncertainty is given in arcseconds. Note these sim- 
ulations do not account for pointing drifts whilst observing. 



Figure A33. Mean positional uncertainty against input flux den- 
sity, for the uniform noise regions of the ELAIS N2 field from the 
"SCUBA 8mJy Survey". 
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Figure A34. Mean positional unecrtainty against input flux den- 
sity, for the uniform noise regions of the Lockman Hole East wide 
area field from the "SCUBA 8mjy Survey". 
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Figure A35. Mean positional uncertainty against input flux den- 
sity, for the uniform noise regions of the Lockman Hole East deep 
strip from the "SCUBA 8mJy Survey". 
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Figure A36. Mean positional uncertainty against input flux den- 
sity, for the non-uniform noise regions of the ELAIS N2 field from 
the "SCUBA 8mjy Survey". 
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Figure A38. Mean positional uncertainty against input flux den- 
sity, for the uniform noise regions of the 03 hour wide area field 
from the "CUDSS". 
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Figure A42. Mean positional uncertainty against input flux den- 
sity, for the uniform noise regions of the 22 hour field from the 
'CUDSS". 
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Figure A39. Mean positional uncertainty against input flux den- 
sity, for the uniform noise regions of the 03 hour deep area from 
the "CUDSS". 




Figure A43. Mean positional uncertainty against input flux den- 
sity, for the uniform noise regions of the SSA13 wide area field 
from the "Hawaii Submm Survey" . 
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Figure A40. Mean positional uncertainty against input flux den- 
sity, for the uniform noise regions of the 10 hour field from the 
"CUDSS". 
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Figure A44. Mean positional uncertainty against input flux den- 
sity, for the uniform noise regions of the SSA13 hour deep area 
from the "Hawaii Submm Survey" . 
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Figure A45. Mean positional uncertainty against input flux den- 
sity, for the uniform noise regions of the SSA17 field from the 
"Hawaii Submm Survey". 
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Figure A46. Mean positional uncertainty against input flux den- 
sity, for the uniform noise regions of the SSA22 field from the 
"Hawaii Submm Survey". 
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Figure A47. Mean positional uncertainty against input flux den- 
sity, for the uniform noise regions of the Lockman Hole deep area 
from the "Hawaii Submm Survey" . 
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Figure A48. Mean positional uncertainty against input flux den- 
sity, for the uniform noise regions of the Hubble Deep Field. 
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APPENDIX B: COMPLETELY SIMULATED MAPS 
Bl Integral completeness tables 

The following tables give the percentage integral completeness for each of the survey fields over the flux density range 
2 — lOmJy, and for signal-to-noise thresholds of > 1.50cr to > 4.00cr. Those fields composed of a deep pencil beam survey 
within a wider-area shallower survey have had these two components treated separately. 



ELAIS N2 field, uniform noise region 


Flux 
Density 








% Integral Completeness 








> 1.50cr 


> 2.00(7 


> 2.50cr 


> 3.00(7 


> 3.50(7 


> 4.00(7 


> 2mJy 


45.7 ± 


0.5 


36.3 ± 0.4 


25.6 ± 0.3 


17.6 ± 0.2 


11.5 ± 


0.1 


7.6 ± 0.1 


> 4mJy 


71.2 ± 


1.7 


68.9 ± 1.6 


60.8 ± 1.4 


50.9 ± 1.2 


38.7 ± 


0.9 


28.8 ± 0.7 


> 6mjy 


79.9 ± 


3.1 


79.6 ± 3.1 


77.3 ± 3.0 


72.8 ± 2.8 


64.7 ± 


2.5 


55.4 ± 2.2 


> 8mJy 


76.4 ± 


4.4 


76.4 ± 4.4 


76.1 ± 4.4 


74.4 ± 4.3 


71.1 ± 


4.1 


65.2 ± 3.7 


> 10 mjy 


76.8 ± 


6.4 


76.8 ± 6.4 


76.8 ± 6.4 


76.8 ± 6.4 


74.6 ± 


6.3 


71.1 ± 6.0 



Table Bl. Percentage integral completeness results for the uniform noise region of the ELAIS N2 field, over the flux density range 
2 — lOmJy, and for signal-to-noise thresholds of > 1.50(7 to > 4.00(7. 



ELAIS N2 field, non-uniform noise region 



Flux % Integral Completeness 

Density > 1.50a > 2.00<t > 2.50a > S.OOct > 3.50a > i.OOa 



> 2 mJy 


12.6 ± 


0.3 


8.2 ± 


0.2 


4.7 ± 


0.1 


2.4 ± 


0.1 


1.0 ± 


0.1 


0.4 ± 


0.1 


> 4mJy 


24.2 ± 


1.1 


18.3 ± 


0.8 


12.7 ± 


0.6 


7.7 ± 


0.3 


4.0 ± 


0.2 


1.9 ± 


0.1 


> 6mjy 


35.6 ± 


2.7 


31.6 ± 


2.4 


24.9 ± 


1.9 


14.7 ± 


1.1 


9.6 ± 


0.7 


5.6 ± 


0.4 


> 8mJy 


37.8 ± 


4.4 


33.8 ± 


3.9 


31.1 ± 


3.6 


23.0 ± 


2.7 


14.9 ± 


1.7 


10.8 ± 


1.3 


> lOmJy 


38.1 ± 


5.9 


35.7 ± 


5.5 


33.3 ± 


5.1 


28.6 ± 


4.4 


19.0 ± 


2.9 


11.9 ± 


1.8 



Table B2. Percentage integral completeness results for the non-uniform noise region of the ELAIS N2 field, over the flux density range 
2 — lOmJy, and for signal-to-noise thresholds of > 1.50<T to > 4.00<T. 



Lockman Hole wide area field, uniform noise region 



Flux % luk\gral ( '()m])l(;kuu^bs 

Density > 1.50(7 > 2.00(7 > 2.50(7 > 3.00(7 > 3.50(7 > 4.00(7 



> 2mJy 


41.8 ± 


0.5 


30.0 ± 


0.4 


19.2 ± 


0.2 


11.9 ± 


0.1 


7.1 ± 


0.1 


4.3 ± 


0.1 


> 4mJy 


72.0 ± 


2.0 


65.3 ± 


1.8 


52.7 ± 


1.4 


38.8 ± 


1.1 


26.9 ± 


0.7 


18.2 ± 


0.5 


> 6mjy 


78.1 ± 


3.5 


77.6 ± 


3.4 


74.4 ± 


3.3 


64.0 ± 


2.8 


51.4 ± 


2.3 


38.2 ± 


1.7 


> 8mJy 


81.3 ± 


5.4 


81.3 ± 


5.4 


80.9 ± 


5.4 


79.6 ± 


5.3 


73.8 ± 


4.9 


60.9 ± 


4.1 


> 10 mJy 


73.0 ± 


6.8 


73.0 ± 


6.8 


73.0 ± 


6.8 


73.0 ± 


6.8 


72.2 ± 


6.7 


67.0 ± 


6.2 



Table B3. Percentage integral completeness results for the uniform noise region of the Lockman Hole wide area field, over the fiux 
density range 2 — lOmJy, and for signal-to-noise thresholds of > 1.50<t to > A.OOa. 



Lockman Hole wide area field, non-uniform noise region 



Flux % Integral Completeness 

Density > 1.50a > 2.00ct > 2.50a > 3.00ct > 3.50a > 4.00(7 



> 2mJy 


19.2 ± 


0.3 


13.1 ± 


0.2 


7.9 ± 


0.1 


3.8 ± 


0.1 


1.8 ± 


0.1 


0.7 ± 


0.1 


> 4mJy 


34.6 ± 


1.2 


25.7 ± 


0.9 


18.0 ± 


0.6 


9.8 ± 


0.4 


5.4 ± 


0.2 


2.5 ± 


0.1 


> 6mjy 


39.6 ± 


2.3 


34.0 ± 


2.0 


25.4 ± 


1.5 


16.2 ± 


0.9 


8.6 ± 


0.5 


4.0 ± 


0.2 


> 8 mJy 


36.8 ± 


3.1 


34.7 ± 


2.9 


29.9 ± 


2.5 


19.4 ± 


1.6 


10.4 ± 


0.9 


4.9 ± 


0.4 


> lOmJy 


30.3 ± 


3.5 


30.3 ± 


3.5 


28.9 ± 


3.3 


21.1 ± 


2.4 


14.5 ± 


1.7 


7.9 ± 


0.9 



Table B4. Percentage integral completeness results for the non-uniform noise region of the Lockman Hole wide area field, over the fiux 
density range 2 — 10 mJy, and for signal-to-noise thresholds of > 1.50cr to > 4.00<T. 
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Lockman Hole deep strip field, uniform noise region 


Flux 


% Integral Completeness 




Density 


> l.SOo- > 2.00CT > 2.50(7 > 3.00CT > 3.50(t 


> 4.00(T 



> 2mJy 57.2 ± 1.5 48.6 ± 1.2 38.0 ± 1.0 29.6 ± 0.8 21.7 ± 0.6 15.0 ± 0.4 

> 4mJy 82.3 ± 4.7 82.3 ± 4.7 80.0 ± 4.6 76.1 ± 4.4 63.3 ± 3.6 50.2 ± 2.9 

> 6mJy 85.2 ± 7.9 85.2 ± 7.9 85.2 ± 7.9 85.2 ± 7.9 83.5 ± 7.8 78.3 ± 7.3 

> 8mJy 84.6 ±11.7 84.6 ± 11.7 84.6 ± 11.7 84.6 ± 11.7 84.6 ± 11.7 84.6 ±11.7 
>10mJy 65.4 ±12.8 65.4 ± 12.8 65.4 ± 12.8 65.4 ± 12.8 65.4 ± 12.8 65.4 ± 12.8 



Table B5. Percentage integral completeness results for the uniform noise region of the Lockman Hole deep strip, over the flux density 
range 2 — lOmJy, and for signal-to-noise thresholds of > 1.50(T to > 4.00(7. 



03h deep area field, uniform noise region 


Flux 




% Integral Completeness 




Density 


> 1.50cr 


> 2.00(7 > 2.50(r > 3.00(7 > 3.50(7 


> 4.00(7 



> 2mJy 58.7 ± 2.6 58.5 ± 2.5 57.2 ± 2.5 53.8 ± 2.3 46.6 ± 2.0 38.4 ± 1.7 

> 4mJy 82.2 ± 7.6 82.2 ± 7.6 82.2 ± 7.6 82.2 ± 7.6 82.2 ± 7.6 82.2 ± 7.6 

> 6mJy 88.6 ±13.4 88.6 ± 13.4 88.6 ± 13.4 88.6 ± 13.4 88.6 ± 13.4 88.6 ± 13.4 

> 8mJy 82.6 ±17.2 82.6 ± 17.2 82.6 ±17.2 82.6 ± 17.2 82.6 ± 17.2 82.6 ± 17.2 
>10mJy 90.0 ±28.5 90.0 ± 28.5 90.0 ± 28.5 90.0 ± 28.5 90.0 ± 28.5 90.0 ±28.5 



Table B6. Percentage integral completeness results for the uniform noise region of the 03h deep area field, over the flux density range 
2 — lOmJy, and for signal-to-noise thresholds of > 1.50(7 to > 4.00(7. 



03h wide area field, uniform noise region 


Flux 




% Integral Completeness 




Density 


> 1.50(7 


> 2.00(7 > 2.50(T > 3.00(7 > 3.50(7 


> 4.00(7 



> 2mJy 43.8 ± 0.6 35.5 ± 0.5 26.5 ± 0.4 19.2 ± 0.3 13.1 ± 0.2 8.9 ± 0.1 

> 4mJy 68.4 ± 2.2 65.8 ± 2.1 60.4 ± 2.0 52.9 ± 1.7 42.3 ± 1.4 32.6 ± 1.1 

> 6mJy 73.8 ± 4.0 73.2 ± 4.0 70.9 ± 3.8 66.8 ± 3.6 60.3 ± 3.3 52.4 ± 2.8 

> 8mJy 71.5 ± 6.0 71.5 ± 6.0 71.5 ± 6.0 71.5 ± 6.0 68.8 ± 5.7 63.9 ± 5.3 
>10mJy 69.1 ± 8.4 69.1 ± 8.4 69.1 ± 8.4 69.1 ± 8.4 69.1 ± 8.4 67.6 ± 8.2 



Table B7. Percentage integral completeness results for the uniform noise region of the 03h wide area field, over the flux density range 
2 — lOmJy, and for signal-to-noise thresholds of > 1.50(7 to > 4.00c7. 



lOh field, uniform noise region 


Flux 
Density 






Vu iiiLcgial (.'oiiii>ioLciicss 






> 1.50(7 


> 2.00(7 


> 2.50(7 


> 3.00(7 


> 3.50(7 


> 4.00(7 


> 2mJy 


44.4 ± 1.7 


39.0 ± 1.5 


33.6 ± 1.3 


26.5 ± 1.0 


19.4 ± 0.8 


14.8 ± 0.6 


> 4mJy 


57.1 ± 5.0 


56.4 ± 4.9 


55.6 ± 4.8 


52.6 ± 4.6 


50.4 ± 4.4 


47.4 ± 4.1 


> 6mjy 


68.1 ± 9.9 


68.1 ± 9.9 


68.1 ± 9.9 


66.0 ± 9.6 


66.0 ± 9.6 


63.8 ± 9.3 


> 8mJy 


70.0 ± 15.7 


70.0 ±15.7 


70.0 ± 15.7 


70.0 ± 15.7 


70.0 ±15.7 


70.0 ± 15.7 


> 10 mjy 


72.7 ±21.9 


72.7 ±21.9 


72.7 ±21.9 


72.7 ±21.9 


72.7 ±21.9 


72.7 ±21.9 



Table B8. Percentage integral completeness results for the uniform noise region of the lOh field, over the flux density range 2 — lOmJy, 
and for signal-to-noise thresholds of > 1.50(7 to > 4.00(7. 









14h field 


, uniform noise region 










Flux 
Density 








% Integral Completeness 










> 1.50(7 


> 2.00(7 


> 2.50(7 


> 3.00(7 


> 3.50(7 


> 4.00(7 


> 2mJy 


51.0 ± 


0.8 


45.1 ± 0.7 


36.2 ± 0.5 


27.5 ± 0.4 


20.6 ± 


0.3 


15.3 ± 


0.2 


> 4mJy 


70.3 ± 


2.3 


69.1 ± 2.2 


65.9 ± 2.1 


61.5 ± 2.0 


55.4 ± 


1.8 


48.8 ± 


1.6 


> 6 mJy 


74.5 ± 


3.9 


74.2 ± 3.9 


73.4 ± 3.8 


70.9 ± 3.7 


69.0 ± 


3.6 


67.4 ± 


3.5 


> 8 mJy 


77.8 ± 


6.0 


77.8 ± 6.0 


77.8 ± 6.0 


76.0 ± 5.9 


75.4 ± 


5.8 


74.9 ± 


5.8 


> lOmJy 


78.3 ± 


9.4 


78.3 ± 9.4 


78.3 ± 9.4 


76.8 ± 9.2 


76.8 ± 


9.2 


75.4 ± 


9.1 



Table B9. Percentage integral completeness results for the uniform noise region of the 14h field, over the flux density range 2 — 10 mJy, 
and for signal-to-noise thresholds of > 1.50(7 to > 4.00(7. 
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22h field, uniform noise region 


Flux 




% Integral Completeness 






Density 


> 1.50(T 


> 2.00CT > 2.50(7 > 3.00CT 


> 3.50(T 


> 4.00(T 



> 2mJy 49.9 ± 2.4 45.1 ± 2.2 36.8 ± 1.8 27.7 ± 1.4 20.5 ± 1.0 14.3 ± 0.7 

> 4mJy 72.2 ± 8.1 72.2 ± 8.1 72.2 ± 8.1 72.2 ± 8.1 65.8 ± 7.4 55.7 ± 6.3 

> 6mJy 80.0 ±16.0 80.0 ± 16.0 80.0 ± 16.0 80.0 ± 16.0 80.0 ± 16.0 80.0 ± 16.0 

> 8mJy 60.0 ±15.5 60.0 ± 15.5 60.0 ± 15.5 60.0 ± 15.5 60.0 ± 15.5 60.0 ± 15.5 

> 10 mJy 75.0 ± 37.5 75.0 ± 37.5 75.0 ± 37.5 75.0 ± 37.5 75.0 ± 37.5 75.0 ± 37.5 



Table BIO. Percentage integral completeness results for the uniform noise region of the 22h field, over the flux density range 2 — 10 mJy, 
and for signal-to-noise thresholds of > 1.50a to > 4.00cr. 





SSA13 deep 


area field, uniform noise region 




Flux 




% Integral Completeness 




Density 


> 1.50(T > 2.00(7 


> 2.50(T > 3.00(7 > 3.50(7 


> 4.00(7 



> 2mJy 51.3 ± 2.0 50.6 ± 1.9 48.4 ± 1.9 47.1 ± 1.8 42.9 ± 1.6 38.8 ± 1.5 

> 4mJy 67.9 ± 5.9 67.9 ± 5.9 67.2 ± 5.8 67.2 ± 5.8 67.2 ± 5.8 67.2 ± 5.8 

> 6mJy 84.1 ±12.7 84.1 ± 12.7 84.1 ± 12.7 84.1 ± 12.7 84.1 ± 12.7 84.1 ± 12.7 

> 8mJy 68.0 ±13.6 68.0 ± 13.6 68.0 ± 13.6 68.0 ± 13.6 68.0 ± 13.6 68.0 ± 13.6 
>10mJy 75.0 ±21.7 75.0 ±21.7 75.0 ± 21.7 75.0 ±21.7 75.0 ±21.7 75.0 ±21.7 



Table Bll. Percentage integral completeness results for the uniform noise region of the SSA13 deep area field, over the flux density 
range 2 — lOmJy, and for signal-to-noise thresholds of > 1.50(7 to > 4.00(7. 



SSA13 wide area field, uniform noise region 


Flux 


% Integral Completeness 




Density 


> 1.50(7 > 2.00(7 > 2.50(T > 3.00(7 > 3.50(7 


> 4.00(7 



> 2mJy 38.6 ± 0.6 29.2 ± 0.5 21.0 ± 0.3 14.2 ± 0.2 9.4 ± 0.2 6.2 ± 0.1 

> 4mJy 61.2 ± 2.1 56.6 ± 1.9 49.5 ± 1.7 39.3 ± 1.4 29.2 ± 1.0 21.8 ± 0.7 

> 6mJy 70.7 ± 3.9 67.0 ± 3.7 63.6 ± 3.5 59.6 ± 3.3 50.6 ± 2.8 42.0 ± 2.3 

> 8mJy 77.0 ± 6.5 76.3 ± 6.5 73.4 ± 6.2 72.7 ± 6.2 66.2 ± 5.6 60.4 ± 5.1 
>10mJy 73.0± 9.2 73.0 ± 9.2 73.0 ± 9.2 73.0 ± 9.2 69.8 ± 8.8 65.1 ± 8.2 



Table B12. Percentage integral completeness results for the uniform noise region of the SSA13 wide area field, over the flux density 
range 2 — lOmJy, and for signal-to-noise thresholds of > 1.50(7 to > 4.00(7. 



SSA17 field, uniform noise region 


Flux 




'/o iliLcgial (.'oiiii)iclciK;ss 






Density 


> 1.50(7 


> 2.00(7 > 2.50(7 > 3.00(7 


> 3.50(7 


> 4.00(7 



> 2mJy 


49.6 ± 


1.2 


42.7 ± 


1.1 


35.3 ± 


0.9 


29.2 ± 


0.7 


22.9 ± 


0.6 


18.3 ± 


0.5 


> 4mJy 


70.9 ± 


3.6 


69.0 ± 


3.6 


68.3 ± 


3.5 


65.1 ± 


3.3 


57.1 ± 


2.9 


50.3 ± 


2.6 


> 6mjy 


78.9 ± 


6.5 


78.9 ± 


6.5 


78.9 ± 


6.5 


78.2 ± 


6.5 


74.1 ± 


6.1 


67.3 ± 


5.6 


> 8mJy 


71.8 ± 


8.5 


71.8 ± 


8.5 


71.8 ± 


8.5 


71.8 ± 


8.5 


71.8 ± 


8.5 


71.8 ± 


8.5 


> 10 mJy 


69.7 ± 


12.1 


69.7 ± 


12.1 


69.7 ± 


12.1 


69.7 ± 


12.1 


69.7 ± 


12.1 


69.7 ± 


12.1 



Table B13. Percentage integral completeness results for the uniform noise region of the SSA17 field, over the flux density range 2 — 10 mJy, 
and for signal-to-noise thresholds of > 1.50(7 to > 4.00(7. 







SSA22 field, uniform : 


noise region 






Flux 






% Integral Completeness 






Density 


> 1.50(7 


> 2.00(7 


> 2.50(7 


> 3.00(7 


> 3.50(7 


> 4.00(7 


> 2mJy 

> 4mJy 

> 6 mJy 

> 8 mJy 

> lOmJy 


50.3 ± 1.2 
65.7 ± 3.6 
73.9 ± 6.8 
69.6 ± 9.3 
78.3 ± 16.3 


49.0 ± 1.2 
65.4 ± 3.6 
73.9 ± 6.8 
69.6 ± 9.3 
78.3 ± 16.3 


46.3 ± 1.2 

64.8 ± 3.5 

73.9 ± 6.8 
69.6 ± 9.3 
78.3 ± 16.3 


43.1 ± 1.1 

64.2 ± 3.5 
73.1 ± 6.7 
69.6 ± 9.3 

78.3 ± 16.3 


39.0 ± 1.0 
62.7 ± 3.4 

73.1 ± 6.7 

69.6 ± 9.3 
78.3 ± 16.3 


34.6 ± 0.9 
60.9 ± 3.3 
73.1 ± 6.7 

69.6 ± 9.3 
78.3 ±16.3 



Table B14. Percentage integral completeness results for the uniform noise region of the SSA22 field, over the flux density range 2 — lOmJy, 
and for signal-to-noise thresholds of > 1.50(7 to > 4.00(7. 
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Lockman Hole deep field, uniform noise region 


Flux 


% Integral Completeness 




Density 


> 1.50(7 > 2.00CT > 2.50(7 > 3.00CT > 3.50(t 


> 4.00(7 



> 2mJy 45.9 ± 1.8 45.5 ± 1.8 43.6 ± 1.7 41.8 ± 1.6 37.9 ± 1.5 33.6 ± 1.3 

> 4mJy 60.4 ± 5.0 60.4 ± 5.0 60.4 ± 5.0 60.4 ± 5.0 59.7 ± 5.0 59.0 ± 4.9 

> 6mJy 68.5 ± 9.3 68.5 ± 9.3 68.5 ± 9.3 68.5 ± 9.3 68.5 ± 9.3 68.5 ± 9.3 

> 8mJy 64.3 ±12.1 64.3 ±12.1 64.3 ±12.1 64.3 ±12.1 64.3 ±12.1 64.3 ±12.1 
>10mJy 70.6 ±17.1 70.6 ±17.1 70.6 ± 17.1 70.6 ±17.1 70.6 ±17.1 70.6 ±17.1 



Table B15. Percentage integral completeness results for the uniform noise region of the Lockman Hole deep field, over the flux density 
range 2 — lOmJy, and for signal-to-noise thresholds of > 1.50(7 to > 4.00(7. 







Hubble deep field, uniform 


noise region 






Flux 






% Integral Completeness 






Density 


> 1.50cr 


> 2.00(7 


> 2.50(T 


> 3.00(7 


> 3.50(7 


> 4.00(7 


> 2 mJy 

> 4 mJy 

> 6mjy 

> 8mJy 

> 10 mJy 


61.6 ± 2.9 

77.3 ± 7.9 
93.0 ± 14.2 

82.4 ± 20.0 
100 ±40.8 


61.6 ± 2.9 

77.3 ± 7.9 
93.0 ± 14.2 

82.4 ±20.0 
100 ±40.8 


61.6 ± 2.9 

77.3 ± 7.9 
93.0 ± 14.2 

82.4 ±20.0 
100 ±40.8 


61.3 ± 2.9 

77.3 ± 7.9 
93.0 ± 14.2 

82.4 ±20.0 
100 ±40.8 


59.8 ± 2.8 

77.3 ± 7.9 
93.0 ± 14.2 

82.4 ±20.0 
100 ±40.8 


57.1 ± 2.7 

77.3 ± 7.9 
93.0 ± 14.2 

82.4 ± 20.0 
100 ±40.8 



Table B16. Percentage integral completeness results for the uniform noise region of the Hubble deep field, over the flux density range 
2 — lOmJy, and for signal-to-noise thresholds of > 1.50(7 to > 4.00(7. 



B2 Count correction tables 

The following tables give the percentage count correction factors for each of the survey fields over the flux density range 
2 — lOmJy, and for signal-to-noise thresholds of > 1.50(7 to > 4.00(7. Those fields composed of a deep pencil beam survey 
within a wider-area shallower survey have had these two components treated separately. 



ELAIS N2 field, uniform noise region 



Flux %■ Count Correction 



Deiisil\' 


> 1.50(T 


> 2.()0a 


> 2..")0(T 


> :',ma 


> :i.50(T 


> l.()0(T 


> 2 mJy 


44.0 ± 0.3 


78.9 ± 0.7 


160.5 ± 


2.2 


3.34.8 ± 6.5 


652.6 ± 17.7 


> 1000 


> 4 m Jy 


14.8 ± 0.1 


20.8 ± 0.2 


36.5 ± 


0.5 


71.1 ± 1.4 


137.3 ± 3.7 


245.5 ± 8.9 


> 6 mJy 


13.3 ± 0.2 


16.6 ± 0.3 


24.8 ± 


0.5 


38.5 ± 0.9 


61.0 ± 1.8 


94.3 ± 3.6 


> 8mJy 


17.1 ± 0.4 


18.7 ± 0.5 


25.4 ± 


0.7 


37.9 ± 1.3 


54.7 ± 2.3 


75.7 ± 3.8 


> 10 mJy 


26.4 ± 1.1 


27.1 ± 1.2 


33.0 ± 


1.6 


45.1 ± 2.5 


59.7 ± 3.9 


75.9 ± 5.6 



Table B17. Percentage count correction results for the uniform noise region of the ELAIS N2 field, over the flux density range 2 — 10 mJy, 
and for signal-to-noise thresholds of > 1.50(7 to > 4.00(7. 



ELAIS N2 field, non-uniform noise region 



Flux 










% Count Correction 






Density 


> 1.50 


(7 


> 2.00 


(7 


> 2.50(7 


> 3.00(7 


> 3.50(7 


> 4.00(7 


> 2mJy 


83.0 ± 


1.5 


158.4 ± 


4.0 


467.8 ±20.5 


> 1000 


> 1000 


> 1000 


> 4mJy 


17.8 ± 


0.3 


33.9 ± 


0.9 


100.2 ± 4.4 


366.2 ±30.7 


> 1000 


> 1000 


> 6mjy 


6.8 ± 


0.1 


11.6 ± 


0.3 


34.1 ± 1.5 


124.6 ±10.5 


491.7 ±81.9 


> 1000 


> 8mJy 


4.0 ± 


0.1 


5.8 ± 


0.2 


14.7 ± 0.7 


52.1 ± 4.4 


205.6 ± 34.3 


616.7 ± 178.0 


> 10 mJy 


4.7 ± 


0.2 


6.6 ± 


0.3 


14.6 ± 0.9 


42.0 ± 4.2 


120.0 ±20.3 


350.0 ± 101.0 



Table BIS. Percentage count correction results for the non-uniform noise region of the ELAIS N2 field, over the flux density range 
2 — lOmJy, and for signal-to-noise thresholds of > 1.50(7 to > 4.00(7. 
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Lockman Hole wide area field, uniform noise region 



Flux % Count Correction 



Density 


> 1.50CT 


> 2.00CT 


> 2.50CT 


> 3.00CT 


> 3.50CT 


> 4.00CT 


> 2mJy 


45.7 ± 0.4 


93.8 ± 1.1 


207.2 ± 3.7 


459.8 ± 12.3 


> 1000 


> 1000 


> 4 mjy 


11.4 ± 0.1 


20.7 ± 0.3 


44.1 ± 0.8 


96.6 ± 2.6 


230.5 ± 9.5 


448.0 ± 25.8 


> 6 mJy 


9.7± 0.1 


12.3 ± 0.2 


20.9 ± 0.4 


41.2 ± 1.2 


91. 7± 3.9 


171.0 ± 9.9 


> 8mJy 


12.3 ± 0.3 


13.5 ± 0.3 


18.1 ± 0.5 


26.6 ± 0.9 


52.2 ± 2.5 


88.6 ± 5.6 


> 10 mJy 


22.1 ± 1.0 


22.7 ± 1.0 


28.8 ± 1.4 


36.4 ± 2.0 


46.9 ± 3.0 


65.0 ± 4.9 



Table B19. Percentage count correction results for the uniform noise region of the Lockman Hole wide area field, over the flux density 
range 2 — lOmJy, and for signal-to-noise thresholds of > 1.50(T to > 4.00iT. 



Lockman Hole wide area fleld, non-uniform noise region 



Flux % Count Correction 



Density 


> 1.50 


cr 


> 2.00 


cr 


> 2.50 


(T 


> 3.00(7 


> 3.50(7 


> 4.00(7 


> 2 mJy 


59.3 d= 


0.8 


102.7 ± 


1.7 


191.8 ± 


4.4 


441.5 ± 15.2 


> 1000 


> 1000 


> 4 mJy 


12.6 ± 


0.2 


21.8 ± 


0.4 


40.7 ± 


0.9 


93.6 ± 3.2 


312.7 ±19.7 


> 1000 


> 6mjy 


5.3 ± 


0.1 


8.5 ± 


0.1 


15.7 ± 


0.4 


36.1 ± 1.2 


120.7 ± 7.6 


673.3 ± 100.4 


> 8mJy 


3.1 ± 


0.1 


4.4 ± 


0.1 


7.6 ± 


0.2 


17.2 ± 0.6 


57.4 ± 3.6 


320.0 ± 47.7 


> 10 mJy 


2.3 ± 


0.1 


2.9 ± 


0.1 


4.5 ± 


0.1 


9.8 ± 0.4 


31.4 ± 2.0 


172.7 ± 26.0 



Table B20. Percentage count correction results for the nonuniform noise region of the Lockman Hole wide area field, over the flux 
density range 2 — lOmJy, and for signal-to-noise thresholds of > 1.50(7 to > 4.00(7. 



Lockman Hole deep strip field, uniform noise region 



Flux 






% Count Correction 








Density 


> 1.50CT 


> 2.00(7 


> 2.50(7 


> 3.00(7 


> 3.50(7 


> 4.00(7 


> 2mJy 


41.5 ± 0.7 


65.2 ± 1.3 


113.3 ± 3.1 


205.5 ± 7.5 


352.7 ± 16.9 


606.0 ± 


38.2 


> 4mJy 


26.5 ± 0.8 


27.2 ± 0.8 


32.5 ± 1.1 


45.3 ± 1.7 


71.6 ± 3.5 


121.0 ± 


7.6 


> 6mjy 


50.2 ± 3.3 


50.2 ± 3.3 


50.2 ± 3.3 


51.8 ± 3.5 


56.1 ± 3.9 


63.9 ± 


4.8 


> 8mJy 


83.9 ± 10.7 


83.9 ± 10.7 


83.9 ± 10.7 


83.9 ± 10.7 


83.9 ± 10.7 


83.9 ± 


10.7 


> lOmJy 


81.2 ± 14.4 


81.2 ± 14.4 


81.2 ±14.4 


81.2 ±14.4 


81.2 ±14.4 


81.2 ± 


14.4 



Table B21. Percentage count correction results for the uniform noise region of the Lockman Hole deep strip field, over the flux density 
range 2 — lOmJy, and for signal-to-noise thresholds of > 1.50(7 to > 4.00(7. 



03h deep area field, uniform noise region 



Flux 






'a C-'ouiiL 


(.'orrcxLioii 






Density 


> 1.50(7 


> 2.00(7 


> 2.50(7 


> 3.00(7 


> 3.50(7 


> 4.00(7 


> 2mJy 


62.3 ± 2.1 


64.9 ± 2.3 


76.5 ± 2.9 


100.4 ± 4.4 


140.1 ± 7.2 


199.2 ± 12.2 


> 4mJy 


58.7 ± 4.1 


59.9 ± 4.3 


62.8 ± 4.6 


64.8 ± 4.8 


67.8 ± 5.1 


72.8 ± 5.7 


> 6mjy 


62.9 ± 7.5 


63.8 ± 7.7 


63.8 ± 7.7 


64.7 ± 7.8 


65.7 ± 8.0 


65.7 ± 8.0 


> 8mJy 


88.5 ± 17.3 


88.5 ± 17.3 


88.5 ± 17.3 


92.0 ±18.4 


92.0 ± 18.4 


92.0 ± 18.4 


> 10 mJy 


71.4 ±19.1 


71.4 ±19.1 


71.4 ± 19.1 


76.9 ± 21.3 


76.9 ± 21.3 


76.9 ± 21.3 



Table B22. Percentage count correction results for the uniform noise region of the 03h deep area field, over the flux density range 
2 — lOmJy, and for signal-to-noise thresholds of > 1.50(7 to > 4.00(7. 



03h wide area field, uniform noise region 



Flux % Count Correction 

Density > 1.50(T > 2.00a > 2.50a > 3.00(7 > 3.50(7 > 4.00(7 



> 2mJy 


52.8 ± 


0.6 


89.2 ± 


1.2 


165.6 ± 


3.2 


308.9 ± 


8.0 


553.0 ±19.3 


919.8 ± 


41.3 


> 4mJy 


20.4 ± 


0.3 


27.7 ± 


0.5 


44.0 ± 


0.9 


71.4 ± 


2.0 


116.9 ± 4.1 


190.5 ± 


8.6 


> 6mjy 


20.7 ± 


0.5 


24.4 ± 


0.7 


33.1 ± 


1.0 


51.9 ± 


2.0 


75.4 ± 3.5 


101.8 ± 


5.6 


> 8 mJy 


24.7 ± 


1.0 


25.6 ± 


1.1 


35.4 ± 


1.8 


53.5 ± 


3.3 


70.6 ± 4.9 


87.8 ± 


6.9 


> lOmJy 


32.2 ± 


2.2 


32.5 ± 


2.3 


34.2 ± 


2.4 


49.3 ± 


4.2 


68.0 ± 6.8 


85.0 ± 


9.5 



Table B23. Percentage count correction results for the uniform noise region of the 03h wide area field, over the flux density range 
2 — lOmJy, and for signal-to-noise thresholds of > 1.50(7 to > 4.00(7. 
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lOh field, uniform noise region 



Flux 
Density 






% Count Correction 








> 1.50CT 


> 2.00CT 


> 2.50CT 


> 3.00CT 


> 3.50CT 


> 4.00CT 


> 2mJy 


77.4 ± 2.6 


100.5 ± 3.9 


147.9 ± 7.0 


231.4 ±13.7 


379.4 ± 28.7 


558.0 ± 


51.2 


> 4 mjy 


58.8 ± 3.9 


63.0 ± 4.3 


73.1 ± 5.4 


88.7 ± 7.2 


104.7 ± 9.3 


129.1 ± 


12.7 


> 6mjy 


64.4 ± 7.5 


65.3 ± 7.7 


78.3 ±10.1 


90.4 ±12.5 


92.2 ±12.9 


94.0 ± 


13.3 


> 8mJy 


66.7 ±12.2 


66.7 ±12.2 


69.0 ±12.8 


76.9 ±15.1 


80.0 ± 16.0 


80.0 ± 


16.0 


> 10 mJy 


122.2 ±40.7 


122.2 ±40.7 


122.2 ±40.7 


122.2 ±40.7 


122.2 ±40.7 


122.2 ± 


40.7 



Table B24. Percentage count correction results for the uniform noise region of the lOh field, over the flux density range 2 — 10 mJy, and 
for signal-to-noise thresholds of > l.SOcr to > 4.00(T. 



14h field, uniform noise region 



Flux % Count Correction 



Density 


> 1.50 


cr 


> 2.00 


cr 


> 2.50 


cr 


> 3.00 


cr 


> 3.50cr 


> 4.00cr 


> 2 mJy 


55.5 ± 


0.6 


83.6 ± 


1.1 


146.4 ± 


2.6 


248.4 ± 


5.8 


394.8 ± 11.6 


584.9 ± 21.0 


> 4 mJy 


34.2 ± 


0.6 


42.3 ± 


0.9 


55.4 ± 


1.3 


71.9 ± 


2.0 


96.3 ± 3.1 


126.9 ± 4.6 


> 6mjy 


39.0 ± 


1.3 


50.2 ± 


1.9 


69.0 ± 


3.0 


80.2 ± 


3.7 


89.5 ± 4.4 


96.8 ± 5.0 


> 8mJy 


45.9 ± 


2.4 


52.8 ± 


3.0 


69.3 ± 


4.5 


84.3 ± 


6.0 


92.8 ± 6.9 


96.5 ± 7.3 


> 10 mJy 


61.1 ± 


5.7 


61.6 ± 


5.8 


71.1 ± 


7.2 


80.2 ± 


8.7 


87.3 ± 9.8 


90.8 ± 10.4 



Table B25. Percentage count correction results for the uniform noise region of the 14h field, over the flux density range 2 — 10 mJy, and 
for signal-to-noise thresholds of > l.SOcr to > 4.00ct. 



22h field, uniform noise region 



Flux 
Density 






% Count Correction 








> 1.50CT 


> 2.00CT 


> 2.50CT 


> 3.00CT 


> 3.50CT 


> 4.00<T 


> 2mJy 


52.4 ± 1.9 


82.6 ± 3.7 


138.7 ± 8.0 


246.5 ± 18.9 


395.3 ± 38.4 


634.8 ± 


78.1 


> 4mJy 


40.7 ± 2.9 


42.9 ± 3.2 


50.3 ± 4.0 


59.8 ± 5.2 


76.0 ± 7.4 


119.7 ± 


14.7 


> 6mjy 


56.8 ± 8.6 


56.8 ± 8.6 


56.8 ± 8.6 


61.0 ± 9.5 


67.6 ±11.1 


78.1 ± 


13.8 


> 8mJy 


115.4 ±32.0 


115.4 ±32.0 


115.4 ±32.0 


115.4 ±32.0 


115.4 ±32.0 


125.0 ± 


36.1 


> lOmJy 


100.0 ± 50.0 


100.0 ± 50.0 


100.0 ± 50.0 


100.0 ± 50.0 


100.0 ± 50.0 


100.0 ± 


50.0 



Table B26. Percentage count correction results for the uniform noise region of the 22h field, over the flux density range 2 — lOmJy, and 
for signal-to-noise thresholds of > 1.50iT to > 4.00<t. 



SSA13 deep area field, uniform noise region 



Flux 
Density 






'a, C-'ouiiL 


(.'oiixxlioii 








> 1.50cr 


> 2.00<T 


> 2.50(7 


> 3.00cr 


> 3.50<T 


> 4.00a 


> 2mJy 


74.6 ± 2.5 


79.0 ± 2.7 


89.3 ± 3.2 


102.0 ± 4.0 


125.8 ± 5.4 


154.4 ± 


7.4 


> 4mJy 


69.1 ± 5.0 


69.4 ± 5.0 


70.5 ± 5.1 


71.7 ± 5.2 


72.8 ± 5.4 


74.0 ± 


5.5 


> 6mjy 


73.3 ± 9.5 


73.3 ± 9.5 


73.3 ± 9.5 


73.3 ± 9.5 


73.3 ± 9.5 


74.6 ± 


9.7 


> 8mJy 


119.0 ±26.0 


119.0 ±26.0 


119.0 ±26.0 


119.0 ±26.0 


119.0 ±26.0 


119.0 ± 


26.0 


> 10 mJy 


100.0 ± 28.9 


100.0 ± 28.9 


100.0 ± 28.9 


100.0 ± 28.9 


100.0 ± 28.9 


100.0 ± 


28.9 



Table B27. Percentage count correction results for the uniform noise region of the SSA13 deep area field, over the flux density range 
2 — lOmJy, and for signal-to-noise thresholds of > 1.50iT to > 4.00(7. 



SSA13 wide area field, uniform noise region 



Flux % Count Correction 



Density 


> 1.50<T 


> 2.00<T 


> 2.50<T 


> 3.00CT 


> 3.50CT 


> 4.00(7 


> 2mJy 


59.7 ± 0.8 


111.6 ± 1.9 


215.8 ± 5.2 


423.2 ± 14.2 


841.6 ±39.7 


> 1000 


> 4mJy 


20.0 ± 0.3 


30.0 ± 0.6 


51.7 ± 1.3 


96.6 ± 3.3 


189.0 ± 8.9 


327.5 ± 20.4 


> 6mjy 


20.4 ± 0.5 


26.1 ± 0.7 


35.9 ± 1.2 


54.4 ± 2.2 


91.0 ± 4.8 


138.5 ± 9.1 


> 8 rnjy 


29.7 ± 1.4 


34.5 ± 1.7 


41.4 ± 2.3 


49.5 ± 3.0 


70.6 ± 5.0 


96.5 ± 8.0 


> lOmJy 


50.0 ± 4.5 


54.8 ± 5.1 


59.4 ± 5.8 


61.8 ± 6.1 


70.8 ± 7.5 


80.8 ± 9.1 



Table B28. Percentage count correction results for the uniform noise region of the SSA13 wide area field, over the flux density range 
2 — lOmJy, and for signal-to-noise thresholds of > 1.50(7 to > 4.00(7. 
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SSA17 field, uniform 


noise region 








Flux 






% Count Correction 








Density 


> 1.50CT 


> 2.00CT 


> 2.50CT 


> 3.00CT 


> 3.50CT 


> 4.00CT 


> 2mJy 

> 4 mJy 

> 6mjy 

> 8mJy 

> 10 mJy 


66.5 ± 1.3 
41.1 ± 1.4 

53.6 ± 3.2 
76.3 ± 7.9 
84.6 ± 13.5 


96.1 ± 2.3 

46.2 ± 1.6 
57.0 ± 3.5 

76.3 ± 7.9 
84.6 ± 13.5 


148.3 ± 4.5 

56.1 ± 2.2 
62.8 ± 4.1 
79.8 ± 8.5 

89.2 ± 14.7 


226.3 ± 8.4 
69.9 ± 3.0 
71.0 ± 4.9 
83.5 ± 9.1 
89.2 ± 14.7 


329.7 ±14.8 
88.7 ± 4.3 
77.0 ± 5.6 
85.5 ± 9.4 
91.7 ± 15.3 


464.7 ± 
113.9 ± 
86.0 ± 
87.7 ± 
91.7 ± 


24.7 
6.2 
6.6 
9.7 

15.3 



Table B29. Percentage count correction results for the uniform noise region of the SSA17 field, over the flux density range 2 — lOmJy, 
and for signal-to-noise thresholds of > 1.50(T to > 4.00cr. 







SSA22 field, uniform 


noise region 








Flux 






% Count Correction 








Density 


> 1.50O- 


> 2.00O- 


> 2.50O- 


> 3.00O- 


> 3.50(7 


> 4.0C 


lo- 


> 2 mJy 

> 4 mJy 

> 6mjy 

> 8mJy 

> 10 mJy 


84.3 ± 1.9 
65.6 ± 2.9 
77.3 ± 6.2 
94.9 ± 12.4 
92.0 ± 18.4 


94.7 ± 2.3 
68.6 ± 3.1 
77.3 ± 6.2 
94.9 ± 12.4 
92.0 ± 18.4 


107.4 ± 2.8 
71.2 ± 3.3 

78.8 ± 6.4 

94.9 ± 12.4 
92.0 ± 18.4 


126.1 ± 3.5 
76.1 ± 3.6 
80.4 ± 6.6 
94.9 ± 12.4 
92.0 ± 18.4 


155.9 ± 4.8 
80.3 ± 3.9 
82.1 ± 6.8 
94.9 ± 12.4 
92.0 ± 18.4 


196.2 ± 
84.5 ± 
85.0 ± 
98.2 ± 
92.0 ± 


6.8 
4.2 
7.2 
13.0 
18.4 



Table B30. Percentage count correction results for the uniform noise region of the SSA22 field, over the flux density range 2 — lOmJy, 
and for signal-to-noise thresholds of > 1.50<t to > 4.00it. 



Lockman Hole deep field, uniform noise region 



Flux 
Density 






% Count Correction 








> 1.50CT 


> 2.00CT 


> 2.50CT 


> 3.00CT 


> 3.50CT 


> 4.00CT 


> 2mJy 


88.2 ± 3.3 


98.0 ± 3.8 


113.1 ± 4.7 


129.5 ± 5.8 


155.6 ± 7.6 


203.4 ± 


11.4 


> 4mJy 


63.2 ± 4.2 


66.7 ± 4.5 


72.4 ± 5.1 


77.8 ± 5.7 


82.8 ± 6.3 


91.1 ± 


7.3 


> 6mjy 


74.0 ± 8.7 


74.0 ± 8.7 


78.3 ± 9.4 


87.1 ± 11.1 


88.5 ±11.3 


88.5 ± 


11.3 


> 8mJy 


96.6 ± 17.9 


96.6 ± 17.9 


96.6 ± 17.9 


96.6 ± 17.9 


96.6 ± 17.9 


96.6 ± 


17.9 


> lOmJy 


113.3 ±29.3 


113.3 ±29.3 


113.3 ±29.3 


113.3 ±29.3 


113.3 ±29.3 


113.3 ± 


29.3 



Table B31. Percentage count correction results for the uniform noise region of the Lockman Hole deep field, over the flux density range 
2 — lOmJy, and for signal-to-noise thresholds of > 1.50(T to > 4.00<t. 



Hubble deep field, uniform noise region 



Flux 






'a, CJouiil 


(.'orrcxLioii 








Density 


> 1.50(7 


> 2.00(7 


> 2.50(7 


> 3.00(7 


> 3.50(7 


> 4.00(7 


> 2mJy 


102.5 ± 4.9 


102.5 ± 4.9 


102.7 ± 4.9 


104.7 ± 5.0 


108.7 ± 5.3 


118.1 ± 


6.1 


> 4mJy 


80.8 ± 7.4 


80.8 ± 7.4 


80.8 ± 7.4 


80.8 ± 7.4 


80.8 ± 7.4 


80.8 ± 


7.4 


> 6mjy 


82.7 ± 11.5 


82.7 ± 11.5 


82.7 ± 11.5 


82.7 ± 11.5 


82.7 ± 11.5 


82.7 ± 


11.5 


> 8mJy 


81.0 ±17.7 


81.0 ±17.7 


81.0 ± 17.7 


81.0 ±17.7 


81.0 ± 17.7 


81.0 ± 


17.7 


> 10 mJy 


85.7 ±32.4 


85.7 ±32.4 


85.7 ±32.4 


85.7 ±32.4 


85.7 ±32.4 


85.7 ± 


32.4 



Table B32. Percentage count correction results for the uniform noise region of the Hubble deep field, over the flux density range 
2 — lOmJy, and for signal-to-noise thresholds of > 1.50(7 to > 4.00(7. 
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B3 Output versus input flux density tables 

The following tables give the percentage of sources recovered from the simulated images, for signal-to-noise thresholds of 
> 1.50(T to > 4.00(7, according to the following criteria: 

1) Fainter. The retrieved flux density was fainter than the input source with which it had been identified (Sin > Sout). 

2) Within error bars. The input flux density lay within the IcTrms error bars of the retrieved value (Sout — errout < Sin < 
Sout + errout). 

3) Boosted. The input flux density was less than the lower error boundary on the output value, but was still within a factor 

of 2 of the measured flux density (Sout/2 < Sin < Sout — errout). 

4) Spurious / confused. The fltted flux density to the peak could not be identified with a source in the input catalogue, located 
within 8 arcseconds and a factor 2 in brightness (Sin < Sout/2). 

The percentage of sources classified as (1), (2), (3) and (4) are given in Columns 3, 4, 5 and 6 respectively. Those 
fields composed of a deep pencil beam survey within a wider-area shallower survey have had these two components treated 
separately. 



ELAIS N2 field, uniform noise region 



Output vs. input flux density statistics 

S/N % noise % % within % % confused/ 

tlircsliold only fainter la errors boosted spurious 



> 1.50O- 


100.7 ± 


0.7 


0.4 ± 


0.1 


15.0 ± 


0.1 


4.1 ± 


0.1 


80.6 ± 


0.6 


> 2.00O- 


81.4 ± 


1.0 


0.6 ± 


0.1 


13.8 ± 


0.1 


7.3 ± 


0.1 


78.4 ± 


0.7 


> 2.50O- 


60.0 ± 


0.8 


0.9 ± 


0.1 


16.4 ± 


0.2 


12.6 ± 


0.2 


70.1 ± 


0.9 


> S.OOo- 


26.6 ± 


0.5 


1.6 ± 


0.1 


21.8 ± 


0.4 


19.8 ± 


0.4 


56.8 ± 


1.1 


> 3.50CT 


14.8 ± 


0.4 


2.1 ± 


0.1 


28.6 ± 


0.8 


27.5 ± 


0.7 


41.8 it 


1.1 


> 4.00<T 


0.0 ± 


0.0 


3.4 ± 


0.1 


35.4 ± 


1.3 


35.1 ± 


1.3 


26.1 ± 


0.9 



Table B33. Output versus input flux density statistics for the uniform noise region of the ELAIS N2 fleld, for signal-to-noise thresholds 
in the range > 1.50cr to > 4.00<T. 



ELAIS N2 field, non-uniform noise region 



Output vs. input flux density statistics 
S/N % noise % % within % % confused/ 

threshold only fainter 1<t errors boosted spurious 



> 1.50O- 


95.8 ± 


1.8 


0.0 ± 


0.0 


3.4 ± 


0.1 


0.6 ± 0.1 


96.0 ± 1.8 


> 2.00(7 


84.8 ± 


2.4 


0.0 ± 


0.0 


2.6 ± 


0.1 


1.2 ± 0.1 


96.2 it 2.5 


> 2.50(7 


38.5 ± 


1.7 


0.0 ± 


0.0 


4.0 ± 


0.2 


3.5 ± 0.2 


92.5 it 4.1 


> 3.00(7 


70.4 ± 


5.9 


0.0 ± 


0.0 


7.7 ± 


0.7 


7.0 ± 0.6 


85.2 it 7.2 


> 3.50(7 


0.0 ± 


0.0 


0.0 ± 


0.0 


13.9 ± 


2.3 


25.0 ± 4.2 


61.1 it 10.2 


> 4.00(7 


0.0 ± 


0.0 


0.0 ± 


0.0 


8.3 ± 


2.4 


50.0 ± 14.4 


41.7 it 12.0 



Table B34. Output versus input flux density statistics for the non-uniform noise region of the ELAIS N2 field, for signal-to-noise 
thresholds in the range > 1.50(7 to > 4.00(7. 



Lockman Hole wide area field, uniform noise region 



Output vs. input flux density statistics 



S/N 


% noise 


% 




% within 


% 


% confused/ 


threshold 


only 




fainter 


1(7 errors 


boosted 


spurious 


> 1.50(7 


89.3 ± 


0.8 


0.3 it 


0.1 


12.2 it 0.1 


2.6 it 0.1 


84.9 it 0.7 


> 2.00(7 


69.8 it 


1.1 


0.5 it 


0.1 


11.7 it 0.1 


5.3 it 0.1 


82.5 it 1.0 


> 2.50(7 


64.3 ± 


1.2 


0.7 it 


0.1 


14.8 ± 0.3 


9.7 ± 0.2 


74.7 ± 1.3 


> 3.00(7 


35.7 ± 


1.0 


1.1 it 


0.1 


20.5 ± 0.5 


15.0 it 0.4 


63.3 ± 1.7 


> 3.50(7 


0.0 ± 


0.0 


2.4 it 


0.1 


32.7 ± 1.4 


23.3 it 1.0 


41.6 ± 1.7 


> 4.00(7 


0.0 it 


0.0 


3.0 it 


0.2 


40.7 it 2.3 


32.5 it 1.9 


23.8 it 1.4 



Table B35. Output versus input flux density statistics for the uniform noise region of the Lockman Hole wide area field, for signal-to-noise 
thresholds in the range > 1.50(7 to > 4.00(7. 
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Lockman Hole 


wide area \ 


Held, non-uniform noise region 










Output vs. 


input flux density statistics 




S/N 


% noise 


% 


% within 


% 


% confused/ 


threshold 


only 




fainter 


Icr errors 


boosted 


spurious 


> 1.50cr 


88.1 ± 


1.1 


0.1 ± 0.1 


3.7 ± 0.1 


0.5 ± 0.1 


95.7 ± 1.2 


> 2.00O- 


94.3 ± 


1.5 


0.1 ± 0.1 


2.6 ± 0.1 


0.9 ± 0.1 


96.4 ± 1.6 


> 2.50(T 


82.9 ± 


1.9 


0.1 ± 0.1 


2.6 ± 0.1 


1.4 ± 0.1 


96.0 it 2.2 


> 3.00O- 


83.4 ± 


2.9 


0.1 ± 0.1 


2.3 ± 0.1 


1.9 ± 0.1 


95.7 ± 3.3 


> 3.50CT 


79.7 ± 


5.0 


0.0 ± 0.0 


3.2 ± 0.2 


3.6 ± 0.2 


93.2 ± 5.9 


> 4.00CT 


0.0 ± 


0.0 


0.0 ± 0.0 


2.2 ± 0.3 


11.1 ± 1.7 


86.7 ± 12.9 



Table B36. Output versus input flux density statistics for the non-uniform noise region of the Lockman Hole wide area field, for 
signal-to-noise thresholds in the range > 1.50cr to > 4.00cr. 



Lockman Hole deep strip field, uniform noise region 



Output vs. input flux density statistics 
S/N % noise % % within % % confused/ 

threshold only fainter la errors boosted spurious 



> 1.50CT 


83.3 ± 


1.4 


0.8 ± 


0.1 


22.7 ± 


0.4 


6.6 it 


0.1 


69.9 it 


1.1 


> 2.00(T 


59.8 ± 


1.7 


1.2 ± 


0.1 


18.8 ± 


0.4 


10.5 it 


0.2 


69.6 it 


1.4 


> 2.50(T 


44.5 ± 


1.2 


1.7 ± 


0.1 


19.9 ± 


0.5 


16.1 it 


0.4 


62.3 it 


1.7 


> 3.00(7 


26.9 ± 


1.0 


2.6 ± 


0.1 


25.2 ± 


0.9 


23.3 it 


0.9 


49.0 it 


1.8 


> 3.50CT 


23.1 ± 


1.1 


3.9 ± 


0.2 


30.7 it 


1.5 


29.8 it 


1.4 


35.6 it 


1.7 


> 4.00CT 


0.0 ± 


0.0 


5.2 ± 


0.3 


36.9 it 


2.3 


35.3 it 


2.2 


22.6 it 


1.4 



Table B37. Output versus input flux density statistics for the uniform noise region of the Lockman Hole deep strip field, for signal-to- 
noise thresholds in the range > 1.50rr to > LOOfr. 



03h deep area field, uniform noise region 



Output vs. input flux density statistics 
S/N % noise % % within % % confused/ 

threshold only fainter la errors boosted spurious 



> 1.50O- 


62.0 ± 


1.7 


3.1 ± 


0.1 


35.9 ± 


1.0 


15.9 it 


0.4 


45.1 ± 


1.3 


> 2.00(7 


61.5 ± 


2.0 


3.4 it 


0.1 


30.5 ± 


1.0 


21.0 it 


0.7 


45.1 ± 


1.4 


> 2.50(7 


41.6 ± 


1.5 


4.0 it 


0.1 


29.1 ± 


1.1 


26.8 it 


1.0 


40.1 ± 


1.5 


> 3.00(7 


0.0 ± 


0.0 


5.1 it 


0.2 


30.9 ± 


1.3 


30.9 ± 


1.3 


33.0 ± 


1.4 


> 3.50(7 


0.0 ± 


0.0 


5.8 it 


0.3 


31.8 ± 


1.6 


35.8 ± 


1.8 


26.5 ± 


1.4 


> 4.00(7 


0.0 it 


0.0 


7.5 it 


0.5 


34.0 it 


2.1 


38.5 it 


2.4 


20.0 it 


1.2 



Table B38. Output versus input flux density statistics for the uniform noise region of the 03h deep area field, for signal-to-noise 

thresholds in the range > 1.50(7 to > 4.00(7. 



03h wide area field, uniform noise region 



Output vs. input flux density statistics 
S/N % noise % % within % % confused/ 

threshold only fainter la errors boosted spurious 



> 1.50(7 


89.4 it 


1.0 


0.7 it 


0.1 


19.5 ± 


0.2 


5.2 ± 


0.1 


74.5 ± 


0.8 


> 2.00(7 


74.2 ± 


1.2 


1.0 it 


0.1 


18.2 ± 


0.3 


9.0 ± 


0.1 


71.7 ± 


1.0 


> 2.50(7 


54.4 it 


1.0 


1.6 it 


0.1 


20.7 it 


0.4 


14.5 it 


0.3 


63.2 it 


1.2 


> 3.00(7 


27.1 it 


0.7 


2.4 it 


0.1 


26.3 it 


0.7 


19.8 it 


0.5 


51.5 it 


1.3 


> 3.50(7 


12.1 ± 


0.4 


4.1 it 


0.1 


30.4 it 


1.1 


25.3 it 


0.9 


40.1 it 


1.4 


> 4.00(7 


20.2 it 


0.9 


5.4 it 


0.2 


34.7 it 


1.6 


30.0 it 


1.3 


29.8 it 


1.3 



Table B39. Output versus input flux density statistics for the uniform noise region of the 03h wide area field, for signal-to-noise 
thresholds in the range > 1.50(7 to > 4.00(7. 
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lOh field, uniform noise region 



Output vs. input fiux density statistics 
S/N % noise % % within % % confused/ 

threshold only fainter la errors boosted spurious 



> 1.500- 


58.7 ± 


1.8 


1.8 ± 


0.1 


35.5 ± 


1.1 


10.5 ± 


0.3 


52.3 ± 


1.6 


> 2.00O- 


59.4 ± 


2.3 


2.1 ± 


0.1 


31.1 ± 


1.2 


15.9 ± 


0.6 


51.0 ± 


2.0 


> 2.50(T 


66.8 ± 


3.2 


2.4 ± 


0.1 


31.4 ± 


1.5 


21.4 ± 


1.0 


44.8 ± 


2.1 


> 3.00O- 


0.0 ± 


0.0 


3.8 ± 


0.2 


32.1 ± 


1.9 


26.8 ± 


1.6 


37.3 ± 


2.2 


> 3.50CT 


0.0 ± 


0.0 


5.7 ± 


0.4 


34.3 ± 


2.6 


29.1 ± 


2.2 


30.9 ± 


2.3 


> 4.00O- 


0.0 ± 


0.0 


6.7 ± 


0.6 


36.1 ± 


3.3 


32.8 ± 


3.0 


24.4 ± 


2.2 



Table B40. Output versus input flux density statistics for the uniform noise region of the lOh field, for signal-to-noise thresholds in the 
range > 1.50ct to > 4.00ct. 



14h field, uniform noise region 



Output vs. input flux density statistics 
S/N % noise % % within % % confused/ 

threshold only fainter la errors boosted spurious 



> l.SOo- 


88.6 ± 


0.9 


1.3 ± 


0.1 


28.5 ± 


0.3 


7.9 ± 


0.1 


62.4 ± 


0.7 


> 2.00(T 


57.9 ± 


1.2 


1.7 ± 


0.1 


27.7 ± 


0.4 


12.9 ± 


0.2 


57.8 ± 


0.8 


> 2.50(T 


28.9 ± 


0.5 


2.5 ± 


0.1 


30.3 ± 


0.5 


19.8 ± 


0.4 


47.3 ± 


0.8 


> 3.00(7 


10.9 ± 


0.3 


3.4 ± 


0.1 


33.7 ± 


0.8 


26.6 ± 


0.6 


36.3 ± 


0.8 


> 3.50O- 


0.0 ± 


0.0 


4.7 ± 


0.1 


38.3 ± 


1.1 


29.8 ± 


0.9 


27.2 ± 


0.8 


> 4.00CT 


0.0 ± 


0.0 


5.6 ± 


0.2 


40.7 ± 


1.5 


33.4 ± 


1.2 


20.3 ± 


0.7 



Table B41. Output versus input flux density statistics for the uniform noise region of the 14h field, for signal-to-noise thresholds in the 
range > 1.50ct to > 4.OO0-. 



22h field, uniform noise region 



Output vs. input flux density statistics 
S/N % noise % % within % % confused/ 

threshold only fainter la errors boosted spurious 



> 1.50(7 


94.8 ± 


3.3 


1.3 ± 


0.1 


28.9 ± 


1.0 


7.8 it 


0.3 


62.0 ± 


2.1 


> 2.00(7 


39.4 ± 


4.2 


1.6 ± 


0.1 


26.6 it 


1.2 


13.0 it 


0.6 


58.8 ± 


2.6 


> 2.50(7 


33.1 ± 


1.9 


2.6 ± 


0.2 


27.8 ± 


1. 


19.2 ± 


1.1 


•50.. 3 ± 


2.9 


> 3.00(7 


0.0 ± 


0.0 


3.5 ± 


0.3 


34.1 ± 


2.6 


22.9 ± 


1.8 


39.4 ± 


3.0 


> 3.50(7 


0.0 ± 


0.0 


5.7± 


0.5 


40.6 ± 


3.9 


23.6 ± 


2.3 


30.2 ± 


2.9 


> 4.00(7 


0.0 it 


0.0 


9.1 ± 


1.1 


45.5 it 


5.6 


25.8 it 


3.2 


19.7 it 


2.4 



Table B42. Output versus input flux density statistics for the uniform noise region of the 22h field, for signal-to- noise thresholds in the 
range > 1.50(7 to > 4.00(7. 



SSA13 deep area field, uniform noise region 



Output vs. input fiux density statistics 
S/N % noise % % within % % confused/ 

threshold only fainter la errors boosted spurious 



> 1.50(7 


122.0 ± 


3.2 


4.3 it 


0.1 


30.6 ± 


0.8 


18.0 it 


0.5 


47.1 ± 


1.2 


> 2.00(7 


40.8 lb 


3.5 


4.2 it 


0.1 


27.6 ± 


0.8 


21.6 it 


0.6 


46.6 it 


1.3 


> 2.50(7 


42.4 it 


1.4 


4.4 it 


0.1 


24.6 it 


0.8 


26.9 it 


0.9 


44.1 it 


1.4 


> 3.00(7 


0.0 it 


0.0 


5.3 it 


0.2 


24.4 it 


0.9 


29.6 it 


1.1 


40.7 it 


1.5 


> 3.50(7 


0.0 it 


0.0 


5.7 it 


0.2 


24.2 it 


1.0 


33.9 ± 


1.5 


36.3 it 


1.6 


> 4.00(7 


0.0 it 


0.0 


6.4 it 


0.3 


23.5 it 


1.1 


37.1 it 


1.8 


33.0 it 


1.6 



Table B43. Output versus input fiux density statistics for the uniform noise region of the SSA13 deep area field, for signal-to-noise 
thresholds in the range > 1.50(7 to > 4.00(7. 
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SSA13 wide area field, uniform noise region 



Output vs. input flux density statisties 
S/N % noise % % within % % confused/ 

threshold only fainter Ict errors boosted spurious 



> 1.500- 


96.1 ± 


1.2 


0.7± 


0.1 


16.7 it 


0.2 


3.9 it 


0.1 


78.7 ± 


1.0 


> 2.00O- 


67.9 ± 


1.7 


1.1 ± 


0.1 


15.7 ± 


0.3 


7.3 it 


0.1 


75.9 ± 


1.3 


> 2.50(T 


51.4 ± 


1.2 


1.7 ± 


0.1 


18.8 it 


0.5 


12.2 it 


0.3 


67.3 ± 


1.6 


> 3.00O- 


11.2 ± 


0.4 


2.1 ± 


0.1 


24.2 it 


0.8 


18.7 it 


0.6 


55.0 it 


1.8 


> 3.50CT 


22.3 ± 


1.1 


3.1 ± 


0.1 


31.4 it 


1.5 


26.7 it 


1.3 


38.8 it 


1.8 


> 4.00O- 


0.0 ± 


0.0 


4.7 ± 


0.3 


36.4 it 


2.3 


32.6 it 


2.0 


26.4 it 


1.6 



Table B44. Output versus input flux density statistics for the uniform noise region of the SSA13 wide area field, for signal-to-noise 
thresholds in the range > 1.50cr to > 4.00(7. 



SSA17 field, uniform noise region 



Output vs. input flux density statistics 
S/N % noise % % within % % confused/ 

threshold only fainter la errors boosted spurious 



> 1.50CT 


91.7 ± 


1.8 


1.7 ±0.1 


29.7 it 


0.6 


10.2 it 


0.2 


58.5 it 


1.1 


> 2.00CT 


51.9 it 


2.2 


2.0 it 


0.1 


25.6 it 


0.6 


16.1 it 


0.4 


56.3 it 


1.4 


> 2.50CT 


9.0 it 


0.3 


2.5 it 


0.1 


26.8 it 


0.8 


22.0 it 


0.7 


48.7 it 


1.5 


> 3.00O- 


13.8 it 


0.5 


3.4 it 


0.1 


29.4 it 


1.1 


27.9 it 


1.0 


39.2 it 


1.5 


> 3.50O- 


20.0 it 


0.9 


4.4 it 


0.2 


30.3 it 


1.4 


31.7 it 


1.4 


33.7 it 


1.5 


> 4.00CT 


0.0 it 


0.0 


5.6 it 


0.3 


32.2 it 


1.7 


35.3 it 


1.9 


26.8 it 


1.4 



Table B45. Output versus input flux density statistics for the uniform noise region of the SSA17 field, for signal-to-noise thresholds in 
the range > 1.50cr to > 4.00rr. 



SSA22 field, uniform noise region 



Output vs. input flux density statistics 
S/N % noise % % within % % confused/ 

threshold only fainter la errors boosted spurious 



> 1.50(7 


89.1 ± 


1.6 


4.8 it 


0.1 


35.1 ± 


0.6 


16.9 ± 


0.3 


43.3 ± 


0.8 


> 2.00a 


62.7 ± 


1.8 


5.4 it 


0.1 


31.7 ± 


0.6 


22.2 ± 


0.5 


40.8 ± 


0.8 


> 2.50a 


11. 0± 


0.3 


5.7it 


0.1 


28.5 ± 


0.7 


28.0 ± 


0.7 


37.7 ± 


0.9 


> 3.00a 


0.0 ± 


0.0 


5.9 it 


0.2 


28.0 ± 


0.8 


32.4 ± 


0.9 


33.7 ± 


0.9 


> 3.50a 


0.0 ± 


0.0 


7.2 ± 


0.2 


26.2 ± 


0.8 


36.6 ± 


1.1 


30.0 ± 


0.9 


> 4.00CT 


0.0 it 


0.0 


7.9 it 


0.3 


26.9 it 


0.9 


39.6 it 


1.4 


25.7 it 


0.9 



Table B46. Output versus input flux density statistics for the uniform noise region of the SSA22 field, for signal-to-noise thresholds in 
the range > 1.50it to > 4.00(7. 



Lockman Hole deep field, uniform noise region 



Output vs. input flux density statistics 
S/N % noise % % within % % confused/ 

threshold only fainter la errors boosted spurious 



> 1.50(7 


52.2 ± 


1.5 


5.1 it 


0.2 


36.3 ± 


1.1 


16.7 it 


0.5 


41.9 ± 


1.2 


> 2.00(7 


65.3 it 


1.7 


5.5 it 


0.2 


33.6 ± 


1.1 


20.9 it 


0.7 


39.9 ± 


1.3 


> 2.50(7 


58.1 it 


2.2 


5.8 it 


0.2 


30.8 it 


1.2 


26.9 it 


1.0 


36.6 it 


1.4 


> 3.00(7 


37.5 it 


1.6 


6.2 it 


0.3 


31.1 it 


1.3 


30.5 it 


1.3 


32.2 it 


1.4 


> 3.50CT 


0.0 it 


0.0 


6.5 it 


0.3 


29.4 it 


1.4 


33.5 it 


1.6 


30.6 it 


1.5 


> 4.00<7 


0.0 it 


0.0 


7.8 it 


0.4 


30.1 it 


1.7 


33.9 it 


1.9 


28.2 it 


1.6 



Table B47. Output versus input flux density statistics for the uniform noise region of the Lockman Hole deep field, for signal-to-noise 
thresholds in the range > 1.50(7 to > 4.00(7. 
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Hubble deep field, uniform noise region 



Output vs. input fiux density statistics 
S/N % noise % % within % % confused/ 

threshold only fainter la errors boosted spurious 



> 1.500- 


81.5 ± 


2.8 


9.8 ± 


0.3 


45.6 ± 


1.6 


26.7 ± 


0.9 


17.9 ± 


0.6 


> 2.00O- 


53.8 ± 


3.0 


9.3 ± 


0.3 


41.9 ± 


1.5 


30.8 ± 


1.1 


18.0 ± 


0.7 


> 2.50(T 


31.1 it 


1.2 


9.9 ± 


0.4 


38.5 ± 


1.5 


34.5 ± 


1.4 


17.1 it 


0.7 


> 3.00O- 


18.6 ± 


0.8 


9.5 ± 


0.4 


35.4 ± 


1.5 


38.0 ± 


1.6 


17.1 it 


0.7 


> 3.50CT 


21. 7± 


1.0 


9.6 ± 


0.4 


33.3 ± 


1.6 


40.9 ± 


1.9 


16.3 it 


0.8 


> 4.00O- 


0.0 ± 


0.0 


10.7 ± 


0.5 


29.4 ± 


1.5 


43.1 ± 


2.2 


16.8 it 


0.8 



Table B48. Output versus input flux density statistics for the uniform noise region of the Hubble deep field field, for signal-to-noise 
thresholds in the range > l.bOcr to > 4.00iT. 
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